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The  purpose  of  this  project  was  to  review  the  state  of  the  art  on  the 
dynamic  response  of  composites  and  to  suggest  possible  future  research 
directions.  More  than  three  hundred  papers  are  reviewed  and  compiled  at 
the  end  of  this  report.  The  existing  theories  of  the  dynamic  response  of 
composites  are  summarized  in  Chapter  I.  The  harmonic  waves  and  transient 
waves  in,  and  the  experiments  on,  composites  are  discussed  in  Chapters  II, 
III  ard  IV,  respectively.  The  differences  between  various  theories  and 
possible  future  research  directions  are  presented  in  Chapter  V. 
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A Ri;vii;w  or 

dynamic:  Riisi’ONsi;  or  coMpasirrs 


lNTR01)Uc:riON 

Tho  JyiKunic  rosponso  of  a oomposito  has  occn  a liotly  pursued 
research  subject  since  the  middle  of  nineteen  sixties.  It  is  a subject 
which  has  wide  technological  applications  and  also  possesses  dial  len.e  ini; 
theoretical  and  experimental  problems.  Vast  literatures  are  now  avail- 
.ible  on  the  subject.  Various  theories  have  been  proposed  to  prcxiict  the 
dynamic  response  of  a composite.  I'he  ultimate  goal  of  research  in  tliis 
area  is  to  obtain  an  approximate  theory  which  is  reasonably  simple  and, 
at  the  same  time,  is  able  to  predict  fairly  accurately  the  res|H>nse  of  a 
composite  stnicture  subject  to  a dni;unic  loading.  Despite  the  voluminous 
research  papers  published,  this  goal  does  not  seem  to  have  been  achieved. 
Ibis  is  not  a reflection  on  the  lack  of  research  ability  in  the  area, 
l^iite  the  contrary,  there  are  several  .sophisticated  theories  which  are 
able  to  predict  accurately  certain  aspects  of  the  dyn;imic  response  of  a 
composite.  This  is  a reflection  on  the  difficulty  of  analyzing  a com- 
posite material.  One  can  have  an  exact  or  nearly  exact  tlteory  which  is 
either  mathematically  intractable  or  practically  unfeasible.  On  the  other 
hand,  one  can  have  a very  simple  approximate  theory  which  is  too  crude  to 
predict  even  the  simplest  dynamic  resjxnise  of  a comiiosite. 

Hte  purpose  of  this  project  was  to  critically  review  the  state  of 
the  art  on  the  subject  of  the  d\ai;imic  resjxntse  of  cixnposites  and  to  suggest 
possible  future  research  directions.  As  we  embarktxl  on  the  project,  it 


soon  became  clear  that  the  task  was  a much  more  difficult  one  than  we 


i 


liave  anticipated.  Tlie  more  tiian  three  liundred  pajiers  compiled  at  tl\e  end  [ 

( 

ol'  tills  report  hy  no  means  exhaust  all  papers  in  the  area.  Papers  wliich  j’ 

deal  with  the  static  response  of  a composite  are  not  included  in  the 

references.  A review  of  the  more  than  tiiree  liundred  papers  reveals  that  * 

most  papers  discuss  either  haniioiiic  waves  or  transient  waves.  Of  course, 

there  are  papers  which  discuss  lioth  harmonic  and  transient  waves.  When  j 

the  material  is  linear,  harmonic  waves  can  he  superimposed  to  obtain  a 
transient  wave.  This  approach  is  not  applicable  for  nonlinear  materials. 

Theri'fitre , papers  which  diseuss  harnu'iiic  waves  invariahli'  as.sume  that  the 
ma  t er  I a 1 is  1 i nea  r . 

Pul'li-.lu'd  paper-,  on  plate-,  .nul  slu-ll-.  which  are  made  of  eomposite 
materials  were  also  rcvieweil.  Although  the  review  of  these  papers  is  not 
presented  here,  the  papers  are  included  in  the  Kefercrccs  at  the  end  of 
this  report.  I.ikewise,  review  of  papers  which  model  comiuisitcs  as  fluids 
IS  not  presented  hut  tlie\  .-irc'  itiv  liided  in  the  Keferences. 

In  I'hapter  1.  we  hrieflv  review  the  existitut  theories  ol'  the  dynamic 
response  of  a composite.  This  is  fi>llowed  by  a review  in  Chapter  11  of 
papers  which  deal  with  li.innonic  waves  in  composites.  Tr.-insient  waves  in 
composites  are  discussed  in  Chapter  111.  Althoiu'.h  the  orijtinal  objectives 
of  this  jiroject  do  not  include  reviewinji  the  experimental  results,  we  felt 
that  some  ex]H'rinietital  results  are  of  sufficient  interest  and  are  relevant 
to  the  theoretical  predictions  that  a few  words  should  be  said  about  them. 

This  is  contained  in  Chaiiter  IV.  I'inally,  in  Chapter  V we  comment  on  the 
differences  betwv'cn  various  theories  and  also  sui^j’est  possible  future 


rc-.e.i  rch  d i iw'c  t i ons  . 


I . A SIMHARY  Ot-  jm  t:XKSTINC7  THliORlt-S 

For  most  realistic  structural  composites,  an  exact  description  of  the 
static  or  dynamic  behavior  is  mathematically  impracticable.  As  an  alterna- 
tive, a number  of  investigators  have  sought  approximate  theories.  The 
representatives  of  such  theories  are  briefly  described  as  follows; 


l.l  tffective  Modulus  Theories 


The  effective  modulus  tlieorics  such  as  tliosc  proposed  by  Postma  [ Pb  J 
and  IVliite  and  Angona  [Will  replace  the  actual  composite  by  a homogeneous, 
generally  anisotropic  medium  whose  material  constants  arc  a geometrically 
weighted  average  of  the  properties  of  the  constituents.  If  '""c  such 

effective  iiunluli  of  the  composite,  this  theory  relates  volume  averages  of 
stresses  to  volume  averages  of  strains  by  a general  anisotropic  linear 


stress-strain  relation  of  the  form 


T . . = C . . . 0 C,  „ , 
ij  l.tkil  kt 


i , j , k , = 1,2,,^ 


Of  course,  the  constants  0*  , „ are  expressions  in  tenns  of  the  material 

ijkt  ‘ 

constants  of  the  constituents  and  the  parameters  defining  the  geometric 

layout  of  the  composite.  The  constants  C?.,.  satisfy  the  relations 

1 jkx. 

I’*  = r*  = ('*  = r*  t 

^ijk£  jik£  ijik  k£ij 

Thus,  of  the  81  constants  only  21  arc  independent.  In  general,  tlie 

number  of  independent  elastic  constants  is  much  less  than  21  because  of  the 
existence  of  sytranctries  in  the  structuring  of  the  material.  In  particular, 
in  the  case  of  a laminateii  medium  consisting  of  alternating  layers  of  two 
isotropic  elastic  materials,  the  number  of  independent  elastic  constants 
reduces  to  only  5. 


While  yioUlin,'.  sat  ist'aetory  results  tor  certain  itoomeiries  under 
static  loads,  the  effective  modulus  tl\eories  exhibit  serious  deficiencies 
for  virtually  all  j^eometries  when  applied  to  wave  propajtat  ion.  Sped f ica  1 1\ , 
these  theories  are  incapable  of  reproducinj’  the  dispersion  and  attenuation 
observed  in  composites.  .Such  effects  become  important  where  dominant 
sii^nal  wave  lenjtths  are  of  tlie  oi\ler  of  the  typical  ciunposite  micrixlimen- 
sion.  Since  dispersion  .ind  attenuation  are  results  of  the  microstnicturo 
with  discontinuous  material  properties,  any  continuum  theory  must  in  one 
wa>-  or  anotlier  take  into  account  tlie  influence  of  microstructure.  five 
followin.c  theories  were  develi'pcd  with  this  purpose  in  mind. 

l.J  d'fectue  Stiffness  I'heories 

Vhe  effective  stiffness  theory  was  tlie  first  continuum  model  for 
l.iminated  media  and  fiber-reinforced  composites  to  account  for  a typicalli- 
diiKunic  effect  such  as  >;eometric  dispersion  and  hence  to  reflect  the  in- 
fluence of  the  microstructuro  of  a Cv'miiosite.  fhe  theories  were  develo\H'd 
in  a series  of  papers  [ S IS.A'.)  .il  It),  tl  1 1 . ,M0|  In  Achenbach,  tirot  , Herrmann 
■ind  Sun.  r.ie.hor  order  theories  of  this  kind  were  derived  by  furhan  I'fl.S). 

the  theories  h.ii  e been  I'l'niml.ited  in  sever. il  dil'ferent  foniis.  but 
the  c.ise  of  the  line.irlx  el.ist  ic  l.iminated  composites  is  I'erh.ips  ,i  tvincal 
one.  Here  we  outline  the  theor'  ot'  el.ist  ic  w.iies  m l.iminated  composites 
1SIS|  briel'l'.  fhe  re  i nt'ore  i iie,  .and  ni.itrix  l.iiers  .ire  hot  l\  assumed  t be 
homojieneous , linear  isotropic  elastic  materials,  lor  elastic  waves  propa 
,C.it  iiu;  in  tile  composites,  this  theorv  app  ro\  i m.i  t es  the  vl  i sp  1 .m-iiK'ii  t s ,if  th>' 
reint'orced  la>er  .iiul  the  matrix  layer  in  the  kth  cell  as 
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where  and  denote  the  displacements  of  the  reinforcing  and  matrix 

layers  respectively,  and  denote  the  displacements  at  the  mid- 
planes of  the  two  corresponding  layers,  and  x™  local  coordinates 

fk  fk 

measured  from  the  corresponding  midplanes,  and  represent  anti- 

f k 

symimetric  thickness  shear  deformations,  and  ^22  represents  symmetric 
stretch  deformation  of  the  kth  reinforcing  layer.  Similar  definitions 
apply  for  >1^23  '^22'  approximate  theory  allows  dynamic 

interaction  of  the  layers  through  the  continuity  of  displacements  at  the 
interfaces.  This  is  obtained  from  Eiq.  (3)  as 


mk,  mk  fk,  fk  . 

u^iCx^,x,  .x,.t)  - u^.CXj.x^  ,X3.t) 

1 , ,fk,  fk  1 , ,mk,  mk 

= y Jf  '^2i*-^l’^2  I ‘^tn  '^2i*-^l’^2 


(4) 


With  the  assumption  of  the  displacement  fields  given  by  Eq.  (3),  one  can 

f Ic  jnk 

obtain  the  strain,  and  consequently  the  elastic  strain  energy  W and  w 
in  the  kth  reinforcing  layer  and  the  matrix  layer,  respectively.  One  also 
obtains  the  kinetic  energy  1^'^  and  Now,  if  the  composite  consists 

' of  n reinforcing  layers  and  n matrix  layers  within  a certain  thickness  i, 
the  total  strain  energy  and  kinetic  energy  are 

= I (5) 

Tj  = I (T^''  + C6) 

The  basic  assumption  in  the  effective  stiffness  theory  is  the  smooth  opera- 

^ JTl 

tion  in  which  W.  is  expressed  in  terms  of  W and  Vv  by 


‘ (W^  + w“')  dx. 


(7) 


where  and  W™  are  now  defined  for  all  X2.  If  the  layering  thicknesses 


(> 


are  small,  W*  anJ  w”*  arc  approximately  eiiual  to  and  within 

each  layer.  rherefore,  the  strain  energy  density  W can  be  defined  as 

W = iw‘^  . l\'"'S/ldj.  * d^J  (.8) 

where  and  are  assumed  to  hold  for  all  \^,  not  just  x'^*^  and 

Simil;»r  smootliinj;  operation  is  applied  to  the  kinetic  enei'Ky.  By 
assumins  the  smallitess  in  the  layer  thicknesses,  hi).  (4)  can  he  ajiproximated 
by  a differential  form; 


S.  (.d,-  ' ‘1  I'"*  ,1'  lx,  .X  ,.x,  .t  1 
I 1 m J 1 I J 


I in 


I .•  j'  I ( X , , X , , X . . t ) d li'"!  ■ I X , X , , X . , t ) U 

t Jl  I J Ill  1 J .1 


(•'inalli-,  one  invokes  Hamilton's  ['rincii'le  in  wliich  t lu'  continuity  conditions 
jtn  .ire  included  b>’  usin^  the  l.a;,;iMni)  i.in  iiui  1 1 i ]'l  i ers  \.; 

iS  [ (.T  - W - - \ ,S,  - A.S,)dV  dt  = 0 (.10) 

♦ 


Ibis  results  in  a system  of  |iartial  differential  ei)uat  ions  for  the  dis)il;ice- 
ments  u . and  'll,.. 

oi  _i 

Sun,  Achenbach  and  Herrmann  ISl.S]  tlien  used  these  disj^lacement  equations 
to  study  tile  jiropai^.it  i on  of  |ilane  harmonic  waves  in  a huiiinated  medium.  His- 
persion  relations  for  liarmonic  waves  )iropai;at  injt  parallel  to  and  normal  to 
the  direction  of  the  layerini^  were  jiresented,  .ind  tlie  ajiproximate  disjiersion 
curves  were  conqiared  with  exact  curves.  Ihe  limit  in>t  phase  velocities  at 
vanishing  wave  members  ayjree  with  the  exact  limits.  The  lowest  ant  i siiiimet  r ic 
moile  for  waves  iiroi>ai;.it  i iii;  in  the  direction  of  the  layeriii)!  sliows  the 
strongest  dispersion  wliich  is  veri’  adequately  described  by  this  theory  over 


— ' n 

I 

7 

a siibstaiitial  raiij^c  of  wave  numbers.  Various  theories  of  iMfoctive  stiff- 
uoss  will  1)0  roviowod  la  tor. 

( 

1..^  Mixtm  o Thooi-  i os 

Aiiothor  approach,  sii>;j;ost  od  hy  l.omprioro  11.7),  is  to  uso  tho  mixturo 
tlioorios  as  modols  of  tlio  dynamics  of  composites.  I'ho  t'liiKlamont a 1 conco|)t 
of  mixturo  was  postulated  liy  Truosdoll  and  Toujiin  |TIJ),  and  furtlior 
developed  hy  I'.reen  and  Nay,hdi  1 1 ll? , I'.S  | , fireen  and  Steel  ll'>(').  Steel  ISM), 
and  others.  In  these  theories,  the  constituents  of  the  structui'al  composite 
are  superimposed  in  space  and  allowed  to  underj’o  individual  deformations. 

The  microstructure  of  the  composite  is  then  simulated  hy  specifying  the 
nature  of  constituent  interactions  and  the  form  of  the  mixture  constitutive 
re lat ions. 

While  general  conservation  taws  j;t)vernin>;  the  mixture  may  easily  he 
tormulated,  the  practical  application  to  ciiinposite  materials  encounters 
difficulties  in  that  it  is  rather  difficult  to  analytically  specify  the 
interactions  between  the  const  i tutent  s on  the  basis  ot"  the  knowledjte  of 
the  j;eometry  and  constitutive  relations  of  the  individual  constituents. 

In  11)71,  hedford  and  Stern  lllll)  first  proposed  a mixture  theory  for  a 
laminated  composite  wherein  the  interaction  par;uneters  were  determined  on 
the  basis  of  results  of  certain  simple  i(uas i -st at  ic  problems,  llien,  in  a 
series  of  papers  hy  Bedford  and  Stern  [SlS.Slo.lU.'^l  , llenemier  and  Nayfeh 
1117),  and  llenemier,  Ciurtman  and  Nayfeh  IMS),  mixture  theories  were  formu- 
lated for  certain  huninated  and  fiber-reinforced  composites  with  varying 
degrees  ot  success.  In  the  tol lowing,  we  tnitline  the  binary  mixture  theory 
for  wave  guide-type  propagation  in  laminated  and  unidirectional  fibrous 
comi)osites  formulated  by  llegemier,  Curtman  and  Nayfeh  in  |H8|. 
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l\'r  .1  juTii'ilu'  .u'lMv  ol.istic,  i ‘iot  ri'pio  .tiul  t»i>mo,i;t‘iu'ou-. 

Im  I am  iM.it  o-i , I>ohkU\I  at  titoii'  i i»t  I'lt  a^v". , tln'-.i'  authors  first  i nt  t'.i;ra  I isl 
t ho  o I as  t »\l\ iiaiii  10  o<|uaiu>ns  ol  im'lion  aiul  ooiist  1 1 ut  iv  o rolatioii-.  lor  isioh 
oi'ii-.t  1 1 uoiil  o\or  tlio  tiuoknoss  I'f  o.ioh  ooiistitnonl  aiul  ilofiiu'il  avoraj;tsl 
strossos  aiul  il  i sp  I aoi'mont  s ovi'r  that  thiokitoss.  By  usinr.  the  oi'iulitioii 
that  o^^  mtist  ho  v’oiit  i iiiioits  aori'ss  i iit  or  t’aoi's  , thov  ohtaiiioil  t lu'  momoiiium 
oipiat  ions  in  tt\o  I'orm 


, lip'  ll''  , . I lal 
.1  i'  ' p ,>-u  r 


I t ^ 


\ \\  j t - 

in  uhioh  I ho  supi' rso  r i pt  s or  snhsoripts  I anJ  rotor  to  I aiul  roust  i t iiont  s 
rospoo  I I V i' 1 \ . t lio  siipiM’so  r I pt  a rofiM's  to  ati'rayp'  \aliK',  aiul  "partial" 


.1  lu'ssos  aiul  vloiisitios  ai'v'  ilol’iiuul  as 

lopl  IvX.a)  Ipl 

o ' n o . t)  - II  (' 

\\  aw  u o a 


uluTi" 


n - h (h,  ' hj 

>\  a 1 


a 1 . J 


IS  a \v'luiiio  t'raotioii  v't  t ho  roust  1 1 iiont  , .iiul  li  is  oiu'  hall  of  I lu‘  lliuk 

v\ 

iioss  of  t ho  a r v'lisi  1 1 iii'iil  . In  1.111,  1’  is  .in  " i iit  oisu  t ion"  toriii  rot  1 or  i i iii', 

nu'iiuMitum  traiistor  I'roiii  i>iio  r i>ns  t i t iionl  to  anotIuM-  vi.i  slusir  i nt  or.ir  t i on 
.irri'ss  laiiiiiuito  i nt  orl'ar  os  . U\  a r.itional  analNsis,  Uoi'.oiiii  or , otr.,  toiiiul 
till'  i nt  orar  I I v'li  toriii  I'  takinr.  tho  tonii 


k /^ilal  ^,iJal\ 

< h , 1 • \ V ^ / 


uhoro 
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'|i||i  . IH|  n , • p .H|  1 


Similarly,  intojii'^it  ion  of  tho  consitutivo  relations  for  the  individual  con- 
stituents followed  hy  a rational  analysis,  He^cinier,  etc.,  found  the 
constitutive  relations  for  the  mixture  as  follows: 


Up)  U;‘) 

' = c,  , 3 lU  ' + c,  -3 

XX  11  X X 12  X X 


(-p)  1 (1‘0  !“•'') 

‘ = c,  ,3  lU  ♦ c , ,3  lU 

XX  12  X X 22  X X 


where 


/ i,  a J a 1^ 

C = I n 1-.  - — I , C „ = — ;r-- 

au  \ a a f.  / al-i  I. 


(.a . B = 1,2;  u / ti) 


in  which 


■ 1040,, 


;ind  \ , p arc  Lame  constants  of  tlie  a-const  i tuent. 

a a 


for  fibrous  composites,  llegemier,  etc.,  approximated  a hexagonal  array 
of  fibers  by  coi\centric,  linear  elastic  cylinders,  witli  perfect  interface 
bonds  and  subject  to  vanishing  shear  stress  and  r.adi.il  displacement  on  the 
outer  bouiidaries  so  that  for  a cylindrical  element,  a = 1 (r  < r^ ) ilenotes 
fiber  and  a = 2 (,r  ^ < r ^ r ,)  denotes  matrix.  fhey  found  that  tlie 
momentum  enuations  for  fibrous  composites  can  be  also  written  in  the  fom 
Ul)  where,  for  this  case. 


SOjUjU, 


X X 


in  wi\ich 


, . (i  . i) 


r‘,  - r-. 
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I'hoy  ;ilso  t'ouiul  that  t hi'  const  i t lit  i vi'  relations  of  the  mixture  still  take 
the  form  (l(il  wliere 


c = (\  + -’ll' In  - XVI).  c-  , , = X X,,/l).  (u/(X) 

ixa  It  («  It  It  uH  It  n 

0 - ' I I ' , • e , 1 n , • ( V ' Ii  , in  Ml  .fn  , 1 

n,|  M , 1 1 - . - 1 - 1 


VI  I 


(.'om|ia  r ! son  ol  exact  anO  app  rox  i ma  t e phase  velocity  ilata  for  lanunatoil  aiul 
fibrous  composites  imlicate  that  the  theories  just  ilescriheil  provide  ^oikI 
ajtreement  for  wavelenitths  p.reater  than  the  typical  composite  m i c roil  imens  ion  . 


1.1  fontitnmm  I'heories  hased  mi_^sympt  ot  ic  hxpansions 

for  laminates  and  d i rec t i ona 1 1 >■  reinforced  fibrous  composites,  con- 
sidertible  success  has  been  achieved  in  the  development  of  continuum  models 
based  upott  asimptotic  expansion  techniques  in  which  the  ratio  of  the 
ch.ir.ic  t er  i st  i c lengths  of  the  structurin';  to  the  wave  I I'lij'.t  Its  is  .issumed 
much  sm.iller  th.in  unity.  One  appro. ich,  uf  ilizi'iv;  direct  .isymptotic  expan- 
sions. has  been  proposed  b\'  Ren-Aiiioz  lltJ.2,HJ.X|  and  is  appropriate  for 
problems  of  the  w;ive-v;uide  type.  Attother  technique,  utili^inj;  sp;iti;il 
and  asNtiqu ot  i c exptins  i ons , h.is  been  proposed  by  llei;emier  and  developed  in 
a series  of  papers  b\  llev;emier  and  Niayfeh  [irj.  lle.yemicr,  llurtman  and 
Ntiyfeh  |ll''^|,  llejtemier  and  Hache  |ll;),  HI  .111  l)|  , and  Oiirtmiin,  etc.  lOl.X)  . I'he 
latter  tipplied  to  problems  of  both  the  wave->;uide  and  wave-reflect  tyites. 
rite  technique  developed  b>  lleijemier,  et  al.,  models  a het  eroj;eneous 
comiHisite  as  a continuum  with  iiiicri'st  ruclure.  In  this  theorv,  the  p.overn 
iiii;  equ.it  ions  .ire  completeli  determined  from  a knowledy.e  ot  the  j'.eomet  ry 
and  constitutive  rel.it  ions  of  the  composite  microcomponents.  In  addition, 
this  theory  provivles  i r.form.i  t i on  on  stress  ;ind  d i s]il  acement  fields  within 
the  microcomponents  ot'  the  composite.  A ti'iiical  example  is  the  elastic 
waves  in  lamin.ited  coiiqMv.ites.  We  will  briefly  |>resent  the  case  of  wave 


pro|iay,at  i on  norm.i  1 to  the  laminate  in  the  followini;. 
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I'or  an  elastic  bi  laminate,  let  and  be  the  half  thickness 

of  the  layers  and  y be  the  distance  perpendicular  to  the  layerinj;. 

Ilejjemier,  et  al.  [ll?)  started  from  the  equations  of  motion  and  constitu- 
tive relations  for  the  individual  layers  and  expanded  the  stress  aijd  dis- 
placement into  i)ower  series  of  tlie  local  coordinates  with  origin  at  tl\e 
centroid  of  the  constituent  layer.  by  imposing  the  continuity  conditions 
of  the  stress  and  displacement  across  the  layer  interface,  tltey  obtained 
the  equations  in  a differential-difference  form.  Finally,  assuming  that 
all  difference  expressions  admit  Taylor  series  expansions  in  the  quantity 
A = b*'*^  ♦ h^“\  which  is  the  half  thickness  of  a unit  cell,  they  converted 
the  d i f ferent ia 1 -d i f ference  equations  to  partial  differential  equations. 

After  some  algebraic  and  trigonometric  manipulations,  the  two  partial 
differenti.il  equations  are  shown  to  satisfy  the  following  global  differen- 
tial equation  for  4'  and  4'  which  are  essentially  the  stress  and  the  strain: 


|coshi2't*‘ * ^ I .l^lcoshi^y'-*'^  i 3^)  ♦ sinhfJy*'^^  e 3^)sinlu2Y'^“'^  c 3^) 

- coshi2  t:  3^)}|J  \ = fl 

where 

= y/^  . r = ct/«  . E = A/ii 


i2.3) 


In  F.q.  (.2.^),  t is  a reference  length,  c^^'^  and  p^^''  are  the  wave  speed 
and  mass  density  of  the  ath  l.iyer,  .iiul  c is  the  wave  speed  of  the  com- 
posite when  t « 0.  In  the  limiting  c.ise  E * 0 . I.q.  (221  reduces  to 


>-r 


c is  obtained  by  letting  the  coefficient  of  3| 


t2-n 


.iiul  hence 


to  be  unit > ; 


(25) 


l;i|iiat  ion  (22),  when  iwpamls  in  iiowers  of  e.  can  bo  written  as.  after  makin.o 
use  of  liq.  (.25), 


where 


|U  > a,c‘3^  ♦ * ■..)'^i 

- (1  ^ b,.‘'3';  . = 0 

a,  = 2V).'  . a^  = 2-Vo! 

I.,  ■ jn'"'  . ,‘-'-'1  - n'"'  - 


(2(.) 


(27) 


y 1 - (Y 


(1  )^’  ^ , (2)‘- 


2 m*'  (2)“'  ■)  (1)^  (2)^ 

^ fr,  ir  * y'  ) * ,,  y'  y' 


This  was  obtainotl  b\'  lle.i^emier.  et  al.,  Ill7|.  However,  b obtaiiun.)  liere  is 

H 

different  from  that  of  [Hill,  and  it  seems  that  li  j of  |II1I|  is  in  error. 

One  could  hate  I'l'tai ned  a genera  1 expression  for  the  coefficients  a, 

2n 

and  1>  ol  l-q.  (.d”)  it'  one  rewrites  lii|.  (22)  in  the  I'ol  lowin.v;  form: 

|o  cosh(2>tiy^)  - (()  - 1 )coslu2iS  c 3| ) - '-■osh( 2 c 3- )|  ^ ^ =0  (2S) 


where 


Y - ^ Y^-'’ 


0 = ( I ♦ ICl  2/. lie 


tS  . 


(20) 


Noticin.i;  tltat  I'.q.  (25)  can  lie  written  in  tlte  form 
0 V ■'  - 1 1'  - I ) ^>1  - I 


(50) 


)quation  (2S)  reduces  to 


l.S 


|cosh(.2  t-  3.")  - [ ‘ ■ 1 — ^ cosh(,2  ^ ^ 


.2  ^ .2 

r - 

lixpaiision  of  each  tenn  in  powoJ's  of  c leads  to 


y Lov.;ti-  \ J'\  ■ — TtV ' J ’) 

jJ(-0  (.'.'1 


Therefore , 


,J.ul 


'‘y..  = 


Jn 


1. , - — •- 

Un+21! 


r’  - «s' 


It  can  be  shown  tltat  T.q.  1^33)  reprixiiices  lyq.  (27)  for  n = 1 aiul  2. 
by  1 et  t i ng 
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{?}’{$}  litvU -Cj^tll 


l3-n 


wltere  c^  ;'nd  k are  tlie  nondimens ional  phase  velocity  aiwl  wave  number, 

lu).  (221  reduces  to  the  exact  frequency  equation  obtained  by  Rytov  [R13) . 

On  the  other  hand,  if  we  substitute  T.q.  (34)  into  l.vi,  (2o)  . then  oite  obtains 

various  approximate  frequency  eqii.at  ii>ns  deiH'iulinj;  on  how  man>'  tenns  in 

2N 

T.q.  (2u)  are  retained.  It  we  keep  all  terms  up  to  c , we  may  call  the 

ap|>rox imat  ion  Nth  order  theory,  llejtemier  llUl]  used  a different  defiiti- 

t ion  for  the  order  of  approximation.  Tor  the  Nth  v>rdcr  theory,  he  uses.! 

2N  2N»  I 

Tq.  (22)  and  kept  the  terms  r and  r in  the  power  series  expansion 

of  cosh(  ) and  sinh(  ).  With  this  definition,  the  first  order  theory 
(N  « 1)  wotild  include  not  only  the  r^  terms,  but  also  simie  (not  all)  terms 
of  i'*  and  I,*'.  Numerical  ex,im|'les  show  tliat  his  first  order  theory  yields 


! 


I 

t I 

! 


better  accuracy  than  several  existing  theories  of  the  same  order. 


IJ 
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It  should  be  mentioned  that  the  binary  mixture  theory  of  Hej;emicr, 
et  al.,  outlined  in  Section  I.S  can  be  obtained  from  a modified  first  order 
theory  of  llejicmier. 


1 . S Vjij'  i a t i ona  l_^Met  hi^s 

I'or  harmonic  waves  in  a composite  wiili  a jieriodic  structure,  a varia- 
tion approaclj  may  bo  employed.  This  can  particularly  become  a very  effective 
tool,  if  one  uses  a variation  statement  in  whiclt  not  only  the  displacement, 
but  also  the  stress  field  is  >;iven  independent  variation.  Moreover,  by 
peniiittinjj  d iscont  inu  i ty  in  tiu'  displacement  and  tlie  stress  test  functions, 
one  can  expect  a more  accurate  reproduction  of  the  local  variation  in  the 
d i-^placement  ami  stress  fields  within  and  across  the  constituent  materials, 
lixamples  of  such  calcul.it  ioi\s  can  be  found  in  a [i.iper  by  kohn,  et  al. 
where  the  theorem  of  stationary  potential  enerjiy  which  leads  to  the  Ua>lei};h 
quotient  for  tlie  ei.nen-freviuency  is  usisJ,  .ind  in  a thesis  by  Wheeler  [WIO)  , 
and  in  another  thesis  by  Wu  IW1"|.  Nemat-Nasser  iNld]  developixl  more  nv'neral 
variation  principles  in  which  the  displacement,  the  stress,  and  the  strain 
in  one  case,  anvl  the  displacement  aiul  the  stress  in  another  c.ise.  are  j^iven 
independent  variations  and  which  include  .ippro|>r i.it e general  lunindary  and 
discontinuity  conditions.  Here  we  illustrate  Nisnat -Nasser' s variation 
principles  by  using  the  one-dimensional  case  .is  follows. 

Tor  waves  propagating  in  an  elastic  medium  whose  properties  vary 
periodically  in  the  direction  of  propagation,  i.e.,  the  x-direction,  let 
a be  the  period  icity- length.  I'hen  one  has 


Pix+a!  = p(.x1 


i 


rux*a)  = nix) 


where  n stands  for  when  dilutiunal  waves  are  considort\l,  and  for 

M when  shear  waves  are  considered.  Consider  harmonic  waves  propajjat in^ 
normal  to  the  layers  of  a composite  consistinji  of  periodically  elastic- 
layers  bonded  together.  Assume  that  a ty^>ical  cell  in  t!\is  composite  con- 
sists  of  two  materials,  S=  l.-?i  where  occupies  tl>e  re>;ion 

- a/ 2 s X < - b/2  and  b/2  s x < a/2,  and  M“  occupies  the  re>;ion 

- b/2  < X c b/2.  Now  consider  the  functional 

a/ 2 

= I I T Hoo*  ♦ -V  plO*UU*  - O ♦ C.C.|  dx 

in  which  the  superscript  star  denotes  the  ctwiplex  conjugate , the  term  c.c. 
stands  for  the  complex  conj\ixate  of  quant ities  which  precede  it,  X is  the 
l.a g rang  i an  mu  1 1 i p 1 i er , and 


1 

fl  ’ 

i-’) 

0 = ao^  ♦ 

(.1-alo^ 

"o* 

= RU") 

and 

i-S?)  , 

and  a is  a 

weight 

<u>-  u^-^u^‘^ 


.1  ^ i s t hen 


6di  = 


{ * p(0^uj6u*  + c.c.  I dx 

*> 

-([o(ii)-x]6u*(§)-[o(^4)-^^'‘‘^‘]^u*  (-4) 

- ["(j)  ' (- y)  vSX*  ♦ c.c. I 

- {<o>  .VV  - <u>  .So* 


1() 

where 

u = (l-a)u^  + au^  , (39) 

( 

The  vanishing  of  6Jj  for  arbitrary  variation  of  the  indicated  quantities 
then  guarantees  the  satisfaction  of  the  field  equation,  the  quasi-periodicity 
conditions,  and  the  continuity  of  the  displacement  and  the  stress  across  the 
two  materials  within  the  cell. 

Hence,  by  choosing  the  appropriate  test  functions  u and  o,  in  the 
form  of  Fourier  series  for  this  case,  one  is  able  to  calculate  the  frequency 
and  so  the  Fourier  coefficients  from  J^,  (36),  by  putting  = 0.  For 

this  case,  numerical  results  obtained  by  Nemat-Nasser  showed  an  extremely 
rapid  convergency  to  the  exact  results. 

I . 6 Lattice-type  Models 

In  a fiber- reinforced  composite  the  fibers  act  as  wave  guides  for  a 
wave  propagating  in  the  direction  of  tiie  fibers.  For  a wave  propagating 
normal  to  the  direction,  ti\e  fibers  undergo  little  deformation  and  essentially 
act  as  obstacles  interacting  witii  eacli  other  and  with  the  surrounding  medium 
in  a manner  wliich  is  similar  to  the  behavior  of  mass  particles  in  a lattice 
system.  These  observations  have  motivated  kinematical  assumptions  regarding 
the  deformations  of  the  reinforcing  elements  and  of  the  matrix  material  which 
arc  analogous  to  those  used  by  the  phycists  in  wave  guides  and  lattice  models 
respectively.  Thus,  the  fibers  arc  considered  as  long  and  slender  structural 
elements  and  the  matrix  is  replaced  by  a system  of  springs.  A three-dimen- 
sional theory  of  this  type  was  first  worked  out  by  TUrhan  [T15] . A lattice 
model  simulating  a periodic  structure  of  laminated  plates  which  are 
formed  by  a redistribution  of  masses  and  stiffnesses  of  fibers  was  first 
presented  by  Drumhcller  and  Sutlierland  [tHO]  . Related  works  will  be  reviewed 


later  on. 


r 


\ 


1.7  Micromor^>hic  Iheory  of  Contiima 

A theory  of  mi eromorphic  cent inua  has  been  developed  by  Uringen,  et  al. 
[I’.d.Tlbl,  and  Habip  [HI]  in  a sf  ies  of  papers.  The  theory  is  intended  for 
the  prediction  of  thermodyn;unic  behavior  of  granular  solids,  anisotropic  and 
polymeric  fluids  and,  in  particular,  composite  materials.  In  this  theory  the 
mechanical  fields  are  considered  .is  distributions.  Partial  differential 
filiations  jioverninj;  the  moments  of  fields  up  to  any  order  have  been  derived. 
In  this  theory,  a smoothiii};  operation  is  also  employed  in  whicli  sums  over 
individual  constituents  are  approximated  by  intei;rat  ions  over  the  entire 
material  volume.  I'his  tlieory  is,  in  effect,  a non-classical  mixture  theory. 
As  witli  the  general  mixture  theories,  the  general  forms  for  interactions 
and  constitutive  relations  are  postulated.  However,  the  unknown  functions 
and/or  constants  involved  must  be  determined  from  experiments. 

1 . 8 The  Neighborhood  Concept 

A neighborhood  concept  related  to  the  differential-geometric  method  in 
the  continuum  theory  of  dislocation  has  been  proposed  by  Ben-Amoz  IB21).  In 
this  theory,  one  avoids  the  difficulties  associateil  with  discontinuous 
material  properties  by  utilizing  a neighborhood  averaging  technique.  Unfor- 
tunately, the  relations  between  the  displacements  of  the  constituents  and 
the  corrcsixinding  neighborhood  averages  must  be  postulated  and/or  deduced 
from  experiments.  So  far  this  technique  has  not  been  further  developed. 

1 . 9 I'heory  of  F.l ji s t um ith  Microstructurc 

In  the  continuum  mixlels  afori'inent  ioned  the  mixlia  dealt  with  all 
possess  a periixlic  microstructure  such  as  media  witli  equally  spact\l  fibers 
or  peruxlic  laminae.  The  treatment  in  such  cases  is  greatly  facilitatcxl 
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by  tiu'  oxistinj;  periodicity  which  eiiahles  tlie  derivation  of  continuum 
titeories  l>ased  on  an  analysis  of  the  m i croinot  ion  in  a unit  cell.  While 
substantial  progress  has  been  made  with  media  possessing  a periodic  micror 
structure,  little  progress  has  been  made  with  media  lacking  periodicity, 
such  as  inclusions  of  arbitrary  geonu'try  embedded  in  a matrix  material, 
for  media  lacking  periodicity,  the  theory  of  elasticity  with  microstructure 
developed  by  Mind  1 in  [Ml  | is  certainl>'  an  effective  tool  provided  the  matri.x 
is  an  isotropic  elastic  material.  In  this  theory,  a set  of  equations  for 
the  macro-motion  that  contain  in  some  measure  the  effects  of  tlie  micro- 
motion has  been  deduced.  Recent  1>’,  Ben-.\moz  [B.SOJ  has  extended  the  theory 
of  elasticity  with  microslructure  to  a lieterogeneous  medium  consisting  of 
inclusions  of  arbitrary  geometry  embedded  in  a matrix  material,  llis  treat- 
meitt  is  structured  along  the  lines  of  Mindlin's  theory  although  there  are 
important  differences.  The  crucial  difference  is  that  the  arbitrary 
material  constants  in  Mindlin's  theory  are  deduced  here  in  teniis  of  known 
constituent  properties.  Specifically,  two  pairs  of  characteristic  constants 
are  identified;  both  length  and  time  scales  associated  with  dilatational 
and  shear  waves.  In  the  following,  the  dpiamic  theory  for  composite 
materials  of  Ben-i\mo:  is  outlined  briefly,  for  wave  propagation  in  a 
heterogeneous  medium  consisting  of  inclusions  of  arbitrary  geometry  em- 
bedded in  a matrix  material,  Ren-Anioz  has  obtained  the  d isplacimient 
e<(uat  ions  of  motion  as 


(40) 


. . , 1 . ^ V')V‘U. 

lo  / I 


\ = \ / ^ \ ♦ Ji*  1 . c ^ c /(\  * ji:  1 

f'*  = (a’V  |/(\  ♦ 

O V V 

■ 1 hmv,.  4,,.u>'.,->va;j 
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In  the  t'oroj^oinj^.  \^.  .nul  roprosont  the  inclvisien  Jonsity  and 

moduli  whereas  represent  the  correspond i nv;  mati'ix  properties; 

V,.  and  \ denote  the  volume  ('iMctions  of'  inclusion  and  matrix  in  a unit 
t m 

cell,  respec  t i ve  1 >■ . The  intei;ral  1^.  taken  o\er  the  inclusion  volume  \' ' 
in  a representat  i\e  volume  V is  the  pol.ir  moment  ot'  inertia  ot'  the 
inclusions  about  the  center  of'  the  representative  volume  and  thus  the 
effect  of  inclusion  distribution  is  contained  in  this  integral  whicii  enter; 
into  the  material  constants  a.  c.  1.  I'roni  ^10'.  ben-Amo;  h.is  extracted 
the  following  two  svstems  for  rotational  aiul  ^1  i 1 at  at  i ona  1 modes; 


0 (I  - d*V-1V\o.  = ll  - ! tl  ♦ V'U'-u'. 

v 1 ' ' 'tile  1 

(1  -ii-'V-)V'e  = (1  - (le>-r*  V'|e*  j fl  ♦ V'lV  e 


I-'-' 


vshorc  uK  :iro  the  rotations  aiul  o is  the  <.1  i I vit  at  ion  aiul 
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The  systems  of'  equations  are  valiil  up  to  wavelengths  ot'  the  order  of  a unit 
cell  dimension  .ind  are  fouiul  to  reduce  under  special  assumptions  to  Mindlin's 
eijuations  in  the  lonj;  wave  approximation.  for  a harmonic  plane  vtavc  propa- 

'v*  J 

^i;atin>;  in  tlve  medium,  the  dispersion  curves  for  a cuhic  array  of  spherical 
particles  ohtaiived  by  this  theorv’  are  rather  similar  to  the  curves  sketclicd 

U. 

in  l^ll]  ftir  the  lowest  acoustic  modes. 

1 . U'  V i s c oe_l  a s t i c Ana  I oj^i^'_s 

model  for  the  predict  iiHi  I'f  the  dispersive  effects  i n laverctl  com[MS  i 1 1' 
materials  baseil  upon  a viscoelastic  .in.iUyo  lias  been  proposeii  I'v  barker 
[ B(i  I . The  model  civnsists  of  a particular  stress  relaxing  et|uation  of  state 
of  the  Maxwell  type.  I'he  parameters  involved  are  defined  in  terms  of  the 
proi'erties  of  tlve  const  itvient  materials  and  ceomet  r>  of  the  layerv'd  com- 
posite. The  technique,  which  is  sem  i - emp  i r i ca  1 , predicts  almc>st  exact  Iv' 
tlve  averavtc  stress  in  a unit  cell  of  the  laminate.  In  effect,  the  moiiel 
smootlves  out  the  detailed  behavior  arisinj;  from  reverberations  in  the 
layers  of  the  composite. 

1.11  Di.screte  Continuum  fheory 

A discrete  continuum  theory  for  periiKlical ly  layered  composite  materials 

I , 

has  been  jvroposed  by  Chao  and  Tee  IC.s).  The  treatment  is  more  or  less  aloni; 
the  line  of  Achenbach,  et  al.  |A1>1  but  without  usinv;  the  smootlvin^  process. 

In  this  theory,  the  displacement  field  for  each  layer  is  obtained  by  develop 
inn  two-term  truncated  Taylor  series.  The  i>‘’verninn  equations  which 
incorporate  interface  continuity  comlitions  are  deriveil  in  the  form  of  a 
system  of  vl  i f ferent  i a 1 -d  i f fe  rence  equations.  Ai>plication  is  maile  tc>  propa- 
y.at  uni  of  plane  harmonic  u.ivi's  in  an  unbotnuled  layered  medium.  Thicknes  . 

I 

t 


twist  vil>r;ttious  aro  stmliotl.  Numerical  results  predicted  by  this  theory 


a^ree  ijuite  closely  with  the  exact  results.  In  j’eneraJ  , agreements  are 
even  hettei-  than  effective  stiffness  theory,  as  wave  number  gradually 
i nc reuses . 


1.12  Statist  i ca  1 ApproacIt 

Methi)'  1M21,  hose  and  Mai  11141, lfl2|  liave  proposed  a statistical  approacli 
for  1 on,c  i t ud  I na  1 waves  of  botli  compress  i ona  1 aiul  shear  types  in  a fiber 
reinforced  comi'osite  wliere  fibers  are  randomly  di  st  I’i  but  eti  I'ut  of  iilentical 
propei'ties.  In  this  theory,  tl>e  composite  is  considered  to  l)e  statistically 
uiutorm.  I'he  phase  velocit)'  .ind  damping;  of  the  ;ivera};e  waves  are  ol’tained 
by  a st.itistical  cons  i der.it  i on  as  functions  of  the  statistical  and  the 
mechanical  parameters,  (au're 1 a t i ons  in  the  positions  of  the  fibers  is  intro- 
diiceil.  The  theory  leails  to  llashin  and  Rosen's  formulas  llU’l  for  bulk 
modulus  and  shear  modulus  it'  the  correlations  are  ij’.norcd.  Ihe  correlation 
terms  have  a s i y.n  i f i cant  effect  on  the  dampinc.  property  of  the  ci'inposite, 
especially  at  hii'.h  frequencies  and  concentrations.  The  effect  is  to  increase 
the  veUicity  and  decrease  the  specific  damping;  capacity.  .".iettler's  121] 
mean  wave  techniipie  for  laminated  random  media  and  Kriimhansl's  |K‘>|  averay.e 
four  i er- f I oquet  methoil  for  disordered  composites  are  simil.ir  to  this  approach. 

1.1.1  llvil  rodynam  i c Concept 

In  a series  of  papers,  I'sou  and  t'.hou  [ria.l'lll,  I'orvik  [ 11  1 1 , tduni  and 
Waiijt  |1’<)1,  Munson  and  Schuler  iMlS]  developed  a theory  based  on  the  flow 
across  a selecteil  control  volume  of  the  medium  to  predict  the  lluy,oniot  curve 
of  a shock  moviny  in  the  f i her- re  i nforced  and  laminated  composites.  rhe>' 
derived  individual  mass,  momi'nt  urn  anil  enerjty  conservation  equations.  Ry  this 
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theory,  they  were  able  to  determine  not  only  the  average  Hugoniot,  but  also 
the  integral  of  the  interface  shear  stress  over  the  width  of  the  shock.  In 
some  cases,  they  obtained  good  comparisons  between  theoretical  and  experi- 
mental results  on  a variety  of  composite  geometries.  This  hydrodynamic 
approach  seems  to  be  a useful  tool  for  determining  the  range  of  response 
that  may  be  expected  under  compressive  shock  loading.  Related  works  will 
be  discussed  later. 

We  have  thus  summarized  most  of  the  published  theories  in  this  chapter. 
We  will  review  published  work  which  used  these  theories  to  study  harmonic  and 
transient  waves  in  composites  in  the  following  two  chapters. 


A 
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11.  HARMONIC.  WAVCS 

TIk'  wurk  oil  h;i  imon  i c or  s i mi.'^oi  cl;i  1 wavo  propajjat  i on  in  composites  con- 
sists of  wavi'Kiiiilc  analyses,  in  which  tlie  >;eomet  r i ca  1 cross  section  of  tlie 
composite  does  not  vary  in  the  propa>;ation  direction,  and  wave-reflect 
analyses,  in  whicli  tlie  material  jiroiierties  vary  periodically  in  the  direction 
of  iiropak^at  i on . In  the  w;i ve,i;ii i ile  c.ise,  for  wave  propa>;atinR  in  tlie  z 
direct  iiMi,  .1  t>'pic.il  response  function  f is  exiiressed  as 

fl.\,y,z,tl  = I'U.yl  o\p  i(kz-a't)  pM ) 

and  in  the  wave-reflect  case  as 

f(.\,y.:.t)  = l-l.\,y,z')  e.xp  ilK.z-u't)  (-15) 

where  I'  is  the  m wle  shape,  .\,  ,v  and  z are  the  coordinates,  k is  the 
w.ive  mimher,  ui  the  freipiency,  and  t is  time.  In  case  of  wave-reflect 
type  propav;at  i on , 1'  has  the  same  I'erioilicity  as  the  jjeometry  in  the  z 

direction.  The  dispersion  of  the  waves  is  cxpressoil  in  teniis  of  the  rela- 
tionship hetweeii  any  two  of  the  quantities  c.  wi  or  k,  where  c = ui/k  is 
the  phase  velocitv.  Tor  wave-reflect  case,  la  is  periodic  in  k with 
period  -i'/a,  where  .i  is  the  lenyth  of  the  unit  cell  in  the  iiropaj;at  i on 
ilirection.  These  two  types  of  analyses  of  harmonic  waves  in  composites, 
based  on  various  methods,  are  discussed  as  follows. 

II.  1 T..vact  Theories 

The  initial  study  of  harmonic  or  sinusoidal  waves  in  laminates  was 
made  hy  Rytov  iRl.'j.  Ry  use  of  elasticity  theory,  he  ohtaineii  the  ex.ict 
solutions  for  the  case  of  dilatational  waves  propai;;it  i nj;  normal  to  the 
laminates  (w.ive  reflect  iini  prohlem)  and  the  case  of  symmetric  waves 


pro]>ai;at  i HR  parallel  to  the  lavers  (w.ivey.uide  prohlem).  Ile  presented  the 
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pl\;iso  vohK'itN'  s|>cv.'t  nim  for  oach  case.  TIu'  :oro  t'l’c'niK'iii')-  limit  i>l'  tho 
prini.ir>'  mkhIo,  whicit  oiirrosi'oiiJs  to  tho  statio  olastio  solution,  was  also 
ol'tanu'il.  Uytov's  oxaot  solution  is  widoly  ust\i  as  a basis  I'oi’  ostimatinv', 
tl\o  aocuraoy  of  tlio  tlioi'rios  of  oont  inuum  for  compos i t os , 

Sun,  Aclionhach  anil  llorrmann  |SJ(>|  , anil  Aohonliaoli  1A11|  iliscussoil  tlio 
t imo-Marmonic  waros  in  layoroil  oomposito  matorials  propar.at  i ni'.  in  t ho 
iliroction  of  t ho  layorino,.  Tiiox-  cons  iiloroil  a moilium  of  a 1 1 o mat  u\)',  l.iior-- 
of  two  ilift'oront  homoi^onoous  matorials.  Hsiny,  tho  solutions  ol'  t ho  oipia 
tions  I'f  ol.istii'it'  roprosoiu  i tij;  pi, mo  ( e mo- h,i  rmoii  i i-  w.n  os  , fh<'y  ih'rixoil 

o. v.iot  ilisporsion  rol.it  ions  1‘or  hvitli  iM'.os  ot'  simmotrii'  iloformat  ions  .nul 
.mt  i svmmot  r i i‘  iloformat  i ons  . rho  rosults  obtainoil  show  that  for  tho  hu’.li 
valuos  ol'  tl\o  ratio  of  tho  slu'ar  moilul  i . or  t lio  layorini;  stiffnoss,  tho  itis- 
porsion  ourvos  ilopart  sharpli’  from  tho  limit  inp,  phaso  \olooitios  for  "lonp 
ivavos"  at  vor\'  small  wavo  numbors.  Thus,  thov  oono  huloil  that  tho  ap\il  ioabi  1 i t y 
ol'  tho  ol'foctivo  moilulus  thoory  for  wavo  I'ropapat  > i practical  laminates 

IS  \or>  limitoil  sinco  it  cannot  .iccount  tor  il  i si'ors  i on . 

I’uppo,  i'onp  aitil  llaonor  1 1'7 1 .tnaliroil  simisoiilal  wavo  propapal  i on 

p. ir.illol  to  fibor  iliroction  in  a uniili  rod  ional  fiber  roinforcoil  composite. 

I'hoy  moiloloil  tbo  hoxaponal  array  of  circular  fibers  in  a matrix  b\  a 
circular  fibor  with  concentric  cyliiulor  of  matrix.  They  then  solvt'il  tho 
concentric  roil  problem  b\'  the  mothoil  of  o 1 ,ist  i c i t s , It  appears  t'rom  their 
numerical  rosults  tliat  tho  calculation  of  tho  parameter  a,  which  is  tho 
coefficient  of  tho  soconil  term  in  tho  expression  ol'  tho  phase  velocity  in 
terms  ot  wavo  number  w.is  somewhat  ipiost  ionablo. 

Ilottm.m.  {'f  .il.  |ni  '|  .also  diil  the  .in.ilisis  of  the  concentric  rod, 
which  models  the  hex.ic.on.i  I .trr.iv  ol  iircul.ir  libers  in  ,i  m.ilriv  m.ileii.il, 

,ind  comp.iroil  it  with  the  results  for  the  eijuiv.ilenl  l.imin.ito  for  pl.iss  epoxi 
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const  I tuout  s . A vor\  siiiut.ir  1 1 1 .it  i vc-  initurt'  iii.iv  hr  ‘-ren  hut  thr  I'arjimrt  rr 
a for  thr  two  jjromrtries  is  uiKn)ual. 

Kriitrr  (R7|  stiuliovl  thr  il  i sprrs  uxi  of  flrxiiral  waves  in  circular 
himaterial  cylinders.  He  obtained  the  j;encral  form  for  displacements  and 
the  frequency  equation  for  the  first  mode  of  flexural  wave  projiasat ion  in 
an  infinitely  lonj;  circular  himaterial  cylinder.  The  theory  follows  the 
toclinitiue  developed  by  I’ochliammer  and  Three  for  elastic  bars  whicli  is  well 
known  in  the  theorv  of  elasticity.  Thee  presented  sev'er.il  fi  rsf -branch 
dispersion  curves  for  the  first  flexural  movie  for  various  ratios  of  tlte 
constituent  cylinder  radii.  Dispersion  characteristics  significantly 
different  from  those  predicted  hy  the  theory  for  homo.ceneous  cylinders  are 
rea 1 i :ed . 

l.ai,  Dowell  and  I'auchert  [1,11  providevl  a tliorouph  treatment  of 

propa>tation  of  harmonic  waves  in  a composite  elastic  cyliiuler.  I'hey  showed 
rivtorousl)'  lu>w  special  coinlitions  corresi>ond i n^  to  particular  combination 
of  material  properties  c.iu  be  derivevl  from  the  .general  solution.  I'he 
problem  tre.ited  in  detail  pertains  to  a composite  circular  elastic  rod  of 
infinite  lenjtth  consisting  of  two  layers.  live  central  portiotv  is  solivl 
while  the  outer  pvvrtion  is  a cylindrical  shell.  Their  numerical  results, 
in  terms  of  frequency  and  real  wave  number,  were  niven  for  a composite  ri'il 
v>f  a soft  core  with  a stiff  casinj*.  The  results  were  checked  with  the 
asymivtotic  frequency  eqiuitions  at  short  wavelenjtth.  The  explorat  i(>n  of 
the  dispersive  i>henomen.i  aiwl  the  vlevelopment  of  various  simplified  theories 
for  harmonic  waves  in  .1  composite  nnl  can  be  baseil  on  the  results  of  this 
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invest i^at ion. 
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If.uli’.  Of  .1 1 . |HJ|  ooiis  uloroil  flio  H I luli  - t \ pc  pi-op,iv;.it  i on  in  .1 

three  1.1  i mens  I viiKi  1 1 V t'i  her- re  i nt'i>reeii  eoinpositc  of  orthoj;onal  ty)'C.  Mu' 
fihers  in  the  direetion  of  propnif.it  i on  (the  :-il  i rect  i on ) were  assumed  as 
a round  fiber  surrounded  hy  a eirv'ular  slieatlf  consisting;  of  material  for 
tvhich  the  (>rt>[>ert  i es  were  i)hfained  hy  lu>mo>;eni  2 iiyi;  the  resin  m.ttri.x  and  the 
lateral  fihers  rnnninc.  in  the  twii  directions  orthoijonal  to  propac.at  i on . 

I'hei  pert'ormt'd  t))e  I'locli  aii.ilvsis  h>  I'xiiandiny  the  il  i sp  1 acement  s in  a 
thri'e  d 1 mens  I on.i  I li'iirioi'  sei'ies  ,ind  dt'ruiny,  an  infinite  order  mains  for 
li’ie  I o<' f f I V' I en  t s of  the  louncr  st'ries.  Die  detormin.inl  >>f  this  m.itriv 
pro\  ides  the  ilispersion  rel.it  ion  for  the  sinusoid.il  w.ivcs.  I'he  infinite 

m. 'itrix  w.is  solved  in  two  w.iys:  (1)  a perturbation  techniipie  and  (J)  a 

triinc.it  ion  technii|ue  in  which  on  1 \ .1  finite  number  of  terms  in  e.ich  of  the 

sei-ics  was  iiscil.  It  was  found  th.it  the  l.iyeriiii;  in  the  ilirection  of  )'rop.ij;.it  i on 

n. is  little  effect  on  t lie  dispersion.  fa  1 cu  1 .1 1 i ons  usiiii’  onM  three  terms 

in  the  I'oiirier  series  in  the  lateral  directions  .uul  lust  the  :.eroth 
(averai;  i ny ) term  in  the  propa.yat  ion  direction  were  carried  out  for  both 
istitropic  aiul  orthotropic  fiber-bundle  pniperties.  The  isotropic  results 
a.yree  well  with  the  known  axisymmetric  waveguide  solution  while  the 
ortliotropic  properties  i;ive  much  better  agreement  with  experiment. 

Sve  |S4J)  carried  out  an  exact  analvsis  of  time-harmonic  waves 
t ravel  in>:  obliipiely  in  a periodical  Iv  laminatcil  medium.  I'he  analvsis  was 
based  on  two-dimensional  eipiat  ions  ol'  elasticity  and  Bloch  theory.  I'is- 
persion  relation  was  obtained  for  harmonic  waves  |'roi'a>;at  i nr.  in  an  arbitr.in 
direction.  l.imitiiu’,  pliase  velocities  were  presented  for  infinite  w.i  ve  1 enc.t  h 
for  any  .iiiyle  of  pro\'aj',at  ion  in  the  I'orm  of  a fourth  oiwler  determinant.  In 
case  ot'  propay.it  ion  .iloiir.  or  across  the  l.oers,  thi'-  determin.int  reduces  to 


twi'  d«'t  ermi  n.int '.  ol  '>i\  v'lul  I'lvler  th.it  \ield  the  limit  iny,  ph.iso  vi'lovitu's 


ill  root  ly.  Ills  luimorioal  results  iiulicate  clearly  tlio  ilei'eiuienco  of  liis- 

i 

pets  ion  upon  the  anjtle  of  propagation. 


Sve  [S4.^|  also  investigated  t hermoe  last  i c waves  in  a periodical  !>■ 
laminated  medium.  By  an  exact  analysis,  he  studied  the  effect  of  thermo- 
elastic coupling  on  sinusoidal  waves  propagating  in  the  directions  parallel, 
perpendicular  and  ohliipie  to  ttie  planes  of  laminates.  The  effect  of 
thermoelast ici t V is  to  cause  complex  jn'oivigat ion  constants  to  occur  except 
in  the  case  of  shear  wave  perpendicular  to  laminates  in  which  the  response 
IS  un.iffected.  His  numerical  results  indicate  that  thermoelastic  attenuation 
IS  confined  primarilv  to  the  ipias i - 1 ong i tudi na I modes.  Ihe  phase  velocities 
.ind  mode  slia\H's  are  also  influenced  hy  the  attenuation  parameters,  especially 
fi>r  large  freipienc  ies  . 

t’hristensen  |t'.Hl  jn-esented  an  analytical  formulation  of  the  effective 
attenuation  of  harmonic  waves,  of  low  frequency,  through  layered  elastic 
medium.  I'he  effective  attenuation  is  defined  as  the  difference  in  trans- 
mitteil  energy  to  initial  total  energy.  I'his  energy  is  accounted  for  through 
secondary  wave  scattering  effects  resulting  in  pulse  dispersion.  When  the 
layers  have  equal  impedance  aiul  equal  stiffness,  inilse  attenuation  vanishes. 

As  a result  of  the  analvsis,  an  explicit  exjM'ession  for  attenuation  was 
derived  hy  a perturhation  technique.  The  theory  I'resents  a useful  analytical 
result  in  wave  propagation  in  laminated  composites. 

Tee  and  Yang  1 1 -I  | , and  Tee  |T.t|  analyzed  the  harmonic  waves  in  composite 
materials  with  periodic  structure  of  elastic  constants  and  density  variation. 
They  employed  Bloch  or  T'loquet  theory  and  treated  the  propagation  in  terms 
of  T'loquet  waves.  The  theory  was  presented  for  a laminated  compv>site  material 
and  pro|>agat  ion  normal  to  the  lamination.  They  found  that  tive  frequency 
spei'trum  has  ,i  haiuleil  structure,  comprising  pass  or  prvipagatiug  hands  and 


I'.iiuU.  It  v»,is  slu'Mt  tit.it  t t»o  t foquoiti  I at  t ho  hoimvla  i' i o>  ot  t ho 

haiiJs  oorrospi'iiil  ft'  wavo  |'ia>t'ilos  whioli  aro  normal  modes  of  vibration  vit  tlto 

i itil i V i diia  1 rolls  with  fixed  or  free  surfaces,  both  typos  occur  at  each 

limitinj;  fronuencx  . I'hoy  also  interi'reted  properties  of  1 loipiet  waves  in 

terms  of  normal  mv'de  thoor>'  and  interi'reted  the  high  froiiuoncy  limit  t'or 

I loquot  waves  in  terms  of  geometrical  optics  t>i'e  analysis. 

.‘^chooitborg  [SI]  consivlerod  plane  sinusoidal  waves  propagating  througli 

a mOkl  i um  tiuiilo  up  I'f  pi. me  l.i'ors  »'!  .inisotropic  homogeneous  liitearlv  el.istic 

m.iteri.il.  Usiitg  ,i  m.itrix  formu  1 .it  i on , he  found  the  stresses  .uul  ilisplace 

ment  s in  terms  ot  the  bound. irv  coiulitions  on  I'lie  boundarv  , ; ’0.  lie 

calculated  the  goner. iliied  transfer  function  R explicit  1>'.  Ihe  transter 

funct  ion  U can  be  thv'ught  of  .is  .i  tr.insfer  function  between  the  solutions 

.It  i'  .ind  .:  . where  n is  the  la\er  number.  R is  a function  I'f 

n 

t'requency  lO  .uul  the  material  parameters  in  each  of  the  n lavers.  In 
the  special  case  of  normal  incidence,  he  found  that  the  eigenvalues  are 
solutions  to  a bicubic  equation  instead  of  a sextic  equation. 

.Sutherland  |S.'!'l  .ina  1 \ t i ca  1 1 \ predicted  phase  velocit'  and  .it  t enu.it  ion 
tor  two  specific  composites  using  the  t i me- 1 emper.it  u re  superposition  principle 

V 

to  \ ield  the  v i scoe  I .ist  i c portion  of  the  total  dispersion  spectrum.  The 
two  composites,  lie  consivlered,  .ire  qu.irtr,  cloth  embedde(.l  in  .i  phenolic  matrix 
and  stainless  steel  reinforced  epoxx  . lie  specul.itevi  th.it  the  difference 
between  the  tot.il  spectrum  and  tlie  viscoelastic  portion  is  then  .iscrjbed  to 
the  effects  of  internal  geoiiietrv.  lie  determiiuul  the  total  spectrum  experi- 
mentally but  not  an.i  1 \ t i ca  I 1 v . I'he  results  are  ai'plicable  to  harmonic  w.ives 
tr.iveling  perpeiul i cu  1 .ir  to  the  fiber  vlirection  once  the  corresponding  tot.il 


pecti'uiii  li.i-.  been  deteniiined. 
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Kaul  and  Herrmann  [Kl]  considered  free  vibration  of  an  elastic  cylinder 
with  laminated  |)eriodic  structure.  I'he  circular  cylinder  considered  has  a 
j>eriodic  variation  of  elastic  constants  and  density  normal  to  the  ;»xis  of 
the  CN'linder.  They  then  developed  the  theory  of  torsional  vibrations  of  ' 

such  a cylinder  in  terms  of  Tloquet  or  Bloch  waves  which  are  nuasi -per iodic 
waves  and  whose  amplitude  profile  has  the  same  periodicity  as  that  of  the 
material  and  repeats  with  the  periodicity  of  the  cell.  Usinj;  l-UHiuet's 
theory,  the>  obtained  the  dispersion  spectrum  for  time-harmonic  waves  jiropa- 
catini;  in  sucti  a periodically  laminated  cylinder.  It  was  shown  titat  the 
dispersion  spectrum  iias  a band  sti’ucture,  consisting  of  pcissitii;  bands  ,md 
stoppini;  bands.  Motion  in  the  case  of  viracinj’  incidence,  and  motion  at  the 
end  of  the  coites  were  discussed.  It  was  also  shown  tliat  as  the  radius  of  the 
cylinder  teiivls  to  infinity,  tlie  torsional  waves  in  a circular  c>linder 
de.ceiterate  to  SH-waves  ii\  laminated  plates. 

1 I . J t'ar  iyU  icutal  Approaclies 

Kohn,  Knimhansl  aitd  l.ee  [KB]  employed  variational  methods  based  on  the 
Tloquet  or  Bloch  theory  to  stud>’  the  iiropagation  of  harmoitic  elastic  waves 
throujth  composite  materials  of  periodic  structure.  In  these  methods,  varia- 
tional principles  were  developed  in  the  form  of  intejtrals  over  a sitt_cle  cell 
of  the  composite.  The  variational  principles  pi'ovideil  a means  of  determinini; 
phase  velocities  and  stress  distributions  in  Tloipiet  waves  travel  inj;  through 
the  composite  unclianjted  in  form  from  cell  to  cell.  Strain  energ\  principle 
was  used  and  possible  jump  conditions  were  discussed.  Ihe  Kay  lei  sh-K  i t 
proceilure  was  applied  to  the  solution  of  the  variational  eipiation  to  calculate 
dispersion  relations,  the  phase  velocities  and  stress  profiles.  The  one- 


dimensionally  periodic  system  was  evaluateil  and  compared  with  the  exact 


i 

j 


.so 


■solution.  Mh'  rt'suits  slu'w  that  tliis  approach,  h\’  iisinj;  .'iinootli  0 i si'l  acomciit 
test  functions,  proviiles  a sat  i sfactor>'  determination  of  natural  frei|uoncies 
and  pliase  velocities,  Init  is  inadecpiate  for  stress  profile.  In  ]irinci]ile, 
tills  method  can  he  used  to  study  time-harmonic  waves  propajtat  i njt  in  composites 
withrather  i^enera!  structuring;  as  lony  as  the  composites  are  jicriodic  from 
ce 1 1 to  cell. 


Bev  i I acipia , Krumhansl  and  Lee  IB.SS]  made  a generalization  of  the 
previous  work  Based  on  strain  energy  consideration.  The  Bay  le  i gh-ll  i t z 
procedure  was  adojited  in  which  the  displacement  fields  wore  expressed  Ip' 
comiilex  Touricr  series  in  filament  and  by  simple  I'ouricr  series  in  matrix. 
Consequently,  they  obtained  stress  profiles  more  accurately  than  tliose 
when  the  wliole  displacement  field  was  expressed  by  simple  l-ourier  series 
in  the  previous  work. 

Heiil.icqua  .ind  lee  IB.'o]  utilized  different  iiolvnomial  repre- 
sentations in  filament  and  matrix  for  a simple  model  composite  problem 
comprising  a slab  of  each  material  with  fixed  boundar>-  conditions.  This 
simple  model  provided  the  same  difficult!'  of  a d i scont  i nui  t\'  in  strain  at 
interface,  but  provided  a simpler  case  for  assessing  various  appro.iches . It 
was  pointed  out  that  I'ourier  series  representation  is  less  convenient  in 
two-dimensional  and  three-dimensional  cases  when  the  polynomial  form  may 
jirove  to  be  superior. 

I'olwin  (Tlbl  also  made  an  .analysis  of  propagation  of  elastic  waves  in 
composite  materials  b)  variatioiwil  methods.  He  carried  out  a minimum 
strain  energy  calculation  using  a I'ourier  series  expression  for  the  displace- 
ment lield.  In  fite  I'lii"  d i mi'iis  j i^iu  I case,  this  met  ol  t'V'a  I uat  i(>n  is  unsafis- 

l.iclof!  lor  sfio-.s  profiles  .is  shown  in  |k;.|.  It  seems  ih.it  ^.uccessful  two- 
dimension.il  .ind  three  d i iiiens  i on.i  I c.i  I cu  I .it  i ons  b\  this  method  is  t'e.isible 


It'  vl  I scent  i nu  1 1 1 es  ,ire  permitted  in  strain  energ'  c.ilculat  ions  . 


1 


Bevilacqua  aiul  Leo  (Ba7|  presentoil  atiotlier  vai'i  at  i oiia  1 statement  l>aseil 
on  complementary  energ.v  coas  idei'at  i on  so  tliat  the  variational  intejiral  is 
expressed  in  terms  of'  stress.  Starting;  with  a continuous  displacement 
variation  in  I'ourier  series  form,  they  first  determined  the  corresponding 
stress  variation  hy  intei^rat  iny  the  equation  of  motion.  Then  inserting 
tliis  result  into  the  variation  inteyral,  tliey  obtained  a matrix  ei>;envalue 
problem  to  determine  tlie  values  of  fundamental  frequency.  The  correspondirip, 
stress  profile  olitained  exltii'its  tlie  correct  form  of'  st  ress  - j;rad  i ent 
discontinuity  at  tite  interfaces,  and  provides  an  accurate  approximation. 
However,  tlie  frequency  is  less  accuratel>’  ]iredicted  than  b\’  tlieir  improved 
strain  enerjjy  ai'proach  lUeBI. 

i.iriational  methods  were  also  explored  b>'  Wheeler  and  Mura  |W10|,  and 
Iv'u  [Wl?].  Instead  of  usini;  the  jirocedure  of  Ray  le  i i;h-Ki  t : , Wu  employed, 
however,  the  lialerkin  procedure. 

In  a series  of  pa['ers  |NlJ-l(>),  Nemat-Nasser  applied  more  iteneral 
variational  principles  to  the  Tloquot  wave  problem  than  those  considered  by 
Kohn,  Lee,  et  al.  Based  on  He  1 1 inger-Re i ssner  variational  method,  he 
developed  jteneral  variational  I'rinciples  in  which  tlie  displacements,  the 
stresses  and  the  strains  in  one  case,  and  the  displacements  and  the  stresses 
in  another  case,  are  j;iven  liuieiiendent  variations,  and  which  include 
appropriate  yteneral  boundary  and  discontinuity  conditions,  l-rom  the  jteneral 
variation  principle,  he  then  derivcil  a new  quotient,  in  contrast  to  Ra\'leii;h 
quotient,  to  iletermine  the  freipiencies  of  harmonic  waves  in  composite.  IV.ives 
propaKatinj;  normal  to  the  layerini;  in  a laminated  meilium  were  analizeil  in 
detail,  and  dispersion  curves  were  presented.  Numerical  results  show  an 
extremely  rapid  conver.i;ency  to  the  exact  results  and  are  more  accurate  than 
those  obtained  i>y  Kohn,  et  al.  .Mso,  based  on  the  new  quotient,  he  developeii 
quite  accurate  lower  and  upper  bouiuls  for  frequencies  of  the  harmonic  wave  . 


i 


J 


\i‘ni.it  111  ,uul  M i |N1^1  lui'tluM'  i lu  t i i;;i  t t\t  luirmonw 

w.ivos  in  oiu'- , two  aiul  1 1\  I'oo-il  I inons  i ona  1 ci'iniios  i t es  usiiu;  a now  quotiont 
liorival'lo  t'l'oiH  a qonoial  vai-iational  iiiotiuwl  in  wliioli  I'ot  li  st  rossos  aiui 
i s|' 1 aooinont  s aro  \ai'ioil  i lulopoiulent  I y . The  yonoral  onofo,\'  motlioil  was 
tonmilatovl  In  assuming  porioilio  struoturo  with  oontinuous  aiul  oontinuouslN 
li  1 t't'oront  i ah  I o niass  ^.lonsif^  aiui  olastio  oonstants  witliin  tho  unit  ooll.  liio 
final  fosult  is  in  a simiilo  fooiii  aiul  roprosont  s oithor  uppor  or  lower  houiuls 
on  tlu'  f roqiiono  i os . It  was  shown  that  tho  now  ipiotiont  \ioKls  upj'or  houiuls 
It'  t !io  stress  \ariahlos  are  t'roo , aiul  lower  houiuls  if  tho  vl  i sp  1 aoiiu'ii  t laria- 
hlos  are  tree.  rrv'OOi.luros  on  oriu'r  ostimatiiu'.  wore  also  pixoii.  \umorioal 
results  slu'w  tho  suporiorit'  oi'  the  now  qiu'tiont  vomparoJ  to  t lie  Kailoi.yh 
vjiiot  lent  usovl  In  Kohn,  ot  al. 

Noma t -Nasser  aiu!  Minajpiwa  |MS)  ayaiii  ooniparoil  Railoio.h  quotiont  with 
t lie  proposoil  now  quotiont  to  oxamino  reasons  for  tho  astonishing  aoouraoi 
ot'  tho  latter.  iN'iiipar  i son  loil  to  ,i  soliomo  for  ol'tainiipy  i iiiprovoii  tost  I'uiio 
t ions  wiiioli  y,ivo  \or\  .loourato  houiuls  I'or  froquonoios.  I'ho  soliome  sooiiis  to 
ho  o\troiiiol\  oft'ootivo  in  all  oases  of  harnionio  waves  in  hnorod  oompos  i t cs . 

11..''  l.ffootivo  Moihilus  I'hoorios 

Hohrons  |H1S|  troatovl  tho  propap..it  ion  ot'  olastio  waves  of  Hlooh  form 
in  .1  l.imollar  porioilio  ootiiposifo  I'or  lone  wavolon.oth  oases.  The  analvsis 
yioUls  ovpros-ions  for  tho  five  i luU'i'Oiulont  offootivo  olastio  oonstants, 
whioh  ho  oalloJ  axora.yovl  olastio  oonstants.  of  tho  material.  I'he  theor>’  is 
\alul  when  the  w.uolon^th  is  Ion;'.  ooiiipareJ  with  the  i nt  orooiiiponent  spaoinps 
ot  the  I'omposito  aiul  when  t lie  ooinposile  is  .i  porioilio  striiotiire. 

hohrens  |ltl!>|  also  anal'.'Oil  tho  propaontion  of  el.istio  waves  of  Hlooh 
t'orm  in  a filamentary  oomposito  for  lonp  wave  1 on.i;t  h.  oases.  The  I'haso 
\elooitv  .tiul  the  nine  "honuyi'.on  i rovl"  olastio  oonstant  ■ were  v>ht  .i  i noil . \);ain. 


thf  thooi-y  IS  valul  wlioii  the  \»ave  1 on>;th  is  lonn  compared  uith  the  intercom 


ponent  spacin>;s  ot  the  composite. 

Behrens  [BJO]  .i.eain  discussed  lon.c  waves  of  Bloch  form  traveling  in 
filamentarv  composites.  He  sluiwed  t liat  such  stable  lloquet  waves  or  Blocli 
waves  retain  their  form  relative  to  the  periodic  structure  of  the  composite 
as  they  travel.  He  tre.ited  the  problem  using  tlie  "metho«.l  of  long  waies" 
wliich  he  proposed  i^rev  lous  1 v . Iv'hat  he  actually  treated  is  a c\'l  indrical 
inclusion  in  a c\  1 1 lulr  i c.i  1 medium.  .i  very  idealized  condition  for  t>pical 
fiber  comi>osites. 

Using  .in  .ippro.ich  siiii’l.ir  to  that  used  in  studying  the  propagation  of 
electromagnetic  w.nes  through  random  media.  Mok  [M5i  and  Osteii  (OJj  dis- 
cussed the  sinusoid.il  dispersion  in  .i  fiber- reinforced  m.iterial.  They 
considered  isolatoil.  transverse,  cylindrical  and/or  spherical  inclusions  and 
then  combined  tliem  b\-  random  .iveraging  to  obtain  sinusoidal  dispersion 
solution  for  composites.  I'he  effective  modulus  and  densit'  were  then  obtained 
from  tlie  dispersion  solutions.  Both  vif  them  are  complex  numbers  and  ilepend 
on  wave  frequency.  I’his  fact  .ipj'arent  1>  indicates  the  existence  of  dissipa- 
tion and  dispersion  in  the  comi'osite  under  dynamic  loadings. 

Weitsman  [IV"!  presented  an  an.ilytic  treatment  of  wave  propagation  in 
unidirectional  continuous  filament  comiiosite  m.iterial.  I'he  theory  is  based 
on  i next ens ib le  fiber  reinforcement  .ind  therefore  sliould  be  applicable  to 
unidirectional  comiiosites  reinforced  with  closely  spaced  fibers  which  are 
much  stiffer  than  the  matrix  material.  The  tlieory  is  also  restricted  to 
waves  of  sufficient  lengtli  that  geometric  dispersion  is  not  signi  fic.uit  . He 
showed  the  existence  of  three  types  of  harmonic  waves  and  discussed  their 
associated  characteristics  such  as  slv'wness  surfaces,  velocities,  ['ropagation 
directions,  ili  splace'ent  s and  energy  flux.  He  also  considerevl  wave  sc.ittering 


:vi 

.it  .1  i'l.inc  i<t  tiht'r  m i s;t  1 i _i;iiiiK'nt  . (ti  ,m>neni!,  when  ;in>’  at  the  three  waee 
t>iH's  impinjteon  such  a plane,  refleeteJ  atnl  t ran.sni  i t ted  waves  of  all  Throe 
t'pes  are  itenerated.  I'ondition  under  whieli  Kayloijth  surt'aee  wave  may  exist 
was  also  diseussed. 

Ko  i t small  IKS]  t'urtlier  i nvest  i jtat  ed  ttie  reflection  of  liarmonic  waves  in 
I i her- re  i n forced  matoi-ials.  Two  cases  were  considered:  (I)  fibers  are 
inextensible,  IJI  fibers  are  almost  inextonsible.  lie  showed  that  throe 
tipes  ot  hariiRtiiic  waves  exist  for  each  case.  The  reflections  of  each  ty)ie 
ot’  w.ne  I'rom  a pi, me  free  sur  t'ace  - bound  i ny  a iialf-space  were  i lu  est  i yal  Oil , 
and  numerical  result''  were  presenteil.  I'ifferences  aiul  s i m i 1 a r 1 1 i O''  between 
the  nature  of  wave  I'ropagation  of  the  two  cases  were  discussed.  file  results 
might  be  valuable  for  composite  material  that  are  sparsel>-  reinforced  b\’ 
highly  rigid  fibers  although  the  iiractical  significance  of  the  theory  is 
limitei.1  In  the  simplifi  ing  assumi'tions  incorporated  into  the  mat liemat  i ca  1 
iikhIo  1 . 

Datta  lIH)  considereil  a unidirectional  f i her- re  i nforced  composite 
through  which  clastic  shear  waves  propagate  in  tim  direction  perjiendicular 
and  polarization  parallel  to  fiber  direction.  I'ibers  were  assumed  to  be 
randomly  distributed  and  of  elliptic  cross  section  witli  transverse  axes 
ei ther a 1 i gned  or  randomly  oriented.  Wave  equation  was  solved  in  the  limit 
of  long  wavelengths  making  use  of  results  for  single  fiber  in  intinite 
medium.  Since  multiple  scattering  was  neglected,  solution  is  limited  to 
small  fiber-volume  ratios,  for  perfect  orientation  of  fiber  axes,  he 


derived  two  effective  transverse  shear  moduli.  A limit  of  narrow  cavities 
with  van  i sh  i ng  lolume  {I'lritfilh  cr.icks)  was  discussed.  It  is  usetul  in 
mech.inics  ol'  composites. 


Woitsniaii  ai\vl  Ucnvonisto  IK'.)]  coiisiilor  tl\p  pi'opaKat  i on  ol  liannonu- 
w.ivos  in  Ctniipos  1 1 0 material  ro  i n t'orooil  by  fibers  whose  directions  vary 
continuous  in  space.  'I'he>  I'ociised  speci.il  attention  on  the  case  of  sliy.bt 
v.iri.itions  in  1‘iber  directions  aiul  on  plane  w.ives  that  propagate  j'arallel 
to  t!ie  impertiirl'ed  ilirection;  a case  in  which  t lie  solution  can  l>e  e.xpresseil 
completely  in  analytical  form.  they  concluded  with  a solution  to  a specific 
ex.iiiiple  in  which  they  assumed  that  the  f iber  d i rcct  i ons  \ .iry  linearly  with 
the  any.le  I'rom  an  initial  to  a final  position.  In  a discussion  of  the 
specific  solution,  they  pointeil  out  that  this  linear  variation  in  t'iber 
directions  y’ives  rise  to  secondari  w.ives  throiiy.hout  the  composite  body;  aiul 

these  results  t<‘iid  to  iiulic.ite  that  it  m.iv  he  possible  to  detect  infernal 
devi.itions  in  the  fiber  vlirections  by  means  of  dinaiiiic  tests.  Ihe  b.isic 
ide.is  presented  in  their  work  follow  that  classical  formul.ilion  ol  w.ive 
propay.ation  in  anisotropic  materials.  The  f i be  r- re  i nforceil  composite  w.r. 
ri'presented  by  means  of  an  eipiivaleiit  transversely  isotrojiic  medium.  I'he 
solution  is  valid  for  lonp,  wa  ve  1 ei4;t  hs  . 

.Shoenbery.  and  IVeitsm.in  |S11|  further  sluilietl  plane  hariiioiiic  waves 
propayat  i lie,  in  eoiii|iosile  reinforced  by  I'ibers  with  |'er  i od  i ca  I I y variiny 
vlirections.  The  fibers  were  .issumed  to  wobble  per  i ovl  i ca  1 1 >■  about  a 
dominant  ilirecliini  aiul  all  fibers  were  assumeil  to  be  parallel  to  each 
other.  The  t'iber  reinforced  composite  was  represent  eil  by  an  eipiivalent, 
transversely  isotropic  medium  whose  preferreil  direction  coincides  with  the 
direction  of  the  fibers.  The  wobbliness  endows  the  m.iterial  with  a struc- 
tural |>eriodicit>  which  y.enerates  dis|H'rsion  at  all  frequencies  ami 
i IIS  t .’ll)  i I i f y ('or  various  fret|iieiicy  bands.  The  zones  of  i nsl  abi  1 i i >■  were 
analyzed  by  the  perturbation  method.  Uouiuled  solutions  and  unboinulevl 
solutions  were  il  i sciisseil . I'he  existence  of  iinboiiiuled  solution  is  eipiivalent 


tv'  tin'  "pi '■  I ti"vi>ii.mv'('"  oiivv'imt  (' I'OvI  in  v'lht'i'  vl\ii.»miv’  pfv'I'lonis  . Hit' 

.111,1 1 \ s t IS  vtilivi  i>n  1 N'  I'v'f  iv.ivt' U'liv'.t  hs  w1mv.'Ii  .iro  Itii'nt'i'  than  h\  .in  v't'vinr  ni' 
iii.u'.n  i t tivh'  tin'  v.'cll  vl  i iiions  i v'lt  v't'  t ho  ov'iiipt's  i t e' . 

I I.  I I rt'ovt  i VI'  St  jlt'iii'S';  riu'i'i'ins 

As  lilt'll  t I v'lu'vi  I'.irlit'i'.  t lu'  t'l'l'i'vt  i VO  St  I I'l'noss  t’.u'v'vv  ivas  t Ito  first 
vi'iit  1 niniiii  iiK'ilol  t'or  l.imin.itovl  iiiovlia  .iiul  f i her- ro  i n fi'rvovl  Vv'iiipv''.  i t os  tv' 

.ivvv'iint  fv'i  p.ov'iiiot  r i V vl  i spt'i''i  i v'li , I'l'r  .i  l’ i I'or- t o i n 1\' rvovl  Vv'iiipv''- i t o a 

'.iiiiplt'  fv'i'iii  V't'  t ho  thov'i'N  iv.is  viovolv'povi  hy  Avhonl'.ivh  .iiul  llot'i'iiianii 

|ll!l.\ll|  Thv'v  v'liip  I v'\  vvl  Ihv'  thv'v'i'v  tv'  ‘.tiiviv  t In'  pi'v'p.ii',.il  ii'ii  v'f 

pl.iiio  h.i  niiv'ii  I V'  iva\t".  Ill  tin'  vl  i rov' t i v'li-.  v'l  t ho  I'iht'i"'  .iiul  lU'i'iiial  tv'  t lit'  lihor 
vl  I I't'v  I I v'li'. . rho'  I'v'iiiivl  th.il  plant'  t |•,lnsvol■st'  u.itt's  prv'pa.i',,i  t i ii;;  in  t ho 
vli  I't'v't  It'll  v'f  flit'  I'lht'i's  ,ii'('  t!  I 'ipi'j's  n . rht'  vl  i spi'i".  i I'll  viini".  Iv'r  ht'iv'ii 
t'pv'W  .iiiil  I'.lass  t'pv'w  viop.ii't  at  siii.ill  u.ivo  niiiiil'oi's  frv'iii  t ho  vv'ii-.  t .int  ph.iso 
volv'vit'  pi'OvI  ivt  t'vl  l'\  t't't'ovt  i \ o iiiv'vlnlus  thov'i'N.  It  iv.is  .ilsv'  tv'uiul  that  t ho 
vl  i spt'i's  1 v'li  v'l'  1 v'lir,  i t uvl  1 n.i  I iv.ut".  I'v'r  small  wavo  niiiiil'or  is  noy,  I i r,  i h lo  as 

Vv'inp.n'Ovl  tv'  vl  i spt'i's  i v'li  v'l'  t vans  vt' I'st'  iv.uos.  It  .ipi'o.irs  th.it  in  v'l'vlor  tv' 

• iv'v'v'iint  t'v'i'  t ho  vl  i spots  i v'li  v'l'  1 v'liv;  i t tivl  i na  1 iv.itt's  pt'i'p,i>',.i  t i ni;  in  tho  vl  i rov' t i v'li 
v't'  tho  fil't'i's,  iiiv't'O  I'v'iiipit'v  iiit'vh'l  is  lu'ovlotl. 

It'  iiiipi'v't  o t ht'  siiiiplo  iiiv'vlol  I'v'r  f il'ot'-roinfv'rvovl  v'l'iiipv's  i t os  , At'honl'avh 
.invl  Sun  |A1.S|  vt'iis  t t'liv  t OvI  a iiu'ro  .ivviir.ito  lu'iiiv'yonov'us  w'lit  i innim  thooi''  v'l' 

tho  ol'fi'vtivo  St  i I'l'noss  t n po  I'v't'  a I' i hor- t'o  i n I'v't'i'ovl  vv'iiipv's  1 1 o . IIu'n  Ov'iis  i vlt' rt'vl 

.1  I'i  hor- ro  1 n ft'rvovl  v't'iiipv'-.  i l o t'f  iin  ivl  i rovt  i "iia  1 I'ihors  in  a roi' t anyti  1 a i'  arr.n 
t'liil'v'vlvlovl  in  ,1  iii.itriv  iii.i  1 1' I't  .1 1 . rho  iiiv'tlol  iv.is  h.ist'vl  I'li  i'\p.ins  i t'lis  ui  tho 

vl  I ''P  1 .iv't'iiK'iiI  s ati'v''''.  ropri'si'iii  ,it  I \ o v'olls,  I ho  t r.iii'.  1 1 i v'li  If  v'lii  I lit'  avtii.illi 
I nhv'iiii'v'.t'nov'ii''  vt'iiipv's  1 1 o tv'  .i  hv'iiu'r.v'iU'v'iis  t v'lit  imniiii  iv.is  .uhu'tt'vl  I"  i nt  rv'vhiv  i nr, 
oi'iit  i inniiii  ItolvIs  fvir  r.rv'ss  vl  i sp  | .ivomont  .iiul  K't'.il  vio  ft' riiv.i  t i tins  . I'ht'  o I'l't't' t 1 1 o 


stil’t'noss  theory  t'lM-  t'i  i>er- re  i i>t’oreei.l  composites  was  ttien  ileriveil.  I't)!-  the 
transverse  harmonic  waves  propa>;at  i ity.  in  tl\e  direction  of  tlie  fihers,  pliase 
velocities  were  calculated  immer i ca 1 1 y . The  results  show  that  t\>r  wavelengths 
that  are  ot'  the  order  o1'  maitnitiide  of  a t'iher-diameter  or  a distance  hetwcen 
fihers,  the  phase  \elocit\’  is  iii;»rkedly  dependent  on  the  wavelength  if  the 
elastic  constants  of  the  reinforcini;  material  differ  siil>st  ant  i a 1 1 y I'rom  those 
of  t iu'  matrix  material. 

An  analojtous  effective  stiffness  theor\’  for  fiber-reinforced  composites 
with  rectanijular  ly-spaced  rectangular  fihers  w.is  worked  out  hy  ha  rt  ho  1 i>mew 
•iikI  I'oiaik  |1U'>),  and  was  ai'plied  to  discuss  the  harmonic  waves  propa>tat  i nc, 
in  such  composites.  Impressions  for  dispersion  of  loiiit  i t ml  i na  1 and  fle\ur.il 
waves  jiropajtat  i ni;  parallel  to  the  fihers  were  derived.  Trop.iiiat  ion  acia'ss 
the  fihers  was  also  hriefli’  considered.  Aiiain,  a \ery  m.irkeil  ilispersion  w.is 
observed  when  the  elastic  constants  of  the  fihers  differ  subs t .itit  i a 1 I > from 
those  of  the  matrix  material. 

Tor  a laminated  medium,  .Sun,  Achenhach  aiul  Herrmann  | SI S , S 1 7 , AS  ,111  I , 

A1J|  derived  an  effective  stiffness  theory  in  terms  of  macroscopic  displace 
ments  plus  additional  indepeiulent  variables  which  ,v;ive  microiiiechaiiical 
variation  in  the  unit  cells.  They  emploNod  the  governinj;  equations  to 
determine  curves  relatinjj  the  phase  velocity  to  the  wave  luimher  for  harnuniic 
waves  propajtating  parallel  to  and  normal  to  the  layering,  and  they  compared 
these  results  with  the  exact  solutions.  The  lowest  ant  i sx'mmet  ric  moilc  lor 
waves  propajtat  i 11);  in  the  direction  of  the  layeriii);  show  the  st  ronj;est  varia 
t ion  of  the  phase  velocity  over  a substantial  ran);e  of  wave  numbers.  The 
limitiii);  phase  velocities  at  vanishiti);  wave  numbers  a);ree  with  t lio-.e  ot  the 
effective  moiluliis  tlieory  aiul  with  the  exact  limits.  It  was  pointed  out  that 
their  theories  of  the  type  bear  close  resemblance  to  Mindlin's  theorv  ot 
line.ir  elasticity  with  microstructure. 


ni.it  ('I' 1 .1 1 I'.i  iMiiH-t  or-'  t'l'r  i\liu'h  tlif-.o  t'rt\|tK'iK  u's  .irf  '.iil'-.  t ,ini  i ;i  1 1 \ ,it  \ .i  i' i .iikc 

UCM  O I IkI  U'.lt  I'll  . 

1 ho  i;o\onuM>',  oi(ii,i  i i oio-  ot'  tho  o t' t’oo  t i vo  '.tilt'noss  thoor>  \\0t'0  oniploxoil 
In  Sun  jSJll  to  stikh  'Hiil'.ioo  w.ivos  )'roi'.ij'..it  i alono.  t Ito  froo  surt'aoo  of  a 
la>ofOi.l  ha  I f-spaoo.  .Sovor.il  o\,im|'!("-  of  t i mo- harmon  i o siirfaoo  iva\i'  motion 
Koro  cons  i i.loro.1  i no  1 lul  i no,  \\a\0‘i  in  a laminated  laxoi  over  a rip.id  half-spaoo, 
l\a>leiith  i\avos  .ind  lo\o  i%a\os  in  a lamin.itod  laser  oscr  an  olastio  half  spare. 
I'isporsion  oiirs'os  ivero  presented  .iiul  oi'mj'arison  u.is  made  to  results  I'rom  tho 
effect  iso  moilulus  thoors  .iiul  to  ox.iot  solutunis  where  these  are  asailahle. 
These  stuns  that  the  slisporsioii  sairses  .i  ri'  in  itoosl  .ij'.reemont  vs  i t h the  ox.ict 
solutions  and  that  tho  eft'octiso  stift'noss  tlumrs  I’.ises  ,i  hotter  approximation 
th.in  tho  effect  ISO  modulus  thoors  an  ads  .int  .ii'.o  ,r,.iiiu\l  .it  the  expense  s>l  ,i 
more  complicated  iiriu'odiiri'  than  th.it  of  tho  effect  ise  modulus  thos'is. 

riie  effect  ise  stillness  thoors  isas  refiiuul  1 1'  iiicliule  seci'iul  orsler 
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terms  in  the  d i sp  1 .icemen  t exp.insioiis  In  \cheiih.ich  .iiul  llerrm.uiii  | \ I s | . Ii.' 
h.irmoiiic  u.ise  prop. u;. it  ion  in  I .iiii  i ii.i  I lul  mediiiiii  ( iinl'oiinded  I u.is  then  ■ tiidied 
In  me. Ills  ol  this  rel  lin'd  Iheors.  liie  ex.nl  .iiul  .i  pp  ro  x i iii.i  I e »l  i -.pe  r-.  i I'li  (iiise-. 
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I'latc*  was  aiialyzt'il  h)  om|)  1 oy  i iij;  this  rofiitt'd  theory.  Comiiar  i son  of  the 

dispersion  curve  for  flexural  motion  with  an  exact  curve  shows  sinii lar 

i 

qualitative  hehavior.  llie  qualitative  difference  is  due  to  "m.iterial  dis- 
persion” which  is  thus  seen  to  I’e  s i >tn  i f i cant  for  smaller  wavelenj,;ths  even 
in  tlie  presence  of  "geometrical  dispersion.” 

rhe  propavtat  ion  of  t inie-harmon i c waves  traveling  in  an  arhitrari' 
direction  in  a jH'r  i od  i ca  1 1 y laminated  medium  was  examined  by  Sve  |.s.lj|.  In 
.iddition  to  an  exact  .inalysis,  he  .ipplied  t lie  effective  stit'fiiess  tlieorv  to 
the  prol'lem  and  obtained  tlie  corresiiondinn  ap|iroximate  dispersion  relation. 

Numerical  results  also  indicate  the  dependence  of  dispersion  upon  the  anttle 
of  propa.nat  ion.  A comiiarison  with  the  exact  solution  shows  that  a,creement 
is  sat  i s factor)'  for  those  aiiylos  where  the  dispersion  is  the  strongest. 

Hi jther-order  refinements  of  the  effective  stiffness  theory  were 
developed  by  Hrumheller  .ind  Bedford  |nil|  for  a laminated  medium.  I'hey 
presented  improve<.l  il  i sp  1 .acement  and  stress  interface  boutularv  conditions 
suitable  for  higher-order  (higher  th.in  first  order!  theories  so  that  the 
resulting  theori'  is  in  a form  suitable  for  the  solution  of  ihnamic  process 
incUKiing  determination  ot  stresses.  iVo  examples  of'  time-harmonic  wave 
propagation  show  that  the  inclusion  of  stress  boutularv  condition  has 
negligible  effect  on  the  dispersion  curves  while  the  mode  sh.ipes  are 
substantially  modified.  These  indicate  that  an  effective  stiffness  tlieorv 
adequate  for  the  determination  of  stresses  must  include  stress  interface 
boundary  conditiims.  fonsequent  !>■ , they  obtained  more  accurate  dispersion 
curves . 

I he  effective  stiffness  theories  as  describeil  above  are  based  on  the 
construction  of  str.iin  and  kinetic  energv  densities,  and  the  subsequent  use 
Ilf  Hamilton's  principle.  lor  constituent  materi.ils  th.it  are  not  perfect  Iv 
, el.istic  this  .ippro.ich  encounters  il  1 1 t i cu  1 1 i es . An  .ilternative  w.iv  of  deriving 
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t ho  l>;il;nii.'o  t i ons  of  linoor  iiiomontuni  aiul  moment  of  momentum  ilireotly 

from  physical  arj;iiments  l)asecl  on  momentum  considerations  was  jiresented  hy 
Achenhacli  This  alternative  method  has  the  advantajje  that  it  can  lie 

used  for  composites  in  which  the  individual  layers  are  not  elastic. 

tirot  and  Achenhach  (lIlOl  extended  the  alternative  formulation  of  effoc- 
ti\e  stil'fnoss  theor>’  to  a viscoelastic  laminated  composite  including 
temperature  effects.  They  derived  a set  of  constitutive  eipiations  for  a 
l.iminated  medium  composed  of  l.i>ers  of  two  anisotropic  t liermov  i sciie  1 .ist  i c 
solids.  The  special  c.ises  ol'  isotropic  t hermov  i scoe  la.' t i c laiers,  anisotropii' 
thermoelastic  laiers  and  isotropic  t hermoe last i c laiers  were  brieflv  dis- 
cussed. This  nonlinear  titeory  has  tlie  advantaj’e  that  it  can  he  used  for 
composites  which  are  not  elastic. 

An  exteitsion  of  the  effective  stiffness  titeory  to  incliwle  lar.t;e 
deformations  .ind  nonlinear  material  hehavior  was  also  carried  out  liy  (!rot 
and  Achenhach  |t'ill|  for  a laminated  elastic  composite.  The  resultinjt  system 
of  nonlinear  field  enuations,  consistinj;  of  balance  enuations,  constitutive 
eijuations  and  constraint  conditions,  liear  a close  resemblance  to  enuations 
defining  a nonlinear  theory  of  elasticity  with  microstructure.  five  equations 
were  employed  to  study  the  proiiasat ion  of  small  amplitude  time-harmonic  waves 
superimposed  on  a lar^^e  static  deformation.  Tlie  results  are  very  similar  to 
those  in  linear  theory  of  effective  stiffness  except  that  the  coef f i c i ent s t hat 
•ippeareil  in  the  dispersion  relation  depend  not  only  on  the  structurini;  but 
also  on  the  lar.ye  static  deformation. 

An  effective  stiffness  tlieory  for  a comiuisite  of  a cylindrical  geometry 
was  ileveliipeil  b\  t’luni  .ind  .\chenbach  lf7|.  I'hey  presented  field  equations 
;.;o\  era  1 ny  the  mochanic.il  behavior  ol'  laiered  cylinders,  tirot  |II1J1  piMposed 
.in  ol'fectne  st  i fines  . ci'iitinuum  model  for  cun  i I inear  laminateil  comtnisites 
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of  a very  >;eneral  type.  Tlie  model  was  developed  by  expaiidinji  the  variables 
ii\  a series  within  each  layer  and  averajjing  over  the  thickness.  A constitu- 
tive theory  for  nonlinear  elastic  materials  was  then  developed,  and  special 
cases  of  the  theory  (linear  cases)  for  cylindrical  and  splierical  laminae 
were  presented  and  discussed.  It  was  expected  that  if  the  ratio  of  the 
thickness  of  tlie  layerini;  to  the  radius  of  curvature  of  the  layerinj;  is  small, 
the  approximate  theor>'  with  two  term  expansion  in  descril'inj;  the  behavior  of 

the  laminated  body  should  be  as  adequate  as  the  ainiroximate  theories  of  plane 
laminates  presented  earlier. 

brumheller  and  Hedfoi-d  (I11J|  useil  .1  iiioili  ficat ion  of  the  simpler  micro- 
structure theory  developed  earlier  for  elastic  laminates  by  Sun,  Achenbach 
and  Herrmann  1S18|  to  analyze  plane  harmonic  wave  I'ropajiat  ion  in  elastic 
laminates.  This  new  version  incorporates  h i,cher-order  thickness  variations 
in  the  displacement  functions  and  includes  restrictions  on  both  displacement 
and  stress  at  tlie  laminate  interfaces.  To  assess  the  potential  of  tlie 
second-oi'Jer  microstructure  tlieory  for  accurate  model  inj;  of  mechanical 
processes  in  laminates,  dispersion  curves  and  especially  mode  shape  data 
for  botli  displacements  and  stresses  were  obtained  and  compared  to  correspond- 
ini;  exact  solutions,  I'he  comparisons  indicate  that  while  dispersion  curves 
m;i>  be  nearly  identical,  extremely  significant  differences  may  be  observed 
in  the  mode  shape.  it  provides  a strong  evidence  tlwit  dis|’ersion  comparison 
alone  cannot  constitute  a valid  criterion  for  the  assessment  of  a microstruc- 
ture theory. 

Illavacek  (IILS]  used  an  effective  stiffness  theory  for  fiber-reinforced 
composites  with  hexagonal  packing  of  unidirectional  fibers.  In  the  theory 
a typical  hexagonal  fiber  and  matrix  layout  is  replaced  by  circular  cylinders 
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po^so-^s  1 Hi',  tho  s.iiKo  voliinie  of  materials,  lie  tl\en  illustrated  the  essence  of 
this  thoor>'  by  plane  harmonic  waves.  Propagation  in  directions  parallel  and 
normal  to  the  fibers  were  consiilered.  A correspondence  between  tlie  elastic 
constants  ot  an  eipiivalent  transversely  isotropic  medium  and  the  parameters 
of  this  effective  stiffness  model  was  shown  using  various  wave  speeds  at 
in  1 slung  wave  numbers. 

Illav.icek  (lllt'l  also  presented  an  effective  stiffness  model  for  isotropic 
two  phase  el.istic  composites  containing  spherical  inclusions.  His  model 
consists  of  two  concentric  spheres,  the  inner  one  represents  an  inclusion, 
tile  outer  one  the  matrix.  He  assumed  .1  continuous  field  ol'  gross  displace- 
ments ir  which  was  assigned  to  the  interior  of  the  smaller  sphere  and  to 
the  outer  surface  of  the  outer  sphere,  using  a linear  expansion  about  the 
common  center.  In  addition,  he  assumed  a local  homogeneous  deformation  of 
the  matrix  described  by  a continuous  tensor  field.  The  local  displacements 
in  the  matrix  were  then  fixed  by  the  condition  of  continuity  at  the  inter- 
f.ice,  and  bv  assuming  a linear  dependence  on  the  radius.  In  this  way.  he 
arrived  at  an  elastic  potential  density  as  an  expression  in  terms  of  the 
loc.il  deformation  and  the  gradient  of  the  gross  displacements.  Pispersion 
curves  for  harmonic  itlane  waves  were  jn'esented.  The  limit  phase  velocities 
for  vanishing  wave  numbers  yielded  the  effective  moduli  which,  for  a tungsten 
c.irb i de-coba 1 t alloy,  agree  very  well  with  Hashin  and  Shtrikman's  calculations 
IH.-i'. 


11 . Interact  ing  Cont  inuum  Iheories  and  Theories  of  Mixtures 

Bedford  and  Stern  |Blll  proposed  a one-dimensional  mixture  theory  for 
.1  lavered  medium  wherein  the  interaction  parameters  were  determined  on  the 
basis  iH'  results  of  certain  simple  puasi -stat  ic  problems.  I'he  novel  feature 
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of  thoir  theory  is  tluit  e;ich  eonstitiient  is  permitted  to  have  an  individual 

motion,  with  meelianical  coupling  between  tlie  motions  introduced  to  model  j 

constituent  interactions  in  the  actual  composite.  They  emjiloyed  this  theory 

to  study  the  time-harmonic  waves  propa>;at  inj;  parallel  to  the  interfaces  of  a ' 

laminated  composite.  I'lie  dispersion  curve  so  obtained  shows  good  ajireement 

i 

witli  the  exact  solutions. 

Stern  and  Bedford  IS1(<|  developed  a three-dimensional  model  leased  upon 
the  mixture  theory.  Results  were  applied  to  plane  harmoiiic  waves  proji.pcat  iiij; 
in  a layered  medium.  I'ispersion  characteristics  based  on  this  theory  were 
eompared  with  "effective  stiffenss"  tl\eory  and  exact  solutions.  It  was  found 
that  neither  the  mixture  theory  nor  the  "effective  stiffness"  theory  could 
descrilu'  adequately  the  material  response  over  a wide  ranj;e  of  wavelenjjths . 

(Inly  for  disturbances  propajjat  ing  in  the  direction  of  the  layerinj;,  this 
theory  reflects  faii’l>'  accurately  tlie  proper  kind  of  dispersive  behavior 
and  hence  max'  be  reyarded  as  a yood  model  for  such  analyses. 

Bedford  and  Stern  [Bl.'I  proposed  a mu  1 1 i -eont  i nuum  mixture  theory 
based  on  the  development  of  the  equations  by  means  of  two  coupled 
media  in  which  the  fiber  directions  are  natural  axes  of  svmmetry  for  the 
fiber-reinforced  material  and  the  transverse  isotropy  also  corresponds  to 
laminated  material  where  the  directions  normal  to  the  layers  are  the  natural 
symmetry  axes.  The  results  are  the  basic  expiations  for  the  considered  media. 

Martin  [Ml)  used  this  mixture  theory  to  predict  the  phase  velocity  of  pl.ine 
waves  in  a fiber- reinforced  elastic  solid,  lie  found  the  jihase  velocity  as 
a function  of  the  wave  number  and  as  a function  of  the  direction  of  prx>paya- 
tion  relative  to  the  fiber  directiixn.  /^^so,  the  particle  motion  of  each  of 
the  constituents  was  solved  for  at  various  wave  numbers.  It  was  found  that 
there  are  fxnir  possible  propayatiixn  velocities  in  a binary  mixture  of  tliis 
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tNpp.  and  tliat  two  nt'  tlic-^c  waves  correspond  to  nearlv  shear  iiu>t  i on  wliile 
tiie  rcmaininj>  two  correspond  to  nearly  compressive  motion. 

Ilej;emier  and  Nayfch  [I17|  developed  a continuum  theory  for  wave  propaga- 
tion normal  to  the  layering  of  laminated  composites  based  on  an  asvmptotic 
scheme.  Hie  hierachy  of  models  are  defined  by  the  order  of  truncation  of  the 
asimptotic  sequence.  Rased  on  tlie  calculation  of  the  lowest-order  dispersive 
model,  accuracy  superior  to  several  existing  theories  was  observed.  They 
then  cast  the  lowest-order  theory  in  a standard  mixture  form.  For  harmonic 
w.ives  propagating  normal  to  tlie  laminates,  the  jihase  velocity  spectra  of  this 
mixture  tlieory  exhibit  good  correlation  with  exact  results  for  the  first  mode 
to  tile  cutoff  freipiency.  The  sim|ilified  theory  in  binary  mixture  form  has 
practical  applications  in  view  of  its  sim|ilicity. 

Ilegcmier,  Gurtman,  .Vayfeh  and  Rache  [HS.HHJ  discussed  the  waveguide 
t>po  propagation  both  in  layered  and  fiber-reinforced  composites.  Based 
upon  asymptotic  expansions,  again  tliey  developed  a simplified  first-order 
mixture  theory  of  which  the  momentum  interaction  term  is  identical  in  form 
witii  Bedford's  hut  the  constitutive  equations  for  the  stress  in  each  con- 
stituent involve  the  strains  in  both  constituents.  They  emploied  this  first- 
order  mixture  tlieory  to  study  time-harmonic  wave  propagation  in  both  waveguide 
type  laminated  and  fiber- rei nforced  composites.  Dispersion  relations  were 
ol'tained.  .\  comparison  of  the  approximate  and  the  exact  phase  velocity  data 
indicates  that  the  theory  provides  good  first  mode  agreement.  Additional 
accurac)  may  he  obtained  by  taking  more  terms  in  the  asymptotic  expansions 
whicli  result  in  a more  com|ilicated  li i gher-order  mixture  tlieory. 
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Hcgemier  and  Bachc  (HIO]  further  extended  the  first-order  interacting 
continuum  theory  for  wave  propagation  in  one-dimensional  elastic  laminated 
composites  to  the  general  two-dimensional  clastic  case.  The  phase  velocity 
spectrum  of  the  general  two-dimensional  theory  was  investigated  for  one- 
dimensional  harmonic  wave  propagation  at  various  propagation  angles  with 
respect  to  the  laminates.  Comparison  with  the  first  three  modes  of  the 
exact  theory  gives  excellent  agreement.  It  is  seen  that  this  theory  exiiibits 
somewhat  better  correlation  than  the  effective  stiffness  theory,  especially 
at  propagation  angles  close  to  the  normal. 

Nayfeh  and  Curtmun  [NT]  employed  the  first-order  binar\’  mixture  theorv 
to  study  harmonic  shear  wave  motion  in  laminated  composites.  Both  trans- 
versely (SV)  and  horizontally  (,SII)  polarized  waves  were  considered.  Dispersion 
relations  were  derived  and  compared  with  the  results  from  the  exact  and 
other  approximate  theories.  Comparison  indicates  that  the  agreements  with 
the  exact  tiieory  for  both  tyjies  of  shear  waves  are  good,  esj^ecially  for  the 
SM  type  wave,  for  quite  a wide  range  of  wave  numbers,  and  that  the  superiority 
of  this  theory  to  effective  stiffness  is  obvious. 

Gurtman,  Nayfeh  and  Hegemier  [G15]  generalized  the  first-order  theory 
of  interacting  continuum  to  study  the  two-dimensional  wave  propagation 
ill  comoosite  materials.  In  the  generalization,  the  material  effect, 

the  elastic-plastic  behavior  and  two-dimensional  lay-ups  (laminates)  were 
considered.  The  dispersion  relation  of  longitudinal  harmonic  waves  in  the 
elastic  composite  was  put  in  a numerical  dispersion  code.  The  predicted 
response  agrees  well  with  the  experimental  results. 

Gurtman,  et  al.  |G14l  further  applied  the  theory  of  interacting  con- 
tinuum to  wave  propagation  in  tliree-dimensional ly  reinforced  composites. 

The  [»roposed  model  is  an  extension  of  their  two-dimensional  laminate  and 
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ft  ho  f- It*  I II  t'ori'emt’iU  fi>nmi  lat  i oils  w i 1 1)  an  inclusion  of  the  effect-^  of  thoniio 
dynamio';,  finito  Joformations  and  nonlinoai'  constitutive  behavior.  A PlSl' 
code  was  develonod  to  study  tlte  proptijiat  i on  of  liarmonic  lonj;i  tud  i na  1 waves 
in  three-d  iniensiimal  quart  z-pl\ei'.o  I i c composite.  l>ispersion  data  oluainotl 
from  tile  1)!SI'  code  were  compared  with  e.\i>er imental  results.  Aijreement  of 
all  values  witliin  the  reyion  of  validity  (for  frequency  < 1 Mkzl  of  the 
theori  is  quite  yood . 

\.i\  fell  |\.>|  u<ed.  tlie  theor>’  of  continuum  leased  on  asimptotic  expansions 
to  investigate  t i me-lia  rmon  i c waves  propagat  i njt  normal  to  tlie  lax’ors  of  multi- 
laiered  periodic  media.  I'v  retaining  all  terms  in  the  asynqitotic  sequences, 
lie  oi'tained  the  disj'ersion  relation  which  lields  the  exact  phase  velocit>' 
siiectrum.  file  known  relations  for  the  homogeneous  and  tlie  hi  laminated  media 
arc  obtainable  as  special  cases  from  this  theory.  Also,  it  is  noted  that  the 
dispersion  of  the  com[>osife  increases  witli  eacli  additional  layer  in  eacli 
cell  of  the  ci'mposite. 

\.i>  loll  IN'' I proposed  .1  continuum  mixture  theorv  of  lie.it  conduction  in 
laminated  waveguides.  Theory  leads  to  simiile  governiiv’.  equations  for  the 
actual  composite  which  retain  the  integrit.v  of  the  diffusion  process  in  each 
constituent  but  allow  them  to  coe.xist  under  some  defined  interactions.  The 
resulting  equations  were  used  to  study  tlie  harmonic  temperature  pulse.  The 
I'esults  were  fouiKl  to  correlate  well  with  some  existing  exact  solutions. 

It  w.is  also  mentioned  before  that  ben-Atnoz  |Rdd,Bd.il  propo'ied  a con- 
tinuum theory  for  waveguide  type  propagation  of  stress  waves  in  fiber- 
reinforced  composites,  using  direct  asynijitotic  expansions.  based  on  asxmptotic 
expansion  technique,  he  1 bJl | also  developed  a continuum  theory  for  wave 
pr('p.i)’,.it  i on  in  tavered  ci'iiiposi  t es  . Waves  ]’rop.igat  i ng  normal  f i'  the  planes 
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I't'  I, imin.it  ion  were  cons  i ilereil . \ pain  id'  t'iher  aiul  matrix  laminae  were 
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solci'teii  as  tlie  t'luuiainctita  1 unit  of  the  iierioJic  iiiedium.  Starting  from 
tiiree-dimcnsional  eijuations  of  motion  with  order  of  magnitude  arguments, 
a comh inat ion  of  tliiokness  expansion  and  crystal  lattice  techniques  was 
used  to  derive  continuum  equations  (of  zero  order").  The  theorv  was  then 
used  to  stud)'  hanaonic  wave  [iropagation  in,  the  laminated  composite.  Tlie 
jili.ise  velocit)  of  wa\ es  piropagat  i ng  t ransversel)'  to  the  laminate  was 
calculated.  No  numerical  e.xampile  nor  dispersion  results  compiared  witii 
otlier  tlieories  were  reported. 

Ren-.\moz  |B2Sl  considered  l'\eat  conduction  in  laminated  compiosites 

b)'  use  of  tlie  asymptotic  expansion  technique,  li.xpanding  asymptotically  the 
microstructure  equations  in  terms  of  a small  parameter  given  h>'  the  ratio  of 
d i f fus  i V i t i es  of  tlie  two  const  i tuents , lie  again  developed  a continuum  model 
for  heat  conduction  in  which  microstructure  effects  apiiear  as  ;i  consequence 
of  the  fact  tliat  tlie  current  state  in  a laminated  medium  is  lii story-deiiendent . 

Such  a history-dependence  reflects  the  effects  of  microcondiict ion-occurring 
within  the  microstructure  on  the  micro-time  scale. 

11.(1  -'h'thods 

Based  on  perturbation  theor)'  and  on  the  statistical  consideration  assum- 
ing elastic  iiarameters  of  the  medium  to  be  subjected  to  small  random  fluctuat ions . 
Hudson  (HJ(1|  discussed  tlie  scattering  of  surface  waves  due  to  a randoml) 
inhomogeneous  medium,  lie  derived  exiiressions  for  the  attenuation  of  hove  and 
Rayleigh  waves  when  the  size  of  inhomogenci ty  is  small  compared  with  the 
wavelengths  of  incident  waves.  These  expressions  have  the  same  forms  as  the 
ktuiwn  expressions  for  the  v i scoe  1 ast  i c .if  t enuat  ion  . Comparing  the  two.  he 
fiiund  that  the  attenu.ifion  mech.inism  by  scattering  is  iust  like  a cert.iin  t>pe 
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of  V 1 scoo  I ast  1 c i ry  aiui  that  tho  ^^cat  rof  i iij;  is  more  offcwtivo  at  tti^th  t requeue  i os 

aiKl  viscoelastie  attenuation  is  more  effective  at  low  frequenci^ 

Based  on  the  assumption  that  geometric  dispersion  results  mainly  from 
tite  relativel\'  [periodic  arrangement  of  the  rein torcement  elements  in  the  matrix 
rather  titan  from  the  precise  shape  of  each  reinforcing  elemetit,  nrumheller  and 
Sutherland  iniO]  developed  a lattice  model  which  ignores  the  shape  of  the 
reinforcing  elements  hut  preserves  tlteir  periodic  it,v.  In  the  aitplication  ot 
this  lattice  model,  the  composite  was  initially  treated  as  a nondi spers i ve 
homogeneous  mixture.  I'he  ettective  or  average  iM'operties  ot  the  mixtuie  wete 

iletermined  l\v  eithei-  >.t  ead\ -wave  analtsis  or  apitropriate  exiter  i ment  s . A 
lattice  was  then  formed  hy  redi  st  r i itut  i ng  tlie  mass  witliin  the  mixture  tit 
form  a iterioilic  structure  of  laminates.  This  mass  redistribution  was  carried 
out  in  a manner  wliich  yields  a lattice  witit  theoretical  dispersive  charac- 
teristics tliat  match  the  measured  dispersive  characteristic  of  the  composite. 
Hence,  for  liarmonic  wave  [tropagat  i on . tliis  model  describes  ver.v  well  the 
behavior  of  actual  engineering  composites. 

Nelson  and  Navi  [Nil!  also  used  a lattice  model  to  studv  the  harmonic 
waves  of  plane  strain  propag.iting  in  conqiosite  materi.ils.  Tlie  medium  was 
considered  as  a periodic  assemblage  of  identical  cells.  Using  a finite 
element  technique,  tlie  dynamic  beliavior  of  one  reference  cell  was  represented 
by  a number  of  generalized  displacements.  Thus,  the  continuum  periodic 
structure  of  the  medium  was  idealized  b\'  a discrete  lattice-t>pe  structure 
and  the  problem  was  reduts  il  to  an  algebraic  eigenvalue  problem.  They  deter- 
mined the  eigenvalues  for  several  examples.  The  frequency  spectrum  for  a 
fiber-reinforced  composite  medium  was  compared  with  the  results  of  effective 
modulus  .and  effective  stiffness  theories.  I'h  i s theorx’  predicts  correctl>-  the 
vanishing  of  group  \elocit>  tor  ha  1 f -wave  1 engt h equal  to  the  interfiber 
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li  1 t iiMi'v . riu'  cl'loct  no  •.tiHnoss  tiu'orv  yiolils  ;i  oiitot't'  at  a smallor  wavo 
numbor  win  lo  ttio  ot'footivo  iiuhIuIio.  (I\oor\  fair,  to  proJiot  any  J i spot's  i on . 

Ilonoo,  ( lio  I'osiilis.  wboii  oomparovi  with  av.iilal'lo  aiial>tioal  .aiul  nnmoi'ioal 
rosiilts,  •;|u'w  that  t ho  mothovl  y.ivo'.  an  aooiifato  niodol  whiott  i •.  ahio  to  nioilol 
)'oi'io>lio  stt'ootufos  of  wavolon^th  smallof  than  t lio  lattioo  oonstant. 

Ilabip  |lll|  applioil  I ho  iii  n I'omoipli  i o tl\ooi'\  v>l'  ov'ntimia  t ovaniino  t ho 
timo  hai'inonio  wavos  pt'oipay.at  i ny.  in  a oomposito  with  olastio 

inioro  i no  Ins  ions,  Itv  analo);\  wi‘h  t ho  rosults  »it'  Iriny.on  aiul  Siilnibi  1 1 I .mil 
ot  Miiiillin  |NM|.  ho  oblainoil  ovplioit  il  i spot",  i on  I'olation'.,  i;ovoi'ninK  tho  lowost 
ainl  no\i  hi);hoi'  iiioilo'-  ot  piopa};.i  I i on  of  pl.ino  I oi;  i t ml  i na  I w.i\o-.  in  .in 
iiiil'oiiiuloil  ol.istii  oomposito  sol  111,  in  toniis  ot'  t lio  rol.it  ni'  proport  1 1".  ot  tho 
oonstituont  iii.it  o r i a I s . I'lio  oor  rospoiul  i iij;  r.itio  of  i;r'’"l’  volooiti  to  j'liaso 

Ki'looiti  w.is  liKowiso  ova  I iiat  I'll . Rosults  valivl  tor  a spooi.il  oaso  woro 

ovhihitoil  in  v;r.iphs. 

I'.ioy.lor  1 .' 1 I ooiisiiloroii  tho  mo. in  w.ivos  i n 1 ,iiii  i n.i  t oil  r.iiiilom  ooiii|  .v.  1 1 os  . 

Ill  iisoil  .1  I iiip  I o \orsioii  ot  tho  olt'ootivo  st  i t t'ni'ss  thoori  for  l.imiiiat  Oil 

iiioili.i  111  iloriviny,  -.t  ooh.i-.t  i o il  i sp  1 .loomont  oipi.it  ion  of  motion.  App  1 i oat  i on  of 

tho  port  iirl'.it  1 on  proooiluro  ot  Kollor  [KJl.  moil  i I i oil  hi  in  t roiliio  i n>;  tho  spoo 

I r.i  I iloiisitios  ot  tho  r.iiiilom  ooo  I I i o i out  . t lion  loil  to  ilot  ormi  ii  i st  i o oipi.it  ions 

tor  moan  wavo  I'l'op.iypit  ion.  lio  >;avo  spooi.il  .ittont  ion  tv'  unoouploil  moili's  ot 
moan  w.ivos  prop.iy.at  i nr.  pori'oiul  i on  1 .i  r or  par.iliol  to  tho  ilirootion  ot  tho 
l.ivoriii);.  Ho  o-.t  ah  1 i slioil  ilisporsion  rol.it  ions  tor  thoso  piano  h.irmonio  w.ivo 
t ra  1 11'. . 

Kniiiihaiisl  |li''l  .i  I '.o  ilisonssoil  tho  ramloiiinos  s .nul  .ivor.iy.o  or  mo  in  w.ivo 
propaj;.!  t 1 on  in  i iihomoronooiis  moili.i.  Ho  omplovoil  tho  mothovls  usisl  in  .olul 
st.it  o pin  SI  O',  tor  '.iil'st  .nit  1 .1 1 1 ' vl  i '-oriloroil  .illovs  to  stmlv  tho  subst  .iiit  i .i  1 1 v 

il  I so  I'llo  roil  coiiipo'.  1 1 o . lor  .11111  oompo'.  1 1 1's  , ho  sub',  t 1 1 iit  oil  loiirior  1 loipiot 

iiiollioils  for  lonrior  pi. 1110  w.ivo  iiiothoils.  In  t.iot,  his  ,i|'pro.ioh  is  .1  portnrb.i 


tioii  iiii't  lu>il  iii.ikiiii;  iisi'  nt  1 Kuiui'l  '-vilut  ii'ns  ,1 1 1'l'.nK  lU'vi- 1 o|H'iI  (\>i'  .m  .iv  i' ij 
movlium.  Irv'iu  till-'  th^’oi'v  lu'  pri'tl  i v'l  I'll : (M  thfro  \\  i 1 I lu'  Wf  1 I »li’liiu‘il 
I'U'qui't  w.ivi"’,  I'ut  i ti  ^fiuT.il  t(u‘v  will  lu'  il.iiiipi'il  out  in  tinio; 

(-M  till’  low  t'l'i'qiioiK'N  nvoiMi'.o  wjvos  will  In'  J,im|n'il  t lio  {.■>1  tlio  ills 

I'orsioii  I'lii'fos  .iro  slu  t t Oil  t I'om  thoso  ol  t lu'  porioilio  moiliuiii;  (1)  nwnli'f.it 
hoiiioi’.onoous  rniulomiu'Ss  iloos  not  ilost  I'ln  rioi|iiot  proportios  in  t lu'  avoi'.ip.o 
wavo's;  (M  somo  possil'U  uooiU'il  i n to  iiiia  t i on  . siuh  as  t ho  I'luot  iiat  ions  from 
t ho  a\orav;o  tioUl  and  t Ito  soattorod  liold  and  t ho  looal  i n to  riiia  t i on  , is  lost 
In  t lu’  pi'oposi'd  proi'i'dii ro  . 

■'lotos  lM-1  '-tiiiliod  hariiionio  wi\o  propagation  in  d i 'm' rdo  rod  oomj'i".  1 1 os  . 
Ilo  ohtaiiu'd  a forimilalion  tl\at  is  to  ho  '-atisfiod  hv  t ho  '.tatistioal  avorap.o 
t'iold  i|nantitit's  in  a statr-tioal  samplo  ot'  hot  o ropi'noons , linoarls  s'lastu' 
'.oliils,  \ low  t'roi|uono\  lono.  watt'lonrlh  thoors  thon  was  ovtraotod  I rom  t ho 
vu'noral  t'orinn  1 a t i on . 1 ho  prodiotions  ot  this  tlu'ors  oaii  ho  p,i\on  a piirols 

dot  I'l'iiii  n i st  i o i nt  orprot  at  I on  . Somo  spooial  oa'.os  of  his  thoors  ri'diioos  to 
an  ot  ti'oiivi'  inodiiliis  tiu'ors.  llowovor,  hs  in'tainini',  oorri'ot  ion  torms,  it  was 
shown  that  I'lastio  wato  propagation  will  alwa's  oxhihit  both  disporsion  and 
dooas  ovor  larp.o  onoiu'.li  propagation  d i t anoi's  . 

Itoso  and  Mai  pill  al'-o  omplovod  t ho  statistioal  approaoh  to  disoiiss 
t lu'  1 onp,  i t lid  I na  I shear  watos  in  a lihor  rointorood  oomposito.  IVavos  are 
harmonio,  tihi'i'.  are  raiidoiiilt  d i sf  ri  hut  od  in  paralli'l.  I'lu'  ooiiipos  1 1 »'  is 
t at  I '' t 1 oa  1 I N iinitorm.  I'lios  ohtaiiu'd  t lu'  phase  velooit''  .ind  daiiipiiii',  ot  the 
aii'rapo  watos  as  t'linot  i on ol  the  -.1  at  i s|  i o.i  1 and  the  nn'olian  i oa  1 j'a  ramol  o rs  , 
llio  tlioort  li'ads  to  lla'dnn  and  Rosen's  foriinila  ino|  tor  the  .i\ial  shear 
modulus  It’  the  oor  ro  1 a t i oils  in  the  position'-  ol  fihors  are  ipnorod.  Iho 
oorrol.ilion  tonus  h.ito  i s i r.n  i f i oant  otiool  on  I hi'  d.imp  i nr  proport  v of  Iho 
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V oiiipo'- 1 t o , O'.pi'v  I a 1 1 1 ,it  hirh  I ri'iiuono  1 1".  and  voiu  out  ra  I i on  . 
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the  correlation  terms  is  to  increase  the  velocity  and  decrease  the  specific 
dampinj;  capacity. 

Bose  and  Mai  [B42|  a>;ain  discussed  the  time-harmonic  plane  waves  in  a 
fiber-reinforced  composite  by  statistical  approach.  The  compos i t e cons  i st  s of  a 
homogeneous  isotropic  medium  containing  long  parallel  randomly  distributotl 
circle  fibers  of  identical  properties.  They  discussed  the  general  case  of 
time-harmonic  plane  waves  propagating  perpendicular  to  the  fibers.  Ihe  m.iin 
result  was  to  obtain  tlie  phase  velocity  and  the  damping  of  the  waves  of  both 
corapressional  and  shear  types  into  which  the  average  waves  .ire  to  be  se|'arated. 
In  the  case  ol  thin  fibers  Icompared  to  wavelength),  the  results  concerning 
the  transverse  bulk  modulus  and  transverse  rigidity  were  compared  with  those 
of  llashln  and  Rosen.  Numerical  calculations  were  made  for  boron  filters  in 
an  aluminum  matrix.  It  was  found  again  that  correlation  terms  or  tlie 
correlations  in  the  positions  of  fibers  have  a significant  effect  on  the 
damping  property  of  the  composite,  es)H'cially  at  high  frequencies  and  concen- 
t rat  ions . 

As  mentioned  earlier,  Chao  and  l.ee  [C,^|  developed  a discrete  continuum 
theory  for  periodically  layered  composite  materials.  'l'he\  presented  the 
governing  field  equations  incorporating  interface  continuity  conditions  in 
the  form  of  a system  of  di  fferent  i al -d  i f ference  equations.  Tliey  ai’plied  this 
model  to  discuss  the  propagation  of  I’lane  harmonic  waves  in  an  unbounded 
layered  medium.  Thickness  twist  vibrations  were  studied.  Numerical  results 
were  compared  with  those  of  exact  solution  (harmonic  waves).  It  was  seen  th.it 
tlie  results  agree  witli  exact  results  quite  closely.  Agreements  are  lu'tter 
than  effective  stiffness  tlieory  in  general  as  wave  number  graduall\  increases. 


It  was  also  mentioned  before  that  Ben-Amo:  extended  Miiullin’s 

theory  of  elasticity  with  microstructure  to  a heterogeneous  medium  consisting 
of  inclusions  of  .irbitrary  geometry  embedded  in  a m.itrix  materi.il.  lie  .ipi'lied 
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this  thoory  to  examino  harmonic  plane  wave  propagating  in  such  a medium.  The 

(J 

dispersion  curves  for  a cubic  array  of  spherical  particles  obtained  by  this 
theory  are  ratlter  similar  to  the  curves  sketched  by  Mindlin  [M4]  for  tlie  lowest 
acoustic  branches.  Certainly,  the  theory  developed  by  Ben-Amoz  is  an  effec- 
tive tool  in  treating  wave  propagation  in  media  lacking  periodicity  such  as 
inclusions  of  arbitrary  geometry  embedded  in  a matrix  material. 
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III.  TRANS  1[:NT  WAVHS 

In  this  section,  transient  analyses  of  waves  propagating  in  composite 
mate  Is  based  on  Tourier  synthesis  are  first  discussed.  Then,  in  addition, 
other  techniques  for  obtaining  the  solutions  to  transient  problems  are 
mentioned. 

1 1 1 . 1 Tourier  Synthesis 

In  discussing  a pressure  pulse  propagating  in  lieterogeneous  materials, 
Osten  [02]  used  a computer  code  called  TDRRY.  lie  e.xpresscd  the  applied 
pressure  pulse  at  t = 0 as  a Tourier  series  in  time.  The 

’ amplitude  of  each  component  in  the  Tourier  series  were  then  determitied 

numerically.  To  determine  the  pulse  shape  at  a later  time  t,  each  com- 
ponent of  the  initial  pulse  is  allowed  to  propagate  the  distance  ct,  where 
c = c(.x)  and  is  the  circular  frequency  of  that  Tourier  comi>onent.  The 

I differences  in  c thus  cause  the  pulse  to  disperse.  Osten  also  allowed 

I - -Yt 

tor  attenuation  ot  each  component  in  the  form  e where  y is  real  and 

, also  y = yC^).  Osten  applied  the  TURRY  code  to  the  sinusoidal  dispersion 

solutions  he  obtained  for  a composite  with  randomly  distributed  inclusions 

i of  spheres  and  cylinders.  He  compared  TURRY  code  predictions  with  the 

attenuation  curves  obtained  by  Bjork  [B391 , by  a series  of  2-U  finite 

j difference  code  calculation,  and  found  that  TDRRY  underpredict ed  the 

attenuation.  He  attributed  the  discrepancies  to  absence  of  nonlinear 
attenuation  effects  in  TDRRY  and  speculated  that  nonlinearity,  acting  in 
conjunction  with  geometric  dispersion,  may  tend  to  enhance  the  di sjiers i\ <.■ 
at  tenuat ion. 

j Peck  and  (lurtman  [P.3|  carried  out  Tourier  analysis  of  pulse  propagat- 


ing parallel  to  the  interfaces  of  a laminated  composite.  They  considered 
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a uniform  pressure  pulse  of  step- funct Lon  in  time  applied  to  the 
boundary  of  a half-space.  The  laminates  are  perpendicular  to  the  half- 
space boundary.  Starting  from  the  equations  of  elasticity,  they  applied 
Fourier  integral  transforms  over  time  t and  the  propagation  direction  ,x 
in  the  analysis,  hxact  solutions  were  obtained  in  the  form  of  an  infinite 
series  of  integrals,  each  of  which  is  the  contribution  to  the  transient 
response  from  a mode  of  sinusoidal  wave. 

Peck  and  Gurtman  then  evaluated  the  asymptotic  solution  in  integral 
form  by  using  the  saddle-point  technique  developed  by  Polk,  et  al . 11'2J. 

Of  particular  importance  is  a contribution  from  the  first  mode,  which 
dominates  the  long-time  solution.  This  contribution  is  known  as  the 
hcad-of-the-pulse  apjtroximat ion.  Tor  the  strain  in  the  direction  of 
propagation,  the  head-of-the-pulse  solution  takes  the  form 


where 


lUSl  = j - 


A.(n)dn  , 


C = tt-x/Cjjl/T  , T = iSc^^axl 


l40) 


(47) 


in  which  c^^  is  the  speed  of  wave  propagation,  a is  tiie  crucial  disper- 
sivity  parameter,  is  the  static  strain  under  the  same  load  and  A.  is 

the  well-known  Airy  function.  The  above  expressions  can  be  interpreted  as 
follows;  the  wave  is  roughly  a stop  pulse  arriving  at  = 0,  i.e.,  at 
t = x/Cy,  so  that  the  main  disturbance  is  propagating  at  speed  c^^.  The 
quantity  T has  tlie  dimension  of  time,  and  t - x/c 


is  t ime  after  arrival 
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of  the  wave.  Thus  ^ is  a nondimensional  time  of  arrival.  The  quantity 
T may  be  regarded  as  a characteristic  dispersion  time  of  the  pulse. 

liquation  (4b)  states  that  the  strain  is  uniform  in  the  direction 
of  propagation.  The  stress  in  the  x-direction  is  quite  non-uniform  because 
the  moduli  of  the  layers  are  different.  In  terms  of  the  average  stress 
over  the  cross  section,  a(x,t),  one  has 

o(x.t)  O^H(f;)  (48) 

where  is  the  applied  stress. 

Further  applications  of  the  head-of-the-pulse  approximation  will  be 
mentioned  later  when  it  is  employed  in  other  problems. 

Voelker  and  Achenbach  [VI]  carried  out  a similar  analysis  of  stress 
waves  in  a laminated  medium  generated  by  transverse  forces.  The  laminated 
medium  composed  of  alternating  layers  of  two  homogeneous  isotropic  elastic 
solids  is  suddenly  subjected  to  a spatially  uniform  distribution  of 
transverse  forces,  which  are  applied  in  a plane  normal  to  the  layering. 

'Hie  resulting  two-dimens ional  transient -wave  propagation  problem  is 
analyzed  by  means  of  modal  analysis.  The  normal  and  shear  stresses  at 
the  interfaces  arc  expressed  as  infinite  integrals  that  are  integrated 
for  not  too  large  values  of  time.  For  largo  values  of  time,  the  integrals 
are  estimated  by  the  method  of  stationary  phase.  The  predominant  contri- 
bution to  the  interface  shear  stress  comes  from  the  head-of-thc-pulse 
approximation.  I’he  normal  stress  at  the  interface,  which  is 
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composed  of  several  contributions,  is  oscillatory  (and  dies  out  at  t ' ), 
and  thus  the  interface  bonds  may  be  subjected  to  tensile  stress.  The  solu- 
tion obtained  should  be  useful  in  determining  the  validity  of  the  continuum 
theories. 

Balanis  [B5]  considered  the  transient  one-dimensional  wave  propagation 
in  a semi-infinite  periodic  composite  medium,  which  consists  of  periodic 
array  of  two  elastic  layers  perfectly  bonded  together.  Waves  are  generated 
by  a surface  velocity  input  and  travel  perpendicular  to  the  layers.  He 
applied  the  t-ourier  integral  transform  over  time  t in  the  analysis. 

Fourier  transformed  interface  velocities  and  stresses  are  analyzed.  An 
approximation  for  materials  with  similar  impedences  is  then  presented.  He 
found  that  for  situations  where  the  times  taken  by  the  longitudinal  wave 
to  traverse  the  layer  thicknesses  are  in  an  integer  ratio,  periodicity 
results  in  the  frequency  domain.  He  treated  such  a case  and  olitained 
the  interface  stressed  due  to  a square  wave  input.  The  results  of 
this  example  are  quite  interesting.  However,  the  approximation  scheme 
seems  to  suffer  from  lack  of  generality. 

1 1 1 . 2 F loquet  or  R loch  'n^eories 

Krumhansl  [K8]  analyzed  the  propagation  of  transient  waves  in 
periodic  composites  by  means  of  F'loquct  or  Bloch  tlieory.  He  uscxl  the  quasi - 
periodic  Floquet  waves  to  form  a basis  for  the  analysis  of  transient  problems 
and  expressed  the  solutions  of  the  problems  in  the  form  of  l•ourior-l•toquet 
scries.  A brief  summary  of  this  theory  for  one-dimensional  systems  is  as 


fol lows: 


57 


Consider  the  Tloquet  wave  solutions  t'or  displacement 

= v^(.x;q)  exp  Umx)  (.-Ji)) 

with  time  variation  prescribed  by  the  additional  factor  expitiiat),  where 

^ . (50) 

q is  the  wave  number,  and  v represents  the  vth  mode. 

The  function  v^(x;q)  is  periodic  witli  the  period  of  tlie  composite  con- 

fijiuration.  The  solutions  are  defined  in  V denumerable  rej;ions  of  u).  the 

pass  band,  and  for  -ir/a  q < ^/a  where  a is  the  width  of  a unit  cell. 

Ntultiplyinj;  constant  factors  can  be  selected  for  the  u functions  to 

generate  an  orthonormal  set  satisfying  the  relation: 

CO 

p(x')u*(.x;q)u^,  (x.q)dx  = 6^^,6(q-q'')  (51) 


rhe  band  orthogonality  is  given  by  the  Kronecker  delta  and  the  q space 
orthogonality  by  the  IHrac  6-function.  Using  tliis  relatioi\,  solution  for 
propagation  of  waves  without  applied  forces  can  be  represented  in  the  form: 
/ n/a 


u(x;t)  = I 
V 


■n/a 


+ + c^,(-,q)e 


- iu)\,t 


u^lx;q)dq 


(52) 


with  the  coefficients  c^(±,q)  determined  by  application  of  orthogonal ity : 


iu\,t 


- iu)t 


CO 

p(x)u*(x;q)u(x;  t )d.\ . 

-OO 


iu^t 


- iu>^,t 


p(x)u^(x;q)u(.c;t  )dx 


(S.i) 


equation  (52)  thus  provides  the  solutioit  for  an  initial  value  problem  wIumi 
the  initial  displacement  and  velocity  are  prescribevi.  I'or  initial  disturbance 
in  .1  localized  region,  this  met  liod  enables  one  to  maniinilate  the  solution 


r 
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into  the  form  of  ;i  st.it ionary  phase  integral  for  which  a head-of-the-wave 
analysis  can  be  carried  out.  It  is  also  noted  that  by  separation  of  solu- 
tions into  components  even  and  odd  in  x,  half-space  jiroblems  can  be 
t rented . 

Krumhansl  and  bee  lK10|  furtlier  employed  the  l-our i er-l-loquet  methods 
to  transient  elastic  waves  in  periodic  composites.  I'hey  found  tlvat  tlie 
he;id-of-tlu'-wave  solutions  in  the  far  field  are  dominated  by  the  iiuiximum 
itroup  velocity  contributions  from  each  frequency  band.  for  an  infinite 
periodically  layercs.!  nK\lium  wliich  is  initially  at  rest  and  subjected  to  a 
momentum  impulse  at  x = 0 lit  would  be  convenient  to  arranjte  x = 0 at 
symmetry  plane  of  tive  periodic  compositel,  the  velocity  response  at  far 
field  is  jtiven  in  closed  fonii  by  a set  of  Airy  functions.  fhe  ai-ijiearance 
of  Airy  functions  in  the  far  field  rosjxvnse  is  just  like  wlvat  I'eck  and 
t'.urtman  obtained  in  the  Ivead-of-the-pulse  approximation  for  an  analojjous 
problem.  In  additioiv,  krumhansl  and  bee  found  the  asymptotic  behavior  of 
a pulse  initially  uniform  over  a plane  with  the  followini;  features:  (a)  the 
asvnnptotic  peak  amplitude  decreases  as  , lb)  the  pulse  spreads  out 

and  becomes  less  steepei-  with  rise  time  proportional  to  Ic)  peak 

stress  also  falls  off  at  large  x,  tlvough  it  is  complicated  to  calculate. 
They  also  discussed  the  motion  of  a periodic  layered  composite  vinder  applied 
forces.  The  discussion  includes  both  a formal  representat ioiv  of  the 
(Ireen's  function  and  the  treatment  of  short  pulse  pressure  loading 
on  a surface. 

A four ier- Bloch  technique  for  superposing  the  finite-element  sinusoidal 
Bloch  wave  modes  to  generate  transient  solut ions  to  J-D  buindary  value 
problems  was  also  presented  by  I’eck,  et  al.  |l'4l.  lor  a step-pressure  load 
oti  the  sui’face  of  the  conqiosite,  the  responses  are  calculateil  at  small  aiul 
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large  propagation  distances.  Five  modes  give  good  convergence  at  small 
propagation  distance,  about  1/2  spacing  of  the  longitudinal  layers, 
although  even  more  modes  would  be  desirable.  However,  at  the  larger 
propagation  distance  (about  twice  the  spacing  of  the  longitudinal  layers) 
one  mode  only  provides  a good  approximation  to  the  response.  This  depen- 
dence of  number  of  modes  required  on  the  propagation  distance  provides 
some  insight  into  the  regions  of  validity  of  one-mode  or  two-mode  approxi- 
mations that  are  commonly  obtained  with  the  continuum  theories.  The 
transient  solutions  also  showed  that  the  layering  in  the  propagation 
direction  changes  the  micromechanical  shear  stresses  at  the  edge  of  the 
longitudinal  layer  by  a factor  of  two,  even  though  the  layering  has  little 
effect  on  the  overall  dispersion. 

Kohn  [K6]  considered  the  propagation  of  low-frequency  elastic  disturb- 
ances in  an  infinite  one-dimensional  composite  which  is  of  periodically 
varying  density  p(x)  and  stiffness  n(x) , with  spatial  periodicity  a. 

By  using  Floquet  theory  and  low  frequency  expansions,  he  found  that  in  the 
limit  of  low  frequencies  the  displacement  u(x,t)  in  the  composite  can  be 
written  in  the  form  of  a differential  operator  acting  on  a slowly  varying 
envelope  function  U(x,t)  as 

u(x,t)  = [l+v^(x)3/a^  + ...]U(x,t)  (54) 

U(x,t)  itself  describes  the  overall  long  wavelength  displacement  field.  It 
satisfies  a wave  equation  with  constant  (independent  of  x)  coefficients 
obtainable  from  the  dispersion  relation  w =•  u)(k)  of  the  lowest  band  of 
eigenmodes,  i.e. 


- c^  aVax^  - e aVax"  + ...)U(x,t)  = o 


l.'iS) 


()t1 


whore  c is  ttio  lonj;  wavelength  sonnU  velocity  of  the  composite  aiul  3 is 
the  Itlt-orJer  exp;msion  coefficient  of  the  square  of  the  frequency  of  the 
lowest  hand  of  eigenmodes 

u)‘lk)  = c‘k-  - Bk"  ^ 


V|(\)  in  (.51)  is  given  by 


Vjl-v)  = 


.Au- 


U/  - 


where 


I -1. 


n-1  = - 1 t1  lx)dx 


1-’ 


ISS) 


Infomation  about  the  local  str.iin  on  the  microscale  of  the  composite  l;uiiir\ae 
is  contained  in  tlte  function  v^lx).  kohn  then  appli(.\l  tliis  approach  to  study 
tlie  pulse  propagation  in  a lia  1 f-sjvice . Appropriate  tlreen's  functions  are 
constructed  in  terms  of  Airy  functions.  Witli  tliese  tlreen's  functions,  t lie 
solutioit  of  the  initial  value  prolilem  for  a composite  li.ilf-space  with  free 
boundary  is  explicitly  obtaincil.  It  is  noted  that  this  long  wavelength  tltev'vy 
correctly  descrilu's  the  head  of  .i  resulting  iwilse,  even  if  the  ju'escribed 
initial  pulses  have  step-function  character. 

kohn  |k?l  extended  the  previous  work  to  three-dimensional  composite 
material.  In  tliis  case  the  displacement  solutions  to  equations  of  motion  of 
periodic  elastic  medium  are  expressed  as  a vector  in  tlve  form  of  a differential 
i>perator  acting  on  .i  vector  t'unction  whicli  ilescribes  the  mean  displacement 
I’f  iMch  cell.  Ag.iin,  loc.il  '-ir.nn  infi'rmation  c.in  be  ol'laiiuwl  from  the  sv>ln 
lion,  ibie  iif  the  .ipp  1 i c.il  i ons  otthis  metluHl  is  the  structure  of  t lie  lic.wl-ot- 
a-\iulse  propagating  in  an  arbitral')’  direction.  However,  application  ol 
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solutions  to  fiber- reinforced  materials  is  conceptual  rather  than  explicit. 
His  approach  may  be  regarded  as  a method  of  microscopic  field,  although 
relation  of  his  work  to  other  dynamic  theory  was  not  given. 


111.5  Effective  Modulus  Theories 

Aboudi  and  Weitsraan  [Ad]  discussed  the  two-dimensional  plane  problem 
of  an  impact  on  a half-space  reinforced  by  parallel  elastic  fibers  witli  a 
pulse  of  finite  duration.  They  assumed  that  the  rigidity  of  tlie  fibers  is 
greater  than  that  of  the  matrix.  Based  on  effective  modulus  theory,  as 
described  previously,  they  expressed  the  constitutive  equations  in  terms  of 
the  reinforcement  ratio  and  angles  of  inclination  of  tlie  fibers.  Solutions 
were  obtained  by  a finite  difference  sclieme.  Numerical  results  were  presented 
for  the  case  of  normal  impact  on  a lialf-space.  Comparison  with  ♦•he  isotropic 
case  showed  good  accuracy  of  the  solutions  (.with  maximum  error  less  than 
1/2  percent).  The  investigation  indicates  that  the  vertical  displacement 
iunplitudes  decrease  monotonical ly  with  fiber  orientation.  However,  the 
horizontal  displacement  amplitudes  are  shown  to  vary  with  position  with  maximum 

and  minimum  displacements  less  predictable  as  a function  of  fiber  orientation. 

Aboudi  and  Weitsman  [A5]  extended  their  approach  to  study  the  problem 
of  an  impacted  fiber- reinforced  viscoelastic  half-space.  The  matrix, 
they  considered,  is  an  arbitrary  linear  viscoelastic  homogeneous  and 
isotropic  half-space.  The  embedded  oblique  fibers  are  linear  elastic, 
arranged  in  parallel  order  and  randomly  dispersed.  They  adopted  the 
effective  modulus  theory  to  study  the  dynamic  response  of  the  composite 
half-space  to  a time  dependent  surface  load.  Again,  a numerical 
procedure  of  finite-difference  type  was  applied  to  the  set  of  equations 
resulting  from  the  theory.  Results  were  worked  out  for  an  elastic- 
elastic  half-siKico  and  a viscoelastic  epoxy-elastic  glass  fiber  luilf- 
space.  Stal’ility,  convergence  and  accuracy  of  the  numerical  procedure 
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wore  cons uK'vcd . Numerical  results  for  surface  displacements  and  disjilace- 
ments  at  some  distances  lieneath  the  surface  were  obtainc'd  and  shown 
jirajihical  ly . Since  eft'ective  modulus  theory  only  j^ives  weighted  displace- 
me!\ts  and  stresses,  it  is  highl\’  ])ossihle  that  the  effective  stiffness 
theory  or  interacting  continuum  theor>’,  which  accounts  for  microstructure 
information,  might  give  more  reliable  numerical  results. 

bason  [b.l  1 analyzed  the  propagation  of  waves  resulting  from  the  sudden 
aiiplication  of  a force  to  tite  surface  of  a cylindrical  cavit>'  in  an  infinite 
fiber-reinforced  solid.  The  fibers  were  assumed  circumferential  to  the 
cavity  and  uniforml\’  distributed  so  that  the  material  might  be  regarded  as 
a transverse  isotropic  solid.  He  used  an  effective  modulus  model  to  represent 
this  solid.  The  stress  distribution  within  tlie  solid  then  was  determined 
numerical  !>•.  He  found  that  tlte  maximum  of  the  tangential  stress,  , at  the 
ca\ity  increases  from  a value  which  is  sliglitlv  greater  than  the  suddenly 
apj'lied  force  Q for  the  isotropic  solid  to  a value  of  several  times  of  Q 
for  a highly  anisotropic  solid.  He  also  found  that  the  time  at  which  this 
maximum  is  achieved  decrecises  from  the  time  for  isotroj'iic  solid  to  a fraction 
of  that  time  for  highl_\-  anisotropic  solid. 

Sun,  beng  and  Kol)  |SJS,S.%1|  investigated  the  shear  wave  propagation 
it!  a nonlinear  clastic  fiber-reinforced  composite.  I'he  comnosite  was 
modeled  b\'  a medium  consisting  of  tliin  nonlinear  ( i sot  rojn  c ) m.itrix  layers 
alternating  with  the  effective  linearly  elastic  fibrous  Layers.  Using  the 
effective  modulus  theory,  they  derived  a set  of  three-dimension;il  mnil  inear 
constitutive  etjiiations.  The  model  was  emjiloyed  to  investigate  the  wave 
front  propagation,  the  sfabilit)’  of  the  shock  front,  growth  ainl  decay  of 


the  shock  wave  and  the  distortion  of  initial  !>■  sinusoidal  wave.  Hue  fi>  the 


lu'itl  iiuMi  i ( y , ;i  stiulsoid.il  stuMf  w.tvi'  i'  il  i st  ort  oJ . At  sonio  ilislaiu'i'  a ilis- 


lontiiuuty  may  oiiiir  at  a critical  'imo.  Shock  front  may  form  at  siicli 
iliscont  imiity  from  an  initialU  hanmmic  wave.  It  is  notoii  that  shocks 
can  ho  fv>rmo»i  for  composite  with  hariloninjj  matri\  hut  cannot  ho  forimvl  for 
(l)af  with  softoniny  matrix. 
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1 n . l rffoctivo  Stiffness  Ihoorios 

Svo  IS-UI)  discussed  the  propagation  of  a slu'ar  pvilse  pai'allel  to  the 
interfaces  of  ,i  \ier i od i ca  1 1 y laminated  mediuin.  The  propa.vtat  i tni  of  wave  was 
causoil  h\  impulsive  slu'ar  lo.ids  perpendicular  to  the  l;uHinate.  Me  emjiloyed 
the  etl'ective  stit't'ness  tlieory  to  movlel  the  composite  aiul  applied  t lie  method 
I'f  t lie  luMvl  ot  the-pulse  to  viht.nn  an  approximate  ■^oluiion.  It  turns  out 
that  the  solution,  s.ii'  tor  particle  velocity  nonii.i  1 t (>  the  l.iminates.  is  the 
head-of  the-pulse  I'unction  for  inijiulsive  stress  normal  to  the  lamin.ites  |r.'l 
reversed,  with  oscillations  of  increasinyt  tiiinilitiule  iirecedinj’  the  m.iin  pulse 
lolloweil  hy  a monotonic  rise  to  a iiond  imens  i ona  I amplitiule  tif  unit)'. 

.Svo  and  Whittier  iSdl)  a I si'*  invest  iyat  evl  the  one-dimensional  pulse 
propay.at  i on  in  an  oM  uiuo  laminated  half-space.  Usiny,  the  effective  stiff- 
ness theorv  and  the  head-of- 1 lie-pu I se  approximation  ayain.  they  studied  the 
impulsive  wavt's  in  .in  ohliiiue  lamin;»le  c.nised  hy  a st  ep-]'ressure  loadiiij;  on 
the  houndary.  It  is  interest iny  to  note  that  the  res|x>nse  of  an  ohlinue 
l.imin.ite  is  both  anisotropic  and  yeomet  r i ca  1 1 y dispersive.  I'lie  ani>iOtropy 
leads  to  a twii-step  res|Mnse  to  a step  loadiny.  while  the  >1  i spers  i v 1 1 y 
smooths  ot  f the  two  ste|'s  and  causes  osc  i I I ,i  t i ons  preceiliny  or  t'ollowiny 
the  sti'ps.  I'hey  found  th.it  tlii'  character  of  the  response  chanyes  .is 
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the  aiiyle  ot'  l.imin.itt'.  o.  is  varioil.  When  o 0.  a sinyle  piilst*  propayati' 


(.1 


throiij;h  tlic  compos i to;  ;is  ot  is  iiicfcasoil , tlio  juilso  splits  into  two;  wlu'n 
a = '.U)“  , the  ros;x)nso  is  oiico  aj;ain  a siii^lo  i)iilse. 

Svo  jSlSl  also  aj^pliiHl  t lio  oft'octivo  stifl'iu'ss  theory  to  tlio  response 
jH'ohlem  ot'  a per  iovl  ica  1 ly  laminated  elastic  halt'-S|)ace  snl>iect  to  rapid 
internal  heatinj;.  l'i'opai;at  i on  tlirection  inirallel  to  the  planes  ot'  the  lami- 
nates was  considered.  The  analysis  indicates  that  dispersion  .issumes  an 
important  role  in  thermally  induced  stress  waves.  The  pulse  in  the  far 
field  may  hear  little  resemhlance  to  the  pulse  that  orii;inated  at  the  front 
sur!  aci'. 

.Sve  and  Herrmann  (SI  | also  apiilied  the  effective  stil'fness  theiiry  to 
study  the  dyiumiic  resiwnse  of  a periodically  laminated  half-plane  suhjected 
to  a movii\j;  load.  I'he  laminations  are  parallel  t tlie  surface  of  the  half- 
plane  .md  the  travel  inr,  load  is  .>  step  loavl  of  magnitude  I'  inclined  ,it  an 
ancle  4*  t'-’  tl'c"  surt'ace.  I'he  velocit\'  (.•!'  the  load  is  constant  and  supersonic. 

r!ie\'  obtained  the  formal  solution  with  the  aid  of  l.aptace  transfonns  and  con- 
sti’ucted  a far-field  solution  with  the  head-of- 1 he-pu  1 sc  procedure.  They 
catculate>.l  the  normal  str.iin  at  the  rear  face  for:  1.1)  movinc  normal 

lo.id  ; 4'  n .'1,  nuninc,  shi-.ir  lo.id  (4'  (M  . Ihe  la-su  1 ! -.  indicatf  linn 

the  primary  eft'ecl  v'f  the  l.im  i n.it  i ons  is  to  cre.ite  .t  '^tr.iin  that  has  a linite 
rise  time  and  is  oscillatory  about  the  value  obtained  on  the  basis  of  the 
effective  modulus  tlieors'  which  prcilicts  a two-step  resivnso  tor  4*  d and 
a rectancular  res’ponse  for  4'  = n/2.  Also,  for  nonnat  loadinc  tensile  str.iin 
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!\ortl\  .liter  the  .irrii.il  of  the  slusir  w.ive. 


WlutiK'v  aud  Sun  IKMI  presont  od  a rofini'd  o t'Ccct  i vi-  st  i flnoss  tlu'i>r\’ 
t'oi'  oxtonsional  mot  ion  of  laminalod  compositos.  In  this  tiu'ory,  toi-iiis  of 
sooond  ordof  in  _ tvliioh  is  t lio  ooordinato  ii\  tho  di  root  ion  of  t lie  thickness 
of  1 ami  nat  e are  i nc  I uded  I'or  the  loni;i  t lul  i na  I displacements  wiiile  only  the 
first  order  terms  in  z are  included  for  the  flexural  displacements.  I'he 
analysis  is  appropriate  I'or  oblique  impact  loadiiiit-  In  comparison  with  other 
theories,  tlis[)ersion  curves  appear  to  verify  this  refined  theory.  Witli  some 
mod  i I'i  cat  i on  , the  theor>’  is  aiqilicahle  to  fi  her- re  i n forced  composite  materials. 

Uedford  and  I'rumheller  1IU(>|  further  discussed  the  higher-order  effec- 
tive St  i ft' ness  theory.  .\  set  of  displacement  and  stress  houiulan  conditions 
at  an  external  surt'ace  compatihle  with  the  h i v'.her-o rder  theory  were  presented. 
The  method  developed  can  he  used  to  obtain  boundary  conditions  for  theories 
of  arbitrary  order.  Also,  it  makes  explicit  the  relationship  hetween  the 
microdispl.icement  and  microstress  distributions  and  the  actual  d i spl.icement 
and  stress  distributions  on  the  surface  of  the  composite. 

bedford  and  nrumhel  ler  1IU|  applied  a second-orvler  effective  stit't'ness 
theori'  for  the  iliiiamic  behavior  i>f  elastic  l.iminatt's  to  the  problem  of  a 
laminated  hall-space,  with  interfaces  normal  to  the  boiuiJarv,  subjected  to 
harmonically  time  varyin>;  displacement  and  stress  distribution  at  the 
boundary.  The  finite  numher  of  modes  of  the  microstructure  theory  were 
found  to  be  sufficient  to  model  a unil'orm  normal  d isiUacement  boundan  cvindi- 
tion  but  not  a uniform  normal  stress  bouiklari'  condition.  I'he  solutions 
yieKl  the  constituent  displacement  and  stress  distributions  both  near  the 
boundari’  and  in  the  far  fii-ld  and  peniiit  an  assessment  ol  the  usefulness  of 
the  microstructure  theory  tor  such  boundary  v.ilue  problems. 

I'he  first  nonlinear  theory  of  the  efl'ective  stillness  type  was  obtained 
by  Uausch  lR.i|  who  analyzed  the  tiMiisient  w.ive  propagation  parallel  to  the 
interl'.ices  of  ;i  l.iminate.  In  deriviiii;  the  effective  stiffness  is|uations, 
he  expressed  H.imiltori's  principle  in  terms  of  .i  I'otenti.il  enerfjy,  which 
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rosiilts  in  nonlinear  eonst  i t iit  i ve  |>ro[)ert  ies , aiul  a kinetic  enerj;y  containing, 
lai\k;e  strains  and  disiilacements  to  I'roditce  a k iiieiiiat  i ca  1 l,v  nonlinear  rhei>ry. 
lie  then  soujtht  the  solution  to  the  effective  stiffness  equations  hy  a 
coordinate-i>ertur|)at  it'll  technique  hasetl  on  asymptotic  expansitnis  [112]  in 
stress  amplitude.  fhe  eiiuat  ions  for  the  first-order  perturbations  .ire 
linear  differential  et|uat  ions  which  c.in  be  solved  by  integral  transform 
techniques.  lie  approx  imat  eil  the  transformed  solutions  t'l'r  lar,i;e  times  b.v 
use  of  the  saddle-point  techniqiu-'-  I'hus,  he  obtaineti  a nonlinear  head-of- 
rho-pulse  at'prt'x  imat  ion.  I'he  first-ortler  theory  presented  sliows  lu'w  lu'ii- 
1 inear  effects  cause  the  i^eomet r i cal  1 y dispersetl  wave  to  term  a sliock  but 
tloes  not.  stive  the  second-order  effect  in  wliich  dispersion  limits  tlie 
steepness  of  the  wave.  tAinsequent  1 y , steady  w.ives  of  the  t ipe  found  experi- 
mentally cannot  form.  I'hus  a further  aj'prox  imat  ion  will  require  the 
development  of  a secoiKl-otxler  theorv.  However,  the  results  of  the  first- 
order  theory  provide  an  insij;ht  into  tlie  n.iture  I'f  the  interaction  between 
iteometric  disj'ersion  .iiul  nonlinear  constitutive  effects.  It  appears  both 
nonl  ine.ir  i t y and  ,ceometric  d i spers  i v i t >■  liave  the  effects  to  c.iuse  early 
at  t eniiat  ion . 

lll.S  Interact in^-Cont inuum  Theories  and  Mixture  Theories 

Based  on  the  iiitei’act  in,i;  continuum  concept  in  winch  every  constituent 
has  its  own  motion  but  is  allowed  to  interact  witli  others.  Bed  lord  |B12| 
developed  a basic  nonlinear  tlieory  for  composites  and  explored  tlie  implic.i 
lions  of  this  tlieory  witli  reitard  ti'  shock  wave  propa,cat  ion . He  found  that 
t i«  li  constituent  can  supi'ort  separate  shock  w.ives  so  th.it  multiple  slu’.  k 
waves  can  occur  in  such  a theory.  Such  multiple  shock  waves  mav  be  more  .in 
artif.ict  of  the  theorv  th.in  an  essential  part  of  the  miwlelin^,  and  prob.iblv 
are  simil.ir  to  the  ch. tract  eri  st  ic  w.ive  speed-,  in  the  .ipprox  im.it  e theories 
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of  rods  and  heams,  wl\icl\  aro  not  nocossaril)'  olosol)’  assoriatod  with  tlie 
Oil  I cu  1 a t ion  of  the  niiiin  response.  In  this  theory,  the  iippropriiite  lunindar\- 
conditions  iiud  Uankine-llui^oniot  equiitions  for  ;i  imi  1 1 i -cont  immm  model  of 
composite  m;iteriiil  were  foniui  1 ;it ed . 

Bedford  liur>l  extended  the  mu  1 1 i -cont  i niinm  theory  to  fonnulate  tlie 
jump  conditions  for  mass,  momentum  ;ind  enerjiy  iicross  ;i  sinj;ular  w;ive  sui'- 
face  propiijiat  in^  in  ii  mu  1 1 i -cons  t i tuent  honded-oliist  i c composite,  lie  found 
thiit  the  shock  surfaces,  linked  hy  the  reference  conf  i_i;urat  i on , e.vist  for 
each  constituent.  Tlie  results  ;ire  relev. int  to  the  study  of  iicce  lei'iit  i on 
wiives  iind  sliock  w;ives  and  m;iy  be  adojited  to  ijive  bouiuhiry  cimditions.  The 
method  has  been  ;ip\ilied  successfully  to  sinjtle  cont  iniui  and  c;in  be  applied 
to  ii  broiid  cliiss  of  continuum  theories. 

Iler.emier  and  Niivfeh  lll.'|  devi'lopcil  ii  continuum  theory  biised  on  an 
,1  siiiipf lit  i c expiinsion  for  wiive  prop, u’at  ion  in  l.imin.ited  comiiosites  in  which 
dominiint  sii^nii!  wave  leiiitt  hs  were  .issuiiied  to  be  l.irj;e  comp:ired  to  tyi'iciil 
com|iosite  micro-dimensions.  I'he  ciise  of  \irop;ii;at  ion  norm.i  1 to  the  l.imiiKites 
w.is  considered.  I'liey  obtiiined  .i  hieriirchi’  of  models  defined  by  the  order 
ot'  trunciition  of  the  .isymptotic  senuence.  In  princi\ile,  retention  of  ii  1 1 
terms  in  tlie  iisymptot  ic  sequence  yields  the  exact  theory.  lliey  then  cast 
the  lower-order  dispersive  theori'  in  .i  stand, ird  mixture  f'orm,  Thei'  obt,iineil 
the  transient  pulse  res[ionses  at  the  centroids  of  the  liith  layer  and  IJth 
layer  respectively  due  to  a step-stress  ininit  on  the  boundary.  Comparison 
of  the  result  with  exact  analisis  obt,iined  from  a one-dimensional  numerical 
code  shows  quite  nood  ajireement  , It  is  noted  that  the  one-third  points 
(l/.i  first  pe,ikl  on  the  tiMiisient  waves  tr,ivel  ,it  the  mixture  velocit.v,  A 
similtir  result  was  found  in  tlie  head-ot'-t he-pu  1 se  analysis  of  Peck  and 
Curt  man  (P,'|. 
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llo>;i'ini  or , (.lurtman  aiul  Nayt'oli  lllSj  oxtondotl  t l\o  iiitoraot  iti^  ooiitituium 
thoory  h>'  a rational  i>oiist  root  i ot>  toohiiiqiio  for  ti'i'  const  itutivo  oqnations 
anJ  intoraotinjt  form.  Thoy  cons  itlorod  tlio  caso  of  wavoguido-Typo  profiaya- 
tion  in  laminated  and  f i hor- ro i nforcod  composito.  Utilizinj;  mixed  s]Kitial 
and  asMiiptotic  o.\i>an.s i ons , they  found  that  their  asym\'totic  approach  leails 
to  a first-order  expression  for  the  momentum  interaction  term  identical 
in  form  with  Boilfonl ' s Ihll),  hut  the  constitutive  equations  for  each 
constituent  involve  the  strains  in  both  constituents.  Hence,  the 
result  injt  first-ordei-  continuum  mixture  theory  contains  microstructure, 
rhe  utility  of  this  theoi-y  demonstrated  for  both  laminated  and  fibrous  com- 
posites by  correl.it  iny  theoretical  and  experimental  jlVKS]  transient  jiulse 
(.lata  on  boron-carbon  jihenolic  and  I'horne  1 -carbon  phenolic  laminated,  and 
unidirectional  fibrous  quart ;-pheno 1 i c . In  each  case,  the  input  pulse  is 
a step  function  in  velocity  applied  to  both  constituents  at  the  boundary. 
I'he  response  at  the  rear  surface  is  a steady  rise  of  amplitude  followed  b>' 
oscillations.  In  all  cases,  the  a.i;reements  between  theoretical  and  experi- 
mental results  are  quite  yood. 

Ileyemicr  and  baclie  |lh>l  ayain  discussed  the  transient  pulse  propavta- 
tion,  ]iarallel  to  laminates  in  e last  ic- laminated  composites.  based  upon 
asymptotii.  expansions  they  developed  a modified  first-order  mixture 
theor>-.  lor  a step  velocity  input  on  the  boundary,  the  transient  dat.i 
obtained  t'rom  the  simplified  theory  exhibits  yood  correlation  with  experi- 
mental results.  because  of  its  simi'licity  and  s.it  i s fact  or\  prcvliction, 
this  simplified  theor>'.  tlie\'  hoped,  will  leail  to  m.iximum  util  it'. 

bache  .iiul  lleyemier  |bl|  further  discussed  the  transient  wave  propaga- 
tion in  V i scoel. 1st ic- laminated  composites.  With  minor  mod i f i ca t i ons  of  the 
techniques  used  lor  elast  ic  cases,  the\'  constructed  a yeneiMl  theory  tor 
wave  [)ro|Mgafion  in  v i scoe  l,i  s t i c - 1 ,im  i na  t eil  comi'i'sifes. 


I'be  [ifocedurv's  are 
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quite  straightforward.  They  first  applied  the  Laplace  transformation  to  the 
equations  of  motion  and  viscoelastic-constitutive  relations,  then  they 
deduced  the  continuum  theory  for  the  problem  under  consideration  in  tlie 
l.a{)lace  transform  plane  from  the  well-known  correspondence  principle. 
Finally,  they  took  Laplace  inversions  and  obtained  the  theory  in  the 
pliysical  space.  In  view  of  its  simplicity  and  potential  utility,  they 
derived  the  corresponding  viscoelastic  modified  first-order  tiieory,  and 
applied  this  simplified  theory  to  investigate  the  transient  pulse  nonnal 
to  the  laminates  due  to  a step-velocity  input  on  the  boundary  of  a half- 
space formed  by  the  bi laminates.  Hic  composite  was  laminated  by  a quartz 
fiber  bundle  in  a phenolic  matrix  ^constituent  1)  and  Ironsides  FF-17 
phenolic  resin  Iconstituent  2).  The  quartz  fiber  material  was  modeled  as 
elastic  for  propagation  normal  to  the  fibers  while  the  standard  linear 
solid  model  of  viscoelasticity  was  used  for  the  phenolic  resin  material. 

They  compared  tlie  stress-time  history  at  the  midplane  of  the  13th  layer  to 
that  for  a corresponding  linear  elastic  composite.  It  is  seen  that  the 
general  shape  of  the  former  is  more  or  loss  similar  to  that  of  the  latter 
but  with  the  ;implitudc  of  the  peak  stress  and  the  oscillations  following 
the  peak  d;unpcd  out  somewhat. 

Nayfeh,  I’lurtman  and  llegcmier  [N31  used  the  theory  of  interacting 
continuum  to  study  clastic  wave  propagation  normal  to  the  fibers  of 
unidirectional ly-reinforced  composites.  The  composite  consists  of  a 
periodic  array  of  square  fibers  which  is  perfectly  bonded  to  the 
matrix  material.  Tliey  derived  a "first-order"  ty])e  dispersive  model 
{ for  this  complex  geometry.  The  basic  construction  procedure  is  as 


follows:  (1)  they  modeled  the  fiber-matrix  regions  as  a mixture  using 
the  same  procedure  for  layered  materials;  (2)  they  modeled  the  matrix 


material  rojiions  as  homogeneous  laminates  under  unidirectional  motion  in  tlu* 
direction  normal  to  tlie  layers;  1,3)  they  then  followed  the  same  procedure  as 
that  for  obtaining  wave  propagation  normal  to  layered  composites  to  conctruct 
a ilispersive  continuum  tlieory  for  the  fiber  reinforced  comiuisites.  The  simjili 
fieil  first  order  continuum  model  so  constructed  is  seen  to  he  similar  in  form 
to  the  theory  obtained  by  them  for  [)ropagation  normal  to  the  layers  of  a lami- 
nated composite.  The  only  difference,  in  fact,  is  the  complexity  of  the  con- 
stants in  the  equations.  Utilizing  this  theory,  they  obtained  the  solution 
for  a step-stress  input  at  the  boundary  of  a semi-infinite  medium  of  a uni- 
directional fi  lu'r- re  in  forced  quartz  phenolic  composite.  I'he  step-function 
pressure  was  applied  to  tlie  front  surface  of  the  material  which  was  a fiher- 
matri.x  laminae.  I'he  calculated  rear  surface  velocity,  the  velocity  corre- 
sponds to  that  of  a matrix  laminae,  vs.  time  curve  is  again  a steady  rise  of 
magnitude  followed  by  oscillations.  Comparison  of  experimental  data  reported 
by  the  .Aerospace  Corporation  lllll)  and  theoretical  results  showed  that  agree- 
ment is  fair.  If  the  geometrical  approximation  is  improved  and  the  visco- 
elastic nature  of  the  phenolic  is  accounted  for,  tlie  agreement  should  he 
much  better. 

Ciurtman,  Nayfeh  and  llegemier  [Clo]  generalized  the  theory  of  interacting 
contiiunim  ITINCI  to  study  tl»e  t wo-d imens i oi\a  1 wave  propagation  in  structural 
composite  materials.  In  their  generalization,  a caloric  equation  of  state 
was  assumed  for  mean  stress  vs.  density,  and  an  elastic,  perfect ly-plast ic 
law  of  von  Mises  wa-  assumed  for  the  deviation  stress  tensor.  In  addition, 
two-dimensional  lay-ups  (laminates)  were  considered.  I'hey  then  developed  a 
hierarchy  of  two-dimensional  motlels  of  interacting  continuum  defined  b>  the 
order  of  trunc.it  ion  ot'  the  .isymptotic  ex|iansion.  Again,  the>  focused  tlieir 
.ittention  on  tlie  simplest  "first-order"  tlu'i>ry  hoping  that  its  relative 
simplicity  will  lead  to  maximum  lit  i 1 i t v while  preserving  the  desired 
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micromechanical  information  rei|ulred  by  the  composite  iJesii;ners  and  the 
analysts.  Transient  resixnise  to  a st ep- fiinct  ion  in  velocity  at  the  boundary 
of  either  the  elastic  medium  or  the  elastic,  perfect  ly-plast ic  medium 
correlates  well  with  the  correspond inj;  experimental  results. 

Nayfeh  and  tiurtman  [N4|  extended  the  interacting  continuum  theory  to 
study  the  transient  shear  wave  motions  in  huiiinated  waveguides.  Both  trans- 
versely (.SV)  and  horizontally  (^Slll  polarized  waves  were  considered.  Transient 
solutions  for  both  types  of  waves  were  investigated  via  the  lu'ad-of-the-pulse 
approximations.  Tlie  transient  response  of  tlie  SH  waves  to  a step-boundary 
velocity  w^^  is  a pulse  beginning  with  a steady  rise,  followed  by  oscilla- 
tions abtxit  tlie  input  boundary  velocity,  whi  le  tlie  res(Kinse  of  the  SV  waves 
to  a step-boundary  velocity  v^^  is  a pulse  starting  with  oscillations  about 
zero  followed  by  a smooth  rise  to  its  final  value  (.the  boundary  input).  It 
is  noted  that;  (.1)  the  transient  solution  of  SV  wave  is  identical  to  that 
obtained  by  Sve  [S40|  based  upon  the  effective  stiffness  theory,  and  1.2)  the 
transient  solution  of  Sll  wave  closely  paral lels that  of  the  longitudinal 
wave  solution  of  beck  and  Gurtman  [1’3|  obtained  by  the  exact  theory  and  that 
of  llegemier,  et  al.  lllS)  obtained  by  the  continuum  mixture  theory. 

Gurtman,  Nayfeh  and  llegemier,  et  al.  |t'il4]  applied  the  tlicor)’  of  inter- 
acting continuum  to  stud)’  wave  propagation  in  three-dimensional ly  reinforced 
composites.  The  model  developed  is  an  extension  of  previously  developed  two- 
dimensional  l;ui\inate  and  fiber  I'lNC  formulations  and  explicitly  considers 
the  effects  of  tliemodynamics , finite  defonnations  and  nonlinear  constitutive 
behavior.  The  three-dimensional  material  is  made  up  of  "radial"  and  "lateral" 
fiber  bundles  in  an  orthogonal  array,  with  pockets  between  them  filled  with 
pure  resin.  They  treated  the  fiber  bundles  as  isotropic,  but  using  their 
material  profierties  in  tin'  direction  of  jiropagat  ion.  Thus,  material  properties 
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along  the  fiber  are  assigned  to  tlie  radials,  and  those  transverse  to  the 
fibers  are  used  for  the  laterals.  Basically,  tliey  proposed  modeling  the 
unit  cell  of  three-dimensional  quartz-pheiiol ic  as  a combination  of  cylin- 
drical and  kimii\ate  waveguides.  The  initially  square  cross  section  of  the 
radial  fiber  is  transformed  into  a cylinder  of  equivalent  volume  fraction 
surrounded  by  an  annular  sheath  whose  material  properties  are  dyn;unical  and 
dependent  upon  those  of  tlie  alternating  fiber  and  resin  pocket,  for  the 
cylindrical  waveguide,  the  individual  constituents  are  modeled  as  elastic- 
perfect  ly-plast  ic  material,  witii  a Mie-(u'iineisen  caloric  equation  of  state 
relating  mean  pressure,  density  and  internal  energy,  and  with  a von  Mises' 
yield  criteria  and  associated  flow  nilc  governing  the  stress  deviators. 

Thus,  they  developed  a themodynamic  theory  of  iivteracting  continuum,  of 
waveguide  type,  for  finite  ;miplitude  elastic-plastic  wave  jiropagation  in 
fiber-reinforced  composites.  Tlicy  then  smootlted  the  surrounding  structural 
elements  of  the  actual  material  made  up  of  lateral  fibers  and  resin  pockets 
via  the  TINC  laminate  analysis  into  a hollow  cylinder  slieath.  The  complete 
titermoelast  ic-plast  ic , three-<limensional  TINC  model  then  is  obtained  In-  com- 
bining the  solid  waveguide  fiber  with  the  dispersive  hollow  cyclinder  sheath, 
for  elastic  pulse  propagation  iti  three-dimensional  quart  z-plieno  1 ic , tiiey  used 
the  finite  difference  TlNtl  code  to  compute  the  average  velocity  at  the  rear 
of  a specimen  0.035  cm  thick.  A step- boundary  velocity  is  applied  to  both 
filler  and  sheath.  They  found  that  the  ve  loci  t y- 1 ime  res]H)nse  is  again  a 
stead)'  rise  followed  by  oscillations  about  the  input.  llie  TINC  ci>de  result 
was  compared  with  experimental  data  obtained  at  Ai'Wl,  (Air  forces  Weajions 
Laboratory),  and  with  the  head-of-the-pulse  aiiproximat i on.  The  agreement 
is  excellent.  The  head-of-t!ie-pu  1 se  solution,  wl\ile  effectively  duplicating 
the  peak  velocity  and  rise  time  of  the  w.ive,  tends  to  diverge  from  the 


riM’ results  followini;  the  passaj;e  of  the  front.  Ti\ey  also  used  the  TINC  code 
to  calculate  the  stress  lustories  at  the  rears  of  three  three-dimensional 
quartz-phenolic  specimens  of  different  thickness  which  were  tested  at  the 
Al'WL  usins  a light  gas  gun.  Comparison  of  code  predictions  and  experimental 
results  indicated  that  there  appears  to  he  a tendency  on  tlie  part  of  the  code 
to  slightly  overpredict  the  rise  times,  but  agreement  of  peak  pressure  is 
excellent. 

For  nonlinear  pulse  propagation  in  tliree-dimensional  quart  z-pl\eno  1 ic 
they  used  the  corresponding  I’lNC  code  to  calculate  the  stress  t\istories  at 
the  rears  of  several  oHQl'  specimens  which  were  impact -t ested  at  Al'WL  using 
a high  pressure  gas  gun.  Tlieir  nuiHcrical  results  and  experimental  data 
indicated  again  that  the  TINC  code  tends  to  slightly  overpredict  rise  times, 
but  accurately  calculate  the  peak  pressures,  and  significant!)’  diverge  from 
the  experimentally  detemined  curves  after  the  peak. 

The  theory  so  developed  has  been  shown  to  be  capable  of  predicting  the 
nonlinear  response  of  50QT  to  both  mechanical  ai\d  thermodynamic  loading,  and 
capable  of  yielding  infomation  concerning  the  behavior  of  the  three- 
dimensional  material's  const itucnt s.  It  should  be  a useful  tool  in  both 
design  and  iuialysis  of  three-dimensional  composites. 

Uegemier  [1112]  discussed  the  finite-;unpl  itude  elastic-plastic  wave 
propagation  in  huninatcd  composites.  He  developed  a binary  mixture  theory 
for  waveguide  tyi’>e-propagat  ion  parallel  to  the  layers  of  a two-constituent 
laminated  composite  witii  periodic  microstructure.  Tlie  model  incorporates 
the  effects  of  thermodcTiiunics , finite  defoniiat ions , aitd  nonlineai-  elastic- 
plastic  constituents.  For  a step-boundary  velocity  input,  transient 
displacement,  stress  and  internal  energy  d i st  ri  l>ut  ions  witljin  the  micro- 
components are  produced  to  a certain  degree  of  accuracy. 
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Hegemier  and  Ciurtman  [HlaJ  also  imostigatod  tlic  t'inite-ajiiplitude 
elastic-plastic  wave  propagation  in  fiber-reinforced  composites.  Accordingly, 
they  developed  an  approximate  nonlinear  theory  to  describe  waveguide -ty^ie 
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propagation  in  unidirectional  fibrous  composites.  Tlie  model,  an  extension 
of  the  previously  developed  l;uninate  foimiulat  ion , considers  the  effects  of 
thermodyaiiimics,  finite  deformations,  and  nonlinear  elastic-plastic  constitu- 
tive behavior.  A one-dimensional  binary  mixture  tlieory  is  thus  followed. 

Transient  wave-propagation  solutions  due  to  a step-boundary  velocity  input 
were  obtained  numerically  and  compared  with  another  solution  from  a well- 
known  two-dimensional  finite  difference  code.  nie  agreement  is  judged  to  be 
excellent.  However,  the  authors  concluded  that  the  TINC  calculation  is  much 
more  time-saving  than  that  of  the  finite  difference  code. 

Aboudi  [AoJ  extended  his  previous  work  [Ad]  and  fomiulated  a mixture 
theory  for  a thennoelastic  laminated  medium  composed  of  two  constituents  in 
alternating  layers.  In  this  theory,  each  constituent  lias  its  own  motion  and 
temperature,  but  is  allowed  to  interact  mechanical 1\-  and  thennally  with  the 
others.  The  resulting  system  of  coupled  equations  of  motion  and  heat  conduc- 
tion was  then  used  to  study  the  res{x>nse  of  a l;uiiinated  plate  subjected  to 
meclianical  and  tliermal  impulsive  loadings.  Comparison  between  tlie  results 
of  this  theory  and  those  based  on  tl.c  tltermoe  last  ic  effective  modulus  theory 
which  he  formulated  exhibits  the  pronounced  effect  of  microstructure  and  the 
effect  of  the  reinforcement  volume  on  the  resulting  field  in  the  individual 
constituents-  The  effect  of  the  reinforcement  is  to  decrease  the  stress  in 
the  matrix  and  to  increase  it  in  tlie  reinforcement. 

Ai>oudi  and  Benveniste  [A7|  proposed  a superimposed  mixture  theory  for 
wave  propagation  in  a liiaxially  fiber- reinforced  composite.  I'he  mevlium 
I studied  has  alternating  I avers  in  which  the  angle  of  the  biaxial  fibers 
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alternates  from  layer  to  layer,  liach  constituent  has  its  own  motion  but 
interacts  with  others,  llie  layering  microstructure  is  taken  into  account 
by  the  interface  sliear  stresses.  The  composite  is  then  in  tlte  fonn  of  a 
l;uiunated  medium  composed  of  uniiixial  fiber-reinforced  material  in  alternating 
layers.  They  modeled  each  layer  by  Bedford's  mixture  theory  and  modeled  the 
composite  by  Hegemier's  binary  mixture  theory.  In  other  words,  the  model 
is  a superposition  of  the  two  knowai  theories.  Hiey  calculated  the  response 
of  a material  consisting  of  carbon  phenolic  reinforced  by  boron  fibers  for 
pulse  propagation  in  the  direction  parallel  to  the  layers.  In  contrast  with 
the  uniaxial  ones,  numerical  results  showed  that  the  stresses  in  the  fibers 
in  the  case  of  biaxial  reinforcement  attain  higher  values  as  compared  with 
the  uniaxial  one.  The  normal  stresses  in  the  matrix  are  almost  identical  for 
both  types  of  reinforcements.  It  would  be  of  interest  to  compare  the  pre- 
dicted results  with  experimental  work  and  thus  evaluate  the  validity  of  tlie 
proposed  model. 

Nayfeh  [NS]  treated  the  transient  pulse  propagation  in  porous  composites. 
He  employed  the  interacting  continuum  theory  to  analyze  the  influence  of 
inclusions  and  porosity  on  tlic  clastic  response  of  both  homogeneous  and  l;ijjii- 
nated  composite  media.  Hie  general  model  analyzed  iiy  him  consists  of 
periodic  array  of  two  perfectly-bonded  kiminates;  one  of  which  consists  of 
an  elastic  homogeneous  material  while  the  otlier  is  made  up  of  periodic  array 
of  cylindrical  elastic  inclusions  that  are  distributed  in  another  elastic 
matrix  material.  He  dcduccil  several  s]iociric  models  as  special  cases.  In 

all  cases  porosity  is  simulated  in  the  limit  as  the  properties  of  the  inclu- 
sions identically  vanish.  He  demonstrated  that  porosity  plays  a major  role 
in  the  geometric  dispersion  of  sucit  media.  In  particul.ir  it  increases  the 
pulse  arrival  and  rise  times  ispreading)  of  a transient  pilse.  I'or  the 
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Special  case  of  elastic  inclusions  in  a homogeneous  matrix  media,  the  results 
correlate  very  well  with  existing  experimental  data  and  other  approximate 
analyses . 

Nuyfeh  [N7]  extended  the  interacting  continuum  (.or  continuum  mixture) 
theory  to  discuss  the  heat  conduction  in  lajiiinated  waveguides.  Tlie  resulting 
theory  leads  to  simple  governing  equations  for  the  actual  composite  which 
retain  the  integrity  of  the  diffusion  process  in  each  constituent  but  allow 
them  to  coexist  under  some  interactions,  lie  then  utilized  the  resulting 
eouations  to  study  the  response  due  to  transient  loadings.  Solutions  were 
derived  by  means  of  Laplace  transform  techniques.  Analytical  inversion  of 
transforms  was  carrioti  out  only  for  tlie  limiting  cases  of  "weak"  and  "strong" 
thermal  coupling.  The  limit  of  strong  interaction  leads  to  the  coalescence 
of  both  temperatures;  in  this  limiting  case  tlie  composite  behaves  as  a single 
but  high-order  continuum  according  to  his  investigation.  For  the  general 
coupling  cases,  his  results  were  obtained  by  a direct  numerical  inversion  of 
the  transforms.  Since  his  results  for  harmonic  loadings  correlate  well  with 
some  existing  exact  solutions,  and  those  for  transient  loadings  are  reason- 
able from  physical  considerations,  the  theory  sliould  be  useful  in  studying 
any  transient  temperature  pulses  in  laminated  composites. 

Gurtman  and  llegemier  [G15J  modified  a previously  developed  binary 
mixture  theory  to  account  for  debonding  for  waveguide  propagation  in  laminated 
composites  with  periodic  microstructure.  In  modification,  they  relaxed  the 
interface  shear  stress  boundary  condition  and  proposed  an  interface  model  of 
the  Coulomb  frictional-type.  They  then  applied  the  resulting  theory  to 
analyze  the  response  of  a wavegu  ide-type  l;rninate  to  imjuict  loading  at  the 
boundary,  and  compared  the  results  with  exi)er iment a 1 results.  A fairly  good 


agreement  was  observed. 


/ / 

I'ouviMii  St  0 auil  Ahouil  1 also  ox  ( oikIoJ  t ho  mixtiiro  tiu'orv  to  acoount 

t'or  ilohoiuiinj;  iluo  to  iiiii>ulsivo  loailiii>;  i'or  wave  pfopajiat  i on  in  a laiiiinatoil 
moJiiim.  IVav'o  [tropaxat  ion  is  in  I ho  ilirootion  of  i ho  la\oiin};  o('  a hi 
liuninatoil  moilium  with  tlio  \irosonco  ol'  impoorool  hoiulinp,  at  t ho  intort'aoos. 
llto  ilohoiulin^  is  moiloloil  h)'  a I'loxihlo  boiui,  i.o.,  an  inortia-loss  thin 
olastio  ('ilm,  which  was  or  i i na  1 1 y proposoii  In'  ,lonos  aiul  IVhittior  [.U}. 
rhoroforo,  the  (.lohomi  i ng  moohanism  is  roprosont  oil  by  a iiuulol  whioh  allows 
iinporfoct  hoiulinji  both  in  tho  noniial  aiul  tanponi  ilirocl  ions.  broiii  this, 
t hoy  to  rum  1 at  Oil  .i  nioilil'ioil  mixtiiro  t hoory  ami  apjtlioil  tho  theory  to  traimiont 
wave  propa.yation  in  tho  wavoy.ii  iilo-ty|H'  laiiiinatoit  oomposito.  It  was  roiiiul 
that  tho  ilcbonil  1 111’,  in  t lio  tanjtoiit  ilirootion  is  s i j’.ii  i I' i can  f in  iiioil  i t'y  i iij;  tho 
sha\H'  ami  ain)i  1 i tiiilo  ol'  tho  propaija  t i ni;  piilso. 

As  poiritoil  out  prov  loii-.l  \ , I'on  Anuv,:  ilovolopoil  a oont  imium  thoory  ol' 
ooiiii'osilo  malorials  in  a sorios  ot'  papor>  | It  J J . J.i , JS  . J'.)  | iisiiij;  iliroot 
asMiijitot  io  expansions.  I'lio  iloi  o lopmoni  s ol'  tho  Ihoori'  imlioalo  that  tho 
liynaiiiio  behavior  ol'  the  moih'l  is  "non  I oo.i  I " in  I i mo  a-',  a result  of  Iho  hislorv- 
ilopomlonoo  of  tho  oiirront  state.  from  this  Ihoorv  ho  iloiluooil  a ;oro-orilor 
moilol  for  praotioal  applications,  lie  analy.-.oil  Iho  behavior  ol'  this  zero 
oi'ilor  system  ami  foinui  (hal  iliiriiii;  an  early  phase  tho  mol  ion  is  coni  inoil  to 
.1  Iviiiiiilary  layer  ami  consists  of  hiy.hly  ilampoil  waves.  Uiriny  a later  ph.iso 
Iho  behavior  approaches  that  of  a mac ro'Coii i ca  1 1 y lu'iiioy.onooiis  moiiiiim.  The 
boh.ivior  ihirinj;  both  ph.isos  is  iloscriboil  In  two  ilistinci  systems  of  ilifl'oron 
till  oi|ii.if  I <tis . Ho  then  applioil  this  ihoorv  to  Hot  ermine  Iho  early  phase 
!.«»»■  II  •!  Oil  iniin.itoil  hall  sp.ico  siibioct  to  .1  stop  noriii.il  loml  for 

• ••  i - ii:  I to  the  I.iminitos  |l(.’,'t|,  .iml  (J)  prop.nj.it  i on  norm. 1 1 

. t 'I  III’  in’  .i|'proc  1 .11  Oil  ilifforoncos  from 


I' ! ho  r I hoi'i'  1 os 
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lU’ii  Aiiio;.  I lU' ’ j tui'llu'f  » iMi'' uU' i\'vl  I'l.i'.tii'  pl.r.tiv  w.ivos  i\i  I anu  i\.it  Oil 
(.a’mpi's  1 1 1'S . lia'>i'i\  i'i\  t l\i'  a';'.iim)'t  i oi\  ot  por  t ('v  t 1 \' - 1' I ast  i v'  t ibor  lavi'vs  anit 

0 I as;  ii-  p I ast  u'  malriv  lavt'rs,  a la'ut  iwiium  minlol  wa*.  ilislucisl  foi’  waxo  mot  u>i\ 
oithor  paralU'l  to  oi'  aormal  to  t ho  vlii'ootioii  I’t'  tasi'rs.  I'lu'  vloi  ivation  was 
av’,aii\  hasoil  oi>  as\im'tot\o  o\\',\iis i ons  in  toiins  v>t'  a small  paraiiii'tor  so  tliat 

a oont  lumim  thooi'v  was  ol't.iiiusi  id  wliu'li  m i ».  ros  i nii  t u it'  t'llt'ots  appo.novl. 

A^ain  ho  l\'untl  that  iluiiiip.  tho  ('arl\  phasi'  t lu'  nu>t  u'w  \s  ov'ul  luovl  tt'  a hvHintiaix 
Lui'i'  aiwl  Oi'dsists  i>i  hii'JilN  ilampoJ  wa\t's  whoi'oas  tho  hohax  ior  vhirinv'.  tho 
lator  phasi'  is  pi'ovU'm  i uaut  1 \ tliat  v'l  a mao  rosoop  i o a 1 I \ lu'mt',v',oiu'v'us  movlium. 

I'ods  I iloi' i di;  that  piilso  at  t onuat  i t'li  v'vours  mainl\  ilurinp,  t lu'  oai'l\  phast',  ho 
oiMio  liulotl  that  tho  oar\\  I'h.iso  iiu't  ion  ov'iis  t 1 1 u t os  t lu'  oritioal  pliaso  tM‘ 

mot  ion  as  in  tho  ol.istio  wa\o  oa-o'.  I'ho  iiuwlol  shouKI  ho  adotpiati'  tiir  vom 

positos  that  aro  'itVi'iipU  i nhomop.oiU't'iis , 

1 I I . (<  Othor  Mt’tluwl'>  aiul  I'tohU'iii'-  Kolatovl  1 1'  i\'m|,'i's  1 1 1"- 

Sti'ni,  Itovlt'i'ial  aiul  \i'w  aiui  Itoilioial  aiui  Slorn  (l'li>|  analv.'Ovl 

wa\oy’,uulo  propa)',.it  uMi  in  la'ori'il  aiul  iin  ul  i it'o  I i t'lia  1 I'lhi'i'  rointoioovl  ov'iiipv's  1 1 1"- 

oonsistoii  ot'  t'Kistio  aiul  i i s^■,^o  I a s t i o matorial-,  rho\  a.sumt'vl  a simi'liiioa 

il  1 s[>  I aoomont  lu'l.l  oi'iisistnii;  sM  im'tion  v'lili  in  tho  propao.at  i on  vl  i loo  t i on  . 
riu'  vl  i sp  1 aoomi'iit  s aro  vaiiistant  oior  tho  t'la.tio  ro  i n t’v'ro  i nr,  lan'i',  hut  \ai\ 
vuor  tho  ort'ss  sootum  ot'  tho  matriv,  I'lu'  la'sul  t i nr,  t'lpiatioiis  o rii  i nr,  suoh 
motion  woi'o  solvoil  lor  -<  i iiiiso  i via  1 waxi"-  aiul  lu'Uh'vl  \or\  ,i;i'vut  vt'itn'a  r i sv<n  s 
with  o\aot  si'liit  u'lis  III  tho  oasi'  vi|  laxoroil  v v'liipv's  i t o . Ilu'  s,imo  li'rm  I't 
siiliituni  was  .ipplioil  to  tho  (thi'r  ri’ i n (v'lv  I'd  mali'ri.il  with  v'lili  a yjot'irn'ti'u 
corroot  u'li  tor  tho  sovunul  vU'ilor  torm  in  tin'  K'w  t rovpu'iiox  oxpansuin  of  tho 
0 III  oiirvi',  aiul  iioUioil  r,o\Hl  n'liipa  r i svvii  with  tho  o\av  t svilutiv'ii  in  vault  rast 
with  tho  ollootixv'  stiftnoss  tlu'v'iiv's  whoro  t lu'  ay^rooniv-iit  with  tho  ovavt 
thovirx  IS  host  tor  laiyp'  rat  iv>s  vvt  fil.imv'iit  tv'  malriv  stiitnoss  (sov',  tv'r 
ovamplo,  [SlSli. 


I'aytoii  ll’-’l  iiwost  i l !u'  ilxiKiiiii  i'  bi'iiil  st  n'SN  in  >i  v'ompoi^  1 1 o striK'tin'e 

subjootoil  to  a suvklon  I'rossnco  rise,  I'bo  w'oiniH'S  1 1 o striKluro  t.H>iisists  n't"  twi' 
elastic  somi  - i n t'i  n 1 1 1'  roils  boi\Jod  toi;othor  aloni^  tboir  ,i;oiu''‘at  ors  . B>'  int  lodiu'- 

inv;  an  iittoiMot  ion  toiaii  inti'  tbo  ono  dunons  lona  I wa\i.'  oipiat  i ons , lio  obtainod  t bo 
dynamio  bottd  sboa.rinit  st  foss  toe  tbi‘  of'  a peos 'Jni’o  stop  snddonly  appliod 

ovoe  tbo  ond  of  tlu'  roil . I'bo  bonil  --t  ross  is  oomposovi  of'  a statu'  and  a ilinaiiiio 
part  and  tbo  ('oak  bond  st  i-oss  tluit  ooouri'od  at  tbo  loadod  ond  dominatod  by 
tbo  static  I'ai't  . Kosults  obtainovl  ina>'  peovido  information  to  dosi);nors  in 
minimiiin);  tl\o  possil'ility  of  bond  failuro  from  d'liamic  loading. 

Aobonbacb,  llomann  and  Zioj;lor  IMI]  -.tudiod  tbo  tonsilo  failuro  ot'  intorfaco 
bonds  in  a comi'o-iito  bod\  subioctod  to  coiiipross  i vo  loads,  t'.irrytnj;  out  tbo  ono 
dimonsion.il  .inalvsis  of  u.uo  propa.if.it  ion  norm.il  to  tbo  tntorfacos  of  a laminato, 
tbo>  found  tb.it  a oomprossivo  lo.ul  can  lo.id  to  tonsilo  st  ross  wbon  .1  w.ivi’  is 
t ransiiii  1 1 od  across  .m  intorf.ico  to  .1  lower  impond.inco  m.it  o ri  .1 1 , .ind  tli.it  fr.icturo 
conditions  also  do\'ond  on  tbo  dur.it  ion  of  tbo  lo.id  .ind  tbo  tlucknoss  of  tbo 
layers.  Wbon  .1  lo.id  ot'  dur.it  ion  loss  t b.in  t lie  tr.insit  time  of  oitlior  l.iyor  is 
appliotl  to  tbo  lower  imped. nice  l.i'or,  I'racturo  c.iusod  b>  tbo  tonsilo  wavi'  ft' 
fleeted  from  tbo  second  intorf.ico  will  occur  .it  tbo  first  intorfaco.  Ibis  predi- 
cation was  von  fiovl  bv  ovporimonts. 

I'iiij;  and  l.oo  (IS)  .ipplioil  tbo  ci'i’iiiot  r i c.i  I .ici'iistics  .ippri'.icb  1 1'  dotormino 
tbo  response  of  ombovldod  circiil.ir  .ind  splioric.il  inclusions  to  a I'l.ino  d i 1 .it  .it  i on 
wave  in  tbo  surrounding  modium.  Ono  of  tbo  most  1 n t o rest  1 n>t  results  obt. lined  bv 
tbom  is  tbo  OMStonco  ot'  c.iiistics,  wbicli  is  .in  onvolo\'o  of  r.i>  s , witliin  f lu' 
embedded  modniiii.  .Since  t'.uh  r.iy  c.irrios  .1  finite  s 1 v;n.i  1 . while  tbore  are  in 
finite  rays  moot  in .it  .1  c.iustic,  tbo  rospi>nso  is  infinite  .it  the  caustic.  This 
siiktjjests  that  nonlino.ir  response,  \'ossiblv  including;  fr.icturo,  will  occur. 

rtu'  >;oomot  r ic.il  .icoiist  ics  .ippro.icli  simil.ir  to  that  iisovl  by  I'ln^  .ind  1 1'o  w.is 
applied  by  .Acbonb.icb , Hem. inn  aiui  Zioi;lor  l-XPl  to  iiuosti^.ito  tbo  so\'.irat  ion  .it 
tbo  intorfaco  of  .1  circul.ir  inclusion  .iiwl  tbo  sur rouiuli n>;  iiH'ilium  under  .111  inculont 


HO 

oonipri's.si  VO  wave.  I'lioy  ilisous'^oO  t ho  naturo  ot  f lio  s i ai' i t y assooiatoO  vvith 
tho  caustic.  I'ho  s iiii;ulari  ty  is  of  a squaro  root  OoponOouco  on  Oistanco  to  t ho 
caustic  as  it  is  a\iproaclioO  aiiO  ttion  of  a l(.>,i;ari  thnii  c iloi'ciulonco  on  distance 
liohind  tlio  wave  front  onco  tho  caust  it  is  passoil.  iv'hilo  nonlinear  I'ospt'iiso  will 
modify  tho  results,  they  slunvod  that  vor\'  larr.o  stresses  can  bo  transmitted  to 
tho  botindary  of  tito  inclusion. 

l iny  and  I'hou  [ TOl  further  omiiloyod  this  approach  to  invest  iy.ito  the  piopay;i- 
t ion  of  stress  grailient  through  an  inclusion.  fhey  fouiul  that  t lie  rot'lection  .iiul 
transmission  of  tlie  stress  graiiient  and  the  higher-order  tie  ri  vat  i vi's  of'  stress  do 
depend  on  the  geometries  of  tlie  incident  wave  front  ,ind  the  interface  liouiul.iri 
though  it  is  not  so  for  t lu'  stress.  I'hev  ,ilso  studie>.l  tin'  propagation  ol  the 
stress,  the  stress  gr.ulient  .ind  the  lugju'r  order  dor  i i.tt  i les  i>f  the  st  I’ess  tH'hiiid 
the  K.ive  front.  fhey  showed  that  the  stress  gi-.ulient  does  not  maint.iin  tin'  same 
sign  as  tlie  wave  front  is  propag.it  ed . flie  implication  o('  this  result  is  t h.it  the 
plastic  yield  can  occur  .it  pi. ices  beliiiul  the  w.ive  front  before  it  occurs  at  the 
wave  front  even  if  the  initi.il  wave  has  .1  d 1 scon t 1 iiiu'us  rise  in  stress  .it  the  w.ive 
front  followed  bv  .1  gr.idu.il  decrease  in  iii.ign  1 1 iwle  behind  the  w.iie  I'rinit  . 

.\n  analysis  b.iseil  on  the  theory  ol'  el.isticitx  .iiid  aiiplic.ible  to  composites 
was  carried  out  b\  .■\cheiib.ich  (Air’l  lor  weUh'd  I’l.istu-  qii.irti'r  sp.ue.  fhe  c.isi-  I'f 
impulsive  shear  li'.ul  witli  st  ep- fuiict  1 on  time  dependence  .ipplied  to  the  surt.ice  ot 
the  half-sp.ice  formevi  In  the  uiiiied  qu. liter  sp.-ue  w.is  considered.  fhe  shear  U'.id- 
ing  was  parallel  to  the  interface  of  the  two  ipi.irter  space  so  th.it  onlv  she.ir 
waves  were  gene r.i t eil . Itt'cause  the  stresses  ,it  the  corners  iiia>  be  singular,  he 
used  a technique  in  which  e.icli  of  tlie  w.ives  geiu’r.ited  .it  the  Ci'riier  was  .iiialy.:ed 
se|iarately  and  then  appropri.itet\'  joined.  fhe  sIumv  stress  on  the  iiitert'.ice  w.is 
found  to  h.ive  a log,ari  t hiiiic  s 1 ngii  1 ,1  r i 1 1 ,it  the  surf. ice  ot'  the  half  space.  This 
singularity  of  stress  m.iy  c.iiise  fr.icture  .it  interf.ice  corners. 

Hrock  and  Achenbach  | hi  1 ] then  extended  the  tr.nisieiit  shear  solution  to  the 
ease  of  loads  normal  to  the  surface  of  the  half-sp.ice  formed  by  the  two  qu.irter- 
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s|>;u'i's.  Mixoii  lu'uiukirv  I'oiul  i t i i)ns  won'  to  l>o  ;i  iu)nii;i]  with  zt'ro 

s(u';ir  stress.  In  this  e;ise,  t lu'  shi'ar  stress  ,)t  the  siirfoee  is  :ero  hee;nise 
of  the  inixeil  hoiiiKl;ir>'  ooiklitit)n  .iiul  heoee  tln're  is  no  siiii;iil;ir  s)ie;ir  stress 
at  tiu'  corner  of  the  Joiiua!  ijuai'ter-spaces.  However,  they  fouiul  that  for  step 
loaJinj;.  a sinjjnlar  shear  stress  is  propaj^ated  at  tiie  StoneJe)'  wa\e  \’elocit)- 
u'hili'  (he  iiori:i:il  si  rt'sses  a r<.'  houiuled. 

Hriick  aiivl  Acheuhach  furtlter  carried  out  an  analysis  on  t iie 

effect  of  an  incident  shear  wave  on  the  extension  of  an  interface  fl.iw 
hetween  two  half-spaces.  .Shear  waves  polarineil  par.illel  to  the  cilite  of  the 
fl.iw  and  perfect  i>-plast  ic  bond  i iii;  were  cons  i dered . Ihc)’  found  tliat  in  case 
of  a const. int -ve  loc  i tv  crack  extension  tlie  incidinit  stress  must  he  .i  step 

w. ive.  The  velocit>'  of  the  crack  is  expressed  in  terms  of  the  amplitude  and 
.injtle  of  the  incident  wave,  the  pri>perties  of  the  halt'  spaces,  .ind  the 
>ield  stress  of  the  bond.  Hus  .inal\sis  shixild  have  an  important  value  in 
interpreting  fiMctures  in  comp<.>s  i t es . 

ko  [ k,' I useil  a teclinique  b.ised  I’n  kirchoff's  method  ol'  retarded  poti-n- 
tials  to  an;tl>’::e  the  intc'r.ict  ii'ti  of  pi. me  waives  with  circul.ir  inclusion.  He 
showeil  that  the  siiycul.ir  stresses  i>ccur  inside  the  cxlinder  which  is  of 
hii;her  im[)edance  th.in  the  surroundiiii;  medium.  He  .also  shi'wcil  the  position 
of  caustics  atul  associated  fold  ovi’r  wave  fri>nts  for  inclusu'ii  wave  specvl-. 
both  hi^lier  .ind  lower  th.in  those  (>('  tin'  m.itrix.  He  then  exteiuled  the 

an. ilysis  to  the  mu  1 1 i - i nc  liis  ion  c.ise  to  study  the  successive  sc.it  t eri  iijjs 
(kl|.  The  stiidx  h.is  .qip  I i c.it  i on  to  fiber  reinforced  composites.  He  showci) 
th.it  the  w.ive  front  am|ilitU(.le  foi’  the  f i is  t -^ener.i  t ion  scattered  w.ive 
alonj;  the  center  line  ot'  the  row  ol'  fibers  attenii.ites  x er>  r.ipiillx,  hecomin>; 
no  more  th.in  10  |H'rcent  of  the  initi.il  x’.iliie  in  the  third  inclusion.  Hi* 
made  a comparison  of  this  case  with  the  laminated  case  and  found  much  more 
rapid  wave  front  aftenu.it  ion  in  the  cylindrical  inclusion  c.isi*.  neterminat  ion 
of  the  stress  behind  the  wave  front  can  he,  in  principle,  obfaini'd  by  this 
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technique,  but  in  practice  such  calculation  will  be  extremely  formidable. 
However,  this  analysis  can  be  further  extended  to  the  case  of  a composite 
containing  inclusions  of  arbitrary  shape. 

Transient  wave  propagation  in  a unidirectional  composite  was  also 
analyzed  by  Haener,  I'uppo  and  Pagan  (1121  in  addition  to  the  exact  analysis 
for  the  steady-state  vibrations.  A circular  filler  surromided  b)-  a circular 
matrix  shell  with  the  outer  surface  confined  was  considered  as  the  model. 

Based  on  elasticity  tlieory,  an  approximate  theory  was  formulated  from 
variation  of  Hamilton's  potentials.  I'he  resulting  4tli-order  differential 
equation  was  tlien  solveil  for  the  transient  state  in  addition  to  tlie  steady 
state.  By  comparing  tlie  numerical  results  of  tlie  apjiroximate  theory  with 
tltose  of  tlte  exact  theory,  they  found  that  for  the  dimensions  considered 
in  the  composite  the  approximate  theory  is  valid  up  to  frequencies  of  10^ 
cycles  per  second.  With  thie  approximate  theory,  it  should  be  possible  to 
consider  mi^re  realistic  boundary  conditions  and  to  solve  more  transient 
problems. 

nrumheller  [OB)  made  an  extensive  investigation  of  the  propagation  of 
elastic  waves  normal  to  the  interface  of  a laminated  composite.  He  used  a com- 
puter code  TIC  based  on  tlie  exact  solution  of  the  response  in  laminates  to  study 
the  transient  and  sinusoidal  wave  propagation.  He  noted  that  the  lowest 
mode  has  an  upper  cutoff  fre<(ueney  that  tends  to  cause  the  dominant  response 
to  occur  in  the  first  mode,  and  he  also  noted  that  different  laminates  can  give 
essentially  identical  first-mode  behavior.  In  a successive  investigation  (l)7|, 
he  developed  this  concept  further  and  showed  how  a transient  analysis  of  one 
one-dimensional  laminate  can  be  used  to  rejiresent  the  response  of  a different 


one-dimensional  laminate  and  very  likely  of  materi.ils  which  have  the  same 
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first-mode  low-frcqucncy  boUavior.  Actually,  he  carried  out  an  analysis 

t 

of  l;uninate  usin^  both  the  original  laminate  properties  and  an  equivalent 
liuninate  which  lias  the  s;une  first  mode  s inusoidal-wave  behavior  but  which 
has  strikingly  different  densities.  He  showed  that  the  linear  response  is 
virtually  identical  for  the  two  materials.  Even  more  importantly,  lie 
compared  the  nonlinear  response  for  the  two  huiiinates  by  assuming  bulk- 
modulus  nonlinearity  and  showed  that  the  agreement  was  still  quite  good. 

The  important  virture  of  this  technique  is  that  it  can  be  extended  to  com- 
plicated constitutive  behavior,  such  as  crushing  and  rate  dependence,  with 
reasonable  confidence  and  that  the  basic  macroscopic  response  obtained 
from  geometric  dispersion  is  properly  maintained. 

lliney,  et  al.  |I19|  first  studied  the  stress-wave  effects  in  inhomo- 
geneous and  porous  earth  materials.  llien  Okubo,  Sve  and  IVhitter  [OlJ 
investigated  the  effect  of  porosit)  to  the  dispersion  of  an  elastic  step 
pulse  in  a three-dimensional  quart r phenolic  composite.  Subsequently, 

Sve  [SU>1  made  an  extensive  investigation  on  transient  elastic  wave  propaga- 
tion in  a porous  laminated  composite.  In  Sve's  analysis,  it  is  assumed  that 
wave  propagates  normal  to  the  laminations  and  that  the  iiorosity  is  raiuiomlv' 
distributed  throughout  one  constituent  and  is  coiiqiosed  of  small  spherical 
voids.  Hie  randomly  distributed  porosity  produces  Rayleigh  scattering, 
and  also  reduces  the  wave  speed  in  a constituent  thereby  affecting  the 
geometric  dispersion.  A dissipative  equation  of  motion  is  developed  for 
the  porous  material  and  used  for  a constituent  of  a composite.  A dispersion 
relation  and  a pulse  solution  are  obtained  to  determine  the  significance  of 
porosity  in  a laminated  composite,  from  liis  study,  Sve  concluded  that  the 
Rayleigh  scattering  produces  a small  damping  effect  in  the  far-field  pulse 
shape  and  small-void  porosity  can  be  adequately  simulated  with  an  effective 
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wave  speed.  Hence,  he  further  concluded  tliat  if  the  pores  are  small  and 
randomly  located  within  a constituent  of  the  composite,  it  may  be  possible 
to  use  a continuum  approach  for  low-frequency  response  calculation  for 
pulse  propagation  in  such  a porous  l;uninatod  composite. 

Sve  (,S44J  also  analyztxl  the  propagation  of  pulse  in  a dissi[)ative 
laminated  composite  by  using  techniques  of  modal  analysis  and  by  consider- 
ing the  comple.x  wave  numbers.  The  analysis  provides  a far-field  dispersive 
solution  that  is  valid  only  near  tlie  head-of-the-pu Ise,  including  spatial 
attenuation.  He  modeled  tlie  effects  of  damping  by  assigning  an  imaginary 
part  to  tl\e  wave  number  and  found  that  the  osc  i 1 lat  ions  about  tlie  steady 
value  of  the  propagating  jHilse  arc  reduced  in  amplitude  and  the  rise  time 
is  increased,  in  comparison  with  the  undamped  case.  In  order  to  obtain 
the  solution  away  from  tl\e  head-of-the-pu  Ise,  it  is  a|iparent  that  other 
techniques,  sucli  as  tiie  method  of  stationary  phase,  would  be  required. 

Chen  and  Curtin  [Cd]  discussed  the  propagation  of  one-dimensional 
acceleration  waves  in  elastic  and  viscoelastic  laminated  composites  by  an 
exact  analysis.  The  composites  that  they  considered  consist  of  a periodic 
array  of  alternating  layers  with  plane  boundaries.  I'hey  made  no  assumptions 
regarding  linearit)-,  and  assumed  only  that  the  materials  have  fading  memory 
when  they  treated  viscoelastic  materials.  ihey  derived  an  expression  for 
the  ;unplitude  of  an  acceleration  wave  ]nopagat  ing  nonnal  to  the  layers. 

Wlien  confined  to  the  junction  points  of  the  cells  this  expression 
has  exactly  the  s;uiic  form  as  that  for  a single  Inonlinear)  viscoelastic 
material.  They  used  this  fact  to  derive  effective  moduli  for  composites. 

In  addition,  they  generalised  their  results,  derived  for  laminates  consist- 
ing of  )ieriodic  cells  of  two  layers,  to  cells  of  N layers.  Tins  theory  is 
applicable  to  ]il;uu'  lottgi  tud  inal  motions  in  three  dimensions. 
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Chen  and  Clifton  [C5]  discussed  the  transient  longitudinal  waves  in 
elastic  and  viscoelastic  bilajninates.  The  bilaminates  consist  of  either 
elastic  or  viscoelastic  laminates  of  uniform  thickness  and  infinite  lateral 
extent.  Wave  propagates  in  the  direction  perpendicular  to  the  laminates. 
Using  the  Laplace  transform  technique  and  Floquet  theory,  they  obtained 
both  wave  front  and  late-time  solutions  for  step  loading  by  means  of 
as>'mptotic  techniques.  Combination  of  wave  front  and  late-time  solutions 
shows  that  the  wave  profiles  consist  of  a rapidly  decaying  precursor 
followed  by  a dispersive  transition  to  an  equilibrium  state.  Tlie  trajisition 
region  in  elastic  bilaminates  is  much  smaller  than  for  viscoelascic  bilami- 
nates. 

Seymour  and  Mortell  [S7]  studied  the  propagation  of  one-dimensional 
longitudinal  pulses  and  weak  shocks  in  nonlinear  elastic  and  viscoelastic 
laminated  composites  by  an  exact  analysis.  The  composite  consists  of 
alternating  laminates,  with  parallel  pl;me  boundaries,  whicii  repeat  period- 
ically. They  assumed  that  the  defoniiat  ions  undergone  by  the  medium  are  of 
small  amplitude  and,  in  the  case  of  the  viscoelastic  composite,  they  arc  of 
high  frequency.  Thus,  in  the  analysis,  simple  wave  solutions  were  superposed 
to  deal  with  nonlinear  waves  which  undergo  reflection  at  the  boundary. 

They  showed  that  by  an  appropriate  choice  of  the  width  of  the  laminates 
the  nonlinear  composite  can  appear  to  first  order  cither  as  a linear-visco- 
elastic or  linear-elastic  material  when  signal  is  read  at  cell  interfaces , and 
that  at  cell  interfaces  the  deformation  in  a periodic-elastic  composite  is 
identical  with  that  in  an  appropriate  nonlinear  viscoelastic  material.  The 
critical  acceleration  is  obtained  for  elastic  composite  below  wliich  no  sliock 
forms.  These  results  are  extended  to  a viscoelastic  composite  chiefly  by 
introducing  a lumj^ed  daii\ping  coefficient  which  is  the  prcxluct  of  the 


86 


attenuation  Jue  to  the  mismatcli  of  impedances  cuid  the  rate-dependence  of 
each  layer.  I'inally,  they  analyzed  the  weak  sliock  for  two  particular  sijjnal 
functions;  11)  wlien  the  signal  function  is  antisymmetric  about  the  wavelet 
at  which  a shock  foniis,  ;ind  12)  when  tiie  shock  is  at  the  front  of  tlie  pulse. 

As  mentioned  earlier.  Barker  [Bfa]  pro]iosed  the  use  of  a viscous- 
dispersion  model  of  the  Maxwell  type  to  represent  geometric  dispersion  in 
composite  materials.  He  used  liis  model  to  analyze  tlie  transient  response 
of  a layered  composite  medium.  The  result  shows  that  the  rise  times  are 
generally  larger  than  those  of  the  exact  solution  and  the  oscillations  about 
the  steady  stresses  caused  by  step  loading  arc  not  presented  in  Maxwell's 
model  solutions.  It  set'ins  that  the  attenuation  is  overpredicted  by  this 
model.  In  fact,  this  technique,  which  is  semi-empirical,  averages  out  all 
tlie  microbehavior  of  the  composite  and  hence  gives  only  a me;ui  stress  history. 

A different  type  of  "viscous"  dispersion  model  was  developed  by 
Bade,  et  al . [B2J.  In  Bade's  model  a general  expansion  of  the  d;imping  terms 
caused  by  both  x and  t derivatives  of  the  velocity  was  iiostulated.  The 
two  lowest  derivatives  of  the  expansion  were  used  to  matcli  experimental  data 
qualitatively.  Tliey  applied  this  viscous  dispersion  model  to  study  the 
transient  pulse  propagation  in  a waveguide-type  l;iminated  composite  medium. 
Comparison  of  dispersion  model  calculation  witli  three-dimensional  experi- 
mental data  shows  good  agreement  except  near  the  region  of  peak  pressure. 

They  also  carried  out  Bloch's  analysis  to  deteniiine  the  effect  of  tlie  variation 
of  properties  in  the  propagation  direction.  It  was  found  that  the  layering 
in  the  propagation  direction  has  little  effect  on  the  dispersion. 


Shea,  Reaugh,  et  al.  [SlOj  developed  a theory  very  simil.ir  to  tliat 
of  Bade.  They  used  essentially  the  same  form  as  tlie  dispersivity  tenn  of 
Bade,  combined  with  a stress  deviator  deiiending  upoti  current  value  of  strain 
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and  the  history  of  strain.  The  model  was  incorporated  into  a computer  code. 

A unidirectional  quartz-[ihenol ic  was  chosen  in  this  program.  The  results  of 
the  computation  obtained  by  imposing  a step-velocity  loading  condition  on 
the  boundary  show  quite  reasonable  agreement  with  experimental  data  except 
that  the  detailed  oscillations  are  not  reproduced  by  the  model. 

Curran,  Secunan  and  Austin  [C15j  proposed  the  use  of  artificial  viscosity 
to  compute  one-dimensional  wave  propagation  in  composite  materials.  They 
proposed  a macroscopic  model  utilizing  the  dispersion  properties  in  the 
existing  hydrocodes  of  IVONDY  and  PUl'l-  types.  In  such  hydrocodes,  the  arti- 
ficial viscosity  was  estimated.  Computations  were  made  for  wave  propagation 
in  epoxy-steel  kiminates,  quartz-epoxy,  and  quartz-phenolic  conqiositcs  for 
which  experimental  and  computational  data  existed.  The  results  are  in  good 
agreement  with  other  models  and  experimental  profiles  for  the  main  wave 
shapes.  It  is  an  original  approach  for  computing  large  cuiqilitude,  one- 
dimensional wave  propagation  in  composite  materials.  However,  the  approach 
is  a macroscopic  one  in  which  only  averaged  values  of  stresses,  displacements 
and  particle  velocities  can  be  predicted. 

Christensen  [ CIO]  discussed  wave  propagation  in  layered  clastic 
composite  on  the  basis  of  dielectric  theory.  Both  periodic  layering  and 
random  layering  were  considered.  In  case  of  periodic  layering,  the  propaga- 
tion became  dispersive  and  in  case  of  random  layering,  it  exhibited 
dissipative  behavior.  Regarding  the  medium  as  an  equivalent  anisotropic 
one,  long  wave  approximation  for  propagation  nonnal  to  layering  was  modeled 
by  using  Botzmun  constitutive  law.  I’he  dispersion  relation  was  expanded 
for  long  wavelength  assumption  and  compared  with  known  results  to  determine 
tiu'  unknowns  in  tlie  relaxation  function.  Results  were  compared  with  known 
experimental  d.ita  and  very  good  agreement  was  observed.  The  theory  was 
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applied  to  study  the  prolilem  of  a pressure  pulse  act  iuj;  on  a half-space  of 
the  layered  medium.  I'he  response  was  found  to  he  identiccil  with  that 
obtained  by  ltei;emier  IHS]  . The  theor>'  proposed  may  be  regarded  as  an 
effective  modulus  model  with  viscoelastic  constitutive  law. 

brumheller  and  Sutherland  [HlOj  proposed  a lattice  model  to  study  tlie 
transient  waves  in  composite  materials.  Based  on  the  assumption  that 
geometric  dispersion  results  mainly  from  the  relatively  periodic  arrange- 
ment of  the  re inforc  iitg  elements  in  tlie  matri.x  rather  than  from  the  precise 
shape  of  each  reinforcing  element,  rliey  developed  a lattice  which  ignores 
the  sliape  of  the  reinforcing  elements  but  preserves  their  periodicit''. 

In  the  application  of  this  lattice  model,  the  co?nposite  was  initially- 
treated  as  a nond isper sive  homogeneous  mixture.  The  effective  or  average 
properties  of  the  mixture  were  determined  bv  either  steady-wave  analysis 
or  appropriate  experiments.  .-V  lattice  was  then  formed  by  redistributing 
tlte  mass  within  the  mixture  to  form  a periodic  structure  of  laminated  plates. 
This  mass  redistribution  was  carried  out  in  a manner  whicly  yielded  a lattice 
with  theoretical  dispersive  characteristics  that  matched  the  measured  dis- 
persive characteristics  of  the  composite.  I'lie  model  was  applied  to  composites, 
consisting  of  a regular  array  of  tungsten  fibers  in  an  aluminum  matrix,  sub- 
jected to  a step  loading.  Also,  flver-tyjH'  impact  experiments  were  I'erfonued 
in  the  plastic  range  of  the  composites.  I'he  agreement  between  exjH'riment 
and  calculation  for  the  arrival  time  and  rise  time  of  the  wave  front  and  for 
the  freiiueney  of  the  ringing  behind  the  wave  front  is  good.  It  seems  that 
for  a wide  range  of  engineering  applications,  this  iiunlel  can  be  used  to  pre- 
dict the  behavior  of  actual  engineering  composites. 

Nayfeli  [Nbl  presentevi  .1  discrete  viscous  lattice  model  to  simulate 
transient  motions  in  el.istic  and  v iscoel.ist  ic  cvwposite.  Viscosity  was 
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introduced  via  dashpot  in  the  lattice.  Integral  transform  teclmiques  were 
used  to  solve  the  lattice  problem.  From  the  numerical  results,  he  found 
that  discretization  in  the  model  introduces  oscillation  about  the  continuum 
solution,  and  viscosity  damps  those  oscillations. 
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IV . L-XPHRINttlNTAL  INVl-STIGA'1'lONS 

Experimental  work  on  the  propagation  of  waves  in  d irect ional ly- 
reinforced  composites  has  been  carried  out  by  many  investigators  b\’  means 
of  ultrasonic  techniques;  pulse  transmission  measurements,  and  other 
tecliniques.  In  this  sections,  sunuuaries  of  tlie  experimental  investigations 
are  given. 

IV. 1 Ultrasonic  Techniques 

Asay,  ct  al.  [AJ4]  measured  sinusoidal-wave  disjiersion  using  ultra- 
sonic techniques  developed  by  Asay  for  homogeneous  materials.  They  tested 
carbon-phenolic  laminates  reinforced  with  layers  of  high-modulus  filaments 
spaced  about  O.b  mm  apart.  The  results  showed  a pronoiuiced  variation  of 
the  phase  velocity  with  the  frequency. 

T.auchcrt  and  Uuzelsu  [T5]  investigated  tlie  dispersion  behavior  of 
plane-hamonic  waves  in  a boron-epoxy  composite  using  ultrasonic 
techniques.  They  determined  the  dependence  of  group  velocity  upon  the  fre- 
quency for  longitudinal  and  transverse  waves  propagating  cither  parallel  or 
perpendicular  to  the  fibers.  It  was  found  that  transverse  waves  propagating 
in  the  direction  of  the  fibers  show  a very  pronounced  dispersive  beluivior, 
and  that  the  group  velocity  increases  with  the  wave  number.  This  type  of 
behavior  is  consistent  with  that  predicted  by  the  continuum  theories. 

Bedford,  Sutherland  and  Linge  lB14|  c.irried  out  extensive  investiga- 
tions on  a fiber-reinforced  composite  of  tungsten  wires  unid irect iona 1 1> 
embedded  in  a bOOl  aluminum-alloy  matrix.  Hie  composite  cliosen  for  tlie 
experiments  was  jirepared  for  two  constituent  ratios,  2..Z  and  dJ.l  percent 
by  volume  of  tungsten,  respectively.  Ultrasonic  experiments  were  con- 
ducted for  plane  compression  waves  propagating  normal  to  the  direction 
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of  fibers  by  using  water-bath  techniques  with  wide-band  transducers,  lixperi- 
mental  dispersion  data  for  the  two  lowest  modes  agree  well  with  the  predictions 
based  on  continuum  theory  of  mixture,  except  at  the  highest  frequencies.  The 
dispersion  data  obtained  demonstrate  that  for  propagation  normal  to  the 
direction  of  the  fibers,  fiber -reinforced  composites  behave  as  wave  filters 

j 

which  selectively  transmit  or  reflect  periodic  waves.  j 

Robinson  and  Leppelmeir  [RIO]  experimentally  studied  j^ormal  propagation  j 

of  ultrasonic  shear  waves  in  a layered  steel-copper  composite.  Samples  were  | 

manufactured  by  diffusion  bonding  alternate  layers  of  steel  and  copper  foil  |, 

together.  The  passing  bands  and  stopping  bands  in  the  frequency  curve  were  1 

observed.  The  e.xperimental  dispersion  fits  existing  theory  [S40,K5]  quite  j 

well.  I 

Rose  and  Mortimer  [Rll]  measured  wave  velocities  in  unidirectional 
graphite-epoxy  shells  with  four  fiber  orientations.  Ultrasonic  techniques 

P 

and  drop  mass  impact  were  both  used  to  generate  longitudinal  pulse.  Good  | 

agreement  with  anisotropic  shell  theory  was  found.  | 

Reynolds  and  Wilkinson  [R8]  studied  experimentally  the  ultrasonic  f 

waves  in  CFRP  (carbon  fiber-reinforced  plastics)  composites.  Experiments 
were  conducted  on  two-  and  three-ply  laminates  having  various  orientation 
and  ordering  of  the  layers.  Results  were  presented  in  terms  of  shear  and 
compression  wave  velocities  versus  propagation  angle  to  the  fiber  direction 

of  uniaxial  materials.  Experimental  data  agree  well  with  theoretical  pre-  | 

dictions.  In  addition,  resin  porosity  was  found  to  have  important  influences  j 


on  t Ivo  |H‘o\'a.k>.i  t i DU  - wavo  vi' li'c  1 1 los  Jt'Ci'iMso  witl>  iiurtMst'  nl 
ot'  rosM'.  porosity. 

Ii'lix  |11|  lU'tormiiiod  o\|U'r  .i  1 1 y pl\aso  voUuitN  .ii\d  a t ( oiiu.i  t ion 

for  sovoral  plastios,  I'lu'  ri'sults  I'or  Unv  amplitiulo  woro  roprosontod 
adoipiatoly  witli  a stand. ird  I inoar  v isoool.ist  io  modt'l  u'ltltin  t iio  limited 
froi|iionoy  ranv;o.  Ho  ompltasiiod  tliat,  fi.'r  oompositos  usinj'.  I lie  pl.isius  .i-, 
m.tfrix,  iii.itori.il  at  t<'iiuat  ii>n  as  wi‘11  .is  j;oomolrio  dispersion  must  Ho  eon 
sidorod  in  analysis  o('  pulse  prop, ij;al  ion.  |{oli,iHlt'  piasliv  liou  should  Ho 
oHt. lined  by  usiny,  the  mo.isuri-d  ph.iso  vt'U'oitN  of  the  uusliuiii  toyollior  willi 
a f r<s|uonoy-<lopondont  ,i 1 1 onu.it  ion  for  st  ri'ss  pul-.o  prop, ly, 1 1 ion  in  ,i 
ooiiiposito  iiiatori.il. 

fh.iny,  tlouohman  .ind  Voo  lt!d|  oonduotisl  ii  1 1 r.f.on  u i <>s<)n,nu<'  mo.isuro 
mont  s of  sound  volooitv  in  thin  oomposito  1 ,im  ina  t os.  I'his  oould  Ho  .in 
import, lilt  .idd  1 1 i on  to  tlio  Kiuk  on  an  oxpor  imonl  .i  1 I \ dillioull  proHlom. 
Moasuroiuont  s wore  made  h\  an  ultr.isoniv  pul'-o  ooho  I oohii  i i|uo . lU'woior,  the 
dosoription  of  the  oxpor  iiiiont  s is  somewhat  skotohv.  It  would  Ho  dill  mult 
to  dupltoato  the  oxpt'rimonts  with  the  information  provided. 

Uoso,  Wany  and  I'osk.i  1UIJ|  oxpor  imoni  .i  I 1 1 ilot  I'liii  i nod  the  wave  surf.uo 
in  a unidirootion.il  yr.iphito  epoxy  pl.ito  Hv  iismy  the  ultrasonio  loohnii|uo. 
file  results  are  in  ayroomont  with  the  oomposito  w,i\  o proj'.iyat  ion  tlioon 
proposed  earlier  In  ^.lny,  ot  ,i  I . |VJ|  .iiul  {'Xti'iulod  Hv  W.iny  ,ind  I'lioKmantol  , 
llsl|.  Itooentlv,  Martin  l^llT'l  prosonlod  a now  mot  luwl  Iv'r  the  m(Msuromonl  ot 
ph.iso  volooili  of  ultrasonio  w,i\os  m ol.isiu-  oompo-.ilos. 


IV.  l loxur.il  lU'soti.inoo  I'oohniHuos 

.Soliiill.’.  ,ind  I'sa  i |SJ|  oHt, lined  oxpor  inuuU  .i  1 vl.ita  on  iiuvluli  ,uul  d.impiny 
ratios  ot'  t'i  Her- re  inforoovi  oompositi's  In  sludiiny  the  I roe  .uid  tvnood  1 1 .uis 
verso  viHr.it  ions  of  o.uitilovor  Ho.ims  m.nlo  ot  the  m,i  i t'r  i .i  I , I ho  oomposito 
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exhibits  imisot ropic , linear  viscoelastic  behavior  when  underjjoinn  small 
oscillations.  The  data  obtained  are  useful  to  designers  concerned  with 
vibrations  and  impact  loading  of  filament-reinforced  canposite  structures. 

Using  vibrating  cantilever  beiun  specimen,  Tauchert  and  Moon  [i'lj  per- 
formed experiments  to  determine  the  dependence  on  freiiuency  of  the  complex 
moduli  of  unidirectional  glass-epoxy  and  boron-epoxy  materials.  Consider- 
able data  were  presented,  i’he  damping  tends  to  increase  with  frequency 
whereas  the  modulus  is  relatively  constant.  At  high  frequencies,  where 
shear  deformation  and  rotary  inertia  are  significant,  the  specimens  are 
characterized  as  Timoshenko  beams.  Using  these  data,  they  compared  predicted 
and  measured  velocity  and  attenuation  of  longitudinal  pulses  in  the  rods. 

They  concluded  that  the  linear  theory  of  viscoelasticity  is  adequate  for 
predicting  velocity  and  attenuation  of  longitudinal  pulses  in  comixjsite 
rods  of  unidirectional  glass-epoxy  and  boron-epoxy  materials  over  a large 
frequency  range. 

iauchert  [Tdj  also  obtained  flexural  resonance  data  for  a class  of 
woven-fabric  composites.  The  data  were  used  to  measure  tlie  complex  moduli, 
and  the  results,  extraj^olated  to  simie  degree  to  extend  into  the  frequency 
range  pertinent  to  stress  wave  propagation,  wore  then  applied  to  velocity 
and  attenuation  prediction.  However,  he  did  not  discuss  dispersion  and 
scattering  phenomena  which  may  have  introduced  sharp  changes  in  the  complex 
modu 1 i . 

IV. 3 Sho^k  Tube  Tests 

I 

Whittier  and  Peck  IWIS)  tested  impulsively  a liuninate  of  carbon-phenolic 
layers  reinforceil  witli  layers  of  high-modulus  filaments  in  a shock  tube.  A 
step-pressure  rise  was  applics.1  to  one  side  of  the  specimen  by  reflection  of 


a jjas  (.IvTiamic  shock  wave.  I'ho  spatial  averajic  of  tlte  froo-surfaco  velocity 
over  the  cross  section  of  the  composite  was  measuroil  with  a cajiacitunce 
j;ajjc.  The  experimental  results  a^^ree  well  with  the  lieaj-of-the  pulse  approxi- 
mation obtained  by  t'eck  and  (.iurtman  [I’^J  . The  shock  tube  test  was  found  to 
be  an  effective  method  to  characterize  the  dispersivity  of  a composite  in 
avoiding  the  compl icat ions  of  nonlinearity. 

Sve  and  Okubo  [SIS)  invest  i)iut tnl  the  rear  surface  velocity  of  three 
huninated  composites  witlt  lamination  angles  of  0°,  4S°,  iX)"  to  the  loading 
direction  by  usini^  capacitance  transducers.  A step  pressure  pulse  was 
imparted  to  the  specimens  usint;  a shock  tube,  with  res|K>7ise  confined  to 
the  linear  elastic  range,  lixper imental  results  confirm  tlte  conclusioits  pre- 
sented in  Sve  and  Whittier's  [Sill  theoretical  studies  which  are:  n)  elastic 
response  is  composed  of  two  coittribut  ions  tltat  travel  at  different  speeds; 

(.ii)  the  dispersion  is  also  different  fiir  each  of  the  contributions  aitvl 
depends  on  the  low- freiiuency  behavior  of  the  pliase  velocities;  (.iii)  far- 
field  average  resiionsesof  0“  and  lamination  c».miposites  are  similar,  while 
near-field  responsc'sof  tltem  arc  not  similar. 

IV. 4 Impact  rests 

Washington  |Wh|  made  near-field  me.isuremeii;  of  stress  waves  propagating 
parallel  to  I'lexiglas-aluminum  himinates.  I'he  laminates  consisted  of  alumi- 
num plates  alternating  with  I’lexiglas  plates.  The  stresses  were  generated 
by  aluminum  or  I’lexiglas-f Iyer  plate,  witlt  tlticknesses  such  as  to  produce 
load  durations  of  0.0  to  0.8  yi^'OC.  The  stresses  were  measured  at  several 
depths,  using  Manganin  wires  embedded  in  epoxy,  located  primarily  at  the 
center  of  each  l;uninate.  The  results  showed  strong  attenuation  in  the 


aluminum  layers,  even  for  O.S  psec.  load,  but  little  or  no  attenuation  in  the 
I’lexiglas.  It  seems  that  at  these  short  propagation  distances  one  is 


ohsorv injj , primarily,  tlu'  transfer  of  energ\’  fri>m  the  aluminum  into  the 
I'lexi.ylas  layer,  as  opposed  to  overall  pulse  attenuation. 

Seaman,  et  al.  | So  | i nvest  ii;ated  experimentally  the  heliavior  of  three- 
dimensional  orthogonal  quartz-phenolic  composites  at  moderate  to  hinl> 
stresses.  I'ive  plate- impact  shots  were  made,  four  with  tltick  flyers  for 
obiainini;  llugoniot  data,  and  one  with  a thiii  fl)'er  for  pulse  attenuation 
data.  Hie  resjuinse  was  measured  with  Manganin  wire  embedded  in  epoxy 

behind  the  specimen.  Ilu.uoniot  data  derived  from  the  flyer  velocity  and 
shock  velocity  wi’rc  closely  fitted  by  a st  ra  i j;ht  - 1 ine  in  shock  velocity 
versus  [i.irticle  velocity,  and  fair  ajtreement  of  these  data  with  the  itaite 
stresses  was  also  obtained.  The  rise-time  data  sujtjtest  that  the  pulse- 
spre-nlitu;  effects  of  geometric  d i spers  i v i t _v  were  overridden  by  the  steej-iini; 
effects  of  non  I i near  i t >■  as  tlie  (leak  stress  was  increased.  The  jiulse  attenua- 
tion results  indicate  both  the  hijthly  attenuative  nature  of  the  material  and 
the  inadequacv  of  those  predictions  which  do  not  account  for  jtcometric  dis- 
pers ion. 

falvit  and  IVatson  ltU|,  and  Sutherland  and  faivit  [S.181  performed  a 
series  of  dynamic  experiments  for  fiber- re i nt'orced  viscoelastic  composites 
subjected  to  uniaxial  pulse.  The  composites  they  tested  were  viscoelastic 
resins  reinforced  with  (,  i ) neoprene  filaments,  1 i i ) rovinjt  jtlass  fibers, 
and  (iii)  nylon  (>(>  fibers.  I'he  transient  pulse  was  initiated  by  a 
mechanical  striker  and  the  particle  velocity  was  measured  by  usinj;  the  so- 
called  faraday  transducer.  The  particle  velocity  records  were  converted 
into  phase  velocities  and  iittenuation  coefficients.  I'.xperimental  results 
were  compared  with  the  predictions  of  effective  modulus  theory  for  the 
three  kinvls  oC  re  i nl  orcement  and  with  the  prcilict  ioivs  of  Hedford  and  Stern's 
mivture  theorx  for  I lu'  latter  two  kinds  v't'  reinforcenu'iit  . Thev  concluded 
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tluit;  (il  it'  the  ot't'octs  of  intornal  j^oomotry  :u'o  small,  tlio  offoctivo 
moiUiUis  thoorv  may  lu'  usoil  to  predict  tUo  Jyn;imic  rosixinso  of  a f ihor- roiii- 
forcixl  matorial;  if  tlu'  offocts  aro  i\ot  small,  tho  predict  ion  of  the  mixture 
theor\’  is  more  accurate;  and  (ii)  for  a constai\t  volume  ratio  of  fiber 
material,  the  devi.it  ion  of  the  tlieoretical  iinxtel  from  the  exper  iiiunit  a 1 
results  increases  as  the  number  of  fibers  per  unit  area  decreases. 

In  .»  series  of  papers,  l,under>;.»n  and  Drumheller  1 1 S , p , 10, 1 1 ) invest  v- 
>;ated  the  dispersion  of  di  lat.it  ionat  stress  w.ives  in  a l;iminati\l  composite 
both  experimentally  .ind  analytically.  The  composite  consisted  of  a number 
of  hi  laminar  plates,  and  the  i>ropaj;at  ion  of  the  wave  was  normal  to  the 
plates.  A flat  flyer-plate  .iccoler.iti\l  by  a compressed  yas  yun  was  usixl 
to  induce  a rect.inyular  stress  pulse.  The  p.irticle  velocity  was  measured 
at  the  opposite  end  of  tiu'  comj>osite  by  an  optical  i nt  ert  eromet  er . The 
experimental  results  were  compared  with  results  obt. lined  by  means  of  com- 
puter \n’oyr;uns  in  which  the  layers  are  explicitly  represented,  .iiul  which 
include  nonline.ir  effects.  In  tlie  first  set  t'f  experiments  |l..s|,  la\ers 
of  epoxy  and  steel  were  impacted  at  stress  levels  on  the  order  of  .i  few 
kilobars.  It  was  found  that  the  peak  cmpilitude  of  the  t r.insm  1 1 1 ixi  stress 
wave  deereasi'il  directl>'  with  the  widtli  of  the  input  stress  |hi1so.  ihe  Ixilk 
of  the  reiluction  of  the  stress  was  .it  tributis.1  to  the  reflections  of  tlie 
stress  wave  at  the  extreme  left  .ind  riyht  bouiularies  of  the  emiiHisite. 
Peboiuliny  of  the  hi  laminate  plates  and  dissipation  of  eneryy  in  the  epoxy 
material  wiM'e  also  tbouyht  to  have  contributtxl  to  the  stress  re.lucf  ion. 
Comparisons  of  the  experimental  results  with  numerical  results  provided 
.III  indication  of  the  necessity  of  cons iileriny  the  nonlinear  behavior  ol 
the  materials. 
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In  siibsequont  experiments  (l.i)]  l;uiui\ates  of  I’MMA  and  stainless  steel 
were  vised.  By  varying  the  number  of  huninates  it  was  concluded  tliat  tliere 
was  no  clvan^e  in  rise  time  for  propan^ttion  distances  beyond  a few  unit 
cells  from  tive  loading;  face.  Another  set  of  experiments  IblO]  was 
carrieil  out  on  the  epoxy-steel  l.uninates  to  i nvest  i j;vite  tlie  resjionse  at 
sliithtly  bibber  stress  where  tlve  effect  of  fracture  would  be  more  apparent. 

It  was  found  that  the  fracture  location  depended  on  the  pulse  durativav,  with 
fair  agreement  between  the  predicted  and  experimentally  determ iiuxl  fracture 
locations.  Subsequently,  l.undergan  and  Prumheller  [LI  1 1 invest  igatisJ  the 
res(Hvnse  of  obi  iquel.v -laminated  composite,  both  experimentally  and  analyti- 
cally. Several  vui.ilytical  models  and  a two-dimensional  wave  propagation 
program  were  used  to  predict  the  transmitted  wave  fonn.  Comparisons  were 
made  between  the  various  models  and  the  ex(H'rimenta]  results.  General  agree- 
ments exist  between  the  predictions  of  various  models.  The  experimentally 
determined  first-signal  velocities  and  final  wave  ;uiiplitudes  agreeil  with  the 
calculations  of  the  models,  liowovcr,  tive  remaining  portion  of  the  experi- 
mental l>’  detennined  stress  waves  exhibited  slower  rise  times  than  did  those 
of  the  calculated  waves.  It  was  speculated  that  the  pivysical  conditions 
existing  at  the  inixit  bound.iry  of  tlie  comjvosite  were  not  being  adequately 
incorporated  into  the  models. 

Berkowit;  and  Cohen  [B-'o]  stiulied  experimentally  high  amplitude  stress- 
wave  propagation  in  an  anisotropic  quartz-phenolic  composite.  I'hey  perfornuxl 
Ihigoniot  and  pulse-attenuation  plate-slap  tests  on  v(uartz-phenol  ic  huninates. 
They  tested  composites  with  l.unination  planes  noimial  to  and  parallel  to  the 
propagation  directiini.  I'hey  concluded  tliat  the  composites  exhibitevl  rate- 
i iiileponilen  t [>l.is  t ic  i t y.  I'hen  thi\v  useil  a plasticitv’  mi>del  and  were  able  to 
let  the  attenuation  maich  the  experiments  with  some  discrepancies  in  pulse 


stiapo.  Iho  miml'or  ot'  tests  was  tvio  tow  to  valulato  tlie  moklol  cone  Uis  i\ o I v . 
blit  tlto  valulity  ot'  the  plastic  model  was  lUmionst rat ed. 

Ueed  [Kll.  and  Miitisoii,  Keed  and  l.under>;an  [Mil]  pert'onued  tlun  pulse 
attenuation  evperiments  on  a c lotli- laminat e nuart ^-phenol ic  composite, 
riiin-pulse  tests  for  propaj;.it  ion  normal  to  the  layers  indicated  that: 

(.  i ) much  more  .ittenuation  occurred  th.in  wmild  he  jnedicted  hy  simpli'  usinjt 
.1  hvdrodyn.uiiic.il  iinvlel  tor  qinirt  ;-phenol  ic;  i,iil  .in  empirical  fit  to  the 
Barker  iBn)  v iscous-dispersion  model  relax.ition  time  >;.iv  e veri  j;ood  agree- 
ment with  the  .Ittenuation  curves.  Me.isurement  s vmi  a t liree-d  imens  ii'ti.i  1 
i|u.irt  ;-phenol  ic  composite  witli  f.iirl>'  hiitli  (Hiri'sit)'  showed  th.it  the  s.uiie 
Barker  model  did  not  pri\lict  .is  l.ir>;e  .ittenuation  as  w.is  found  experiment  .il  l>  . 
It  remaiiuxl  a m. it  ter  of  controversy,  in  three-dimensional  quart phenol ic 
composite  test;.,  .is  to  whether  the  iiicre.ise  in  .it  t enu.it  ion  in  the  three- 
d imens  ion. 1 1 m. iteri.il  w.is  c.iused  In  the  crushin,^.  In  the  j^eometric  dis\'crsion. 
or  h\  a comhin.it  ion  ol  the  two  ef  fects.  Other  ex)H'riments  indic.ited  ytiWKl 
compart  sv'iis  with  both  micromecii.inic.il  models  .ind  cent  iinivim  theories  of  B.ule 
t\pe  [BJ]  , witliout  accvHint  Inj;  for  jvrositv. 

Keed  .ind  Munson  (KS)  .ijt.iln  i lu  est  i i;.it  ist  ■'tress  pulse  .it  tenu.it  ion  in 
clotli- huiilnat  e quart  ^-phenol  Ic  hot  1\  .in.ilv  t ic.il  l>  .nid  experiment. il  ly.  I’he 
composite  was  iiw'deled  .is  .i  ixvxiio^eiieous  \ iscoel.ist  ic  m.iteri.il  (M.ixwell  typel 
with  nonline.ir  stor.i^e  uuhIuIus  .ukI  one  rel.ixat  u>n  time,  .'iidtcious  curve 
fittinij  from  one  \'l.it  e- imp.ict  test  pix'diicixl  the  time  constant,  and  suhse- 
quent  comparisons  .inalitic.il  predictions  with  experiiiient.il  results  were 
Based  on  this  v.ilite.  lxperiment.il  results  sliowed  m.irkixl  dispersive  sprea'in.c 
of  the  w.ive  profiles  .ind  vein  st  roiii;  .it  t enu.it  ion  of  stress,  hut  no  evivience 
of  el.ist  i c-pl.ist  ic  effects.  Model  predictions  .i.^reevl  very  well  with  the 
miMsiired  profiles.  .\t  t enu.it  ion  \uevl  ic  t umi  w.is  also  quite  >;ood.  The  vl.it.i 
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and  analysis  siij;gostoil  that  tho  viscoelastic  nature  ot'  t i\e  phenolic  matrix 
plays  a dominant  role  in  the  attenuation  behavior. 

Michael,  Christman  and  Isbell  (M.^l  studied  experimentally  the  pulse 
propajiation  in  composite  materials.  I'he  pulse  w.is  jtoner.itixl  by  t'lyer-plate 
impact  .ind  particle  velociti  w.is  me.isured  by  usinj;  velocity  .k;aj;e.  I'he 
materials  tested  were  iXiMRAb  1 ci'iiiposite  .uhI  huiiinated  ijuar  t ;-)>heno  1 i c . 
Measurements  on  tXWRAb  1,  with  pulse  proj'.icat  ion  nonii.il  to  fibers,  sliowed 
that  the  structure  of  the  pulses  was  approx  im.i  t e I y t lie  s.ime  for  linear  .ind 
nonlinear  response  levels.  Tests  on  l.unin.ited  qu.irt  phenol  ic  i iwl  i c.it  i\l 
that  the  response  to  step  loads  was  somewhat  il  i sp-'rs  i\ <>.  It  .ip|u'.ireii  th.it 
rate  dependence  .ind/or  pl.isticiti  pl.iied  .i  role  .is  larj;e  or  l.iijjer  th.iu 
that  played  b>’  the  geometric  disperion. 

Holmes  .ind  Tsou  iHlS]  invest  igatovl  the  ste.uly  slu'ck  w.i\  es  in  fiber- 
reinforced  composite  materials.  The  shock  w.ive  was  getier.it  ..xi  b\  .i  pl.in.ir 
impact  of  a flyer  plate  and  prop.igated  .i  tong  the  direction  of  the  fibers. 

The  siuiiple  used  in  experiments  was  m.ide  of  unid i rect  ion.i  1 .iluminum  t ibers 
cast  in  an  epoxy  matrix,  both  shock-wave  velocitv  .ind  free-surf.ice  velociti 
were  me.isured  b>'  me. ms  of  optical  t ecliii  ii|ues.  The  shock  front  in  the  com- 
posite was  found  to  be  steady.  The  results  of  the  me.isured  Itugoiiiot  .ilso 
gave  a satisfactory  comparison  with  those  obtained  from  an  analysis  propost\l 
earlier  by  Tsou  and  t'luni  ITI.^].  The  experimental  justification  of  steadi 
shock  front  is  quite  significant  since  the  assumption  of  steadv  shock  can 
bo  utilized  to  solve  problems  involi ing  composites  of  v.irious  configurations. 

Barker,  l.iinderg.m,  idien  and  i>urt  in  [Bl'l  studied  experimentally  the 
shock  waves  and  acceleration  w.ives  in  lamiiiatixl  coiin'osit  es.  To  verify  the 
prediction  of  Barker's  model  IIU'],  they  consulerixl  .i  ci.'Hiipv's  1 1 e consisting 
of  altern.it  mg  l.iyers  of  I'ol  xiiiet  hi  I -met  li.icrv  l.it  e and  .iluminum.  Not  onlv  vlul 
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the  oxpor imont s ostablisli  tl\o  existence  ot'  shock  waves  in  low  stress  rejiion, 
but  the  nonlinear  viscoelastic  .iiodel  predicted  quite  accurately  their  gross 
structure.  Also,  to  verify  Chen  and  Curtin's  results  IC4] , they  considered 
a composite  of  alternating  laxers  of  polxmethyl-methacrylate  ;ind  crt^ss-cut 
quart;.  I'hey  generated  in  this  composite  acceleration  waves  witi\  iniKit 
.unplitudes  both  above  and  below  the  critical  ;unplitude  a^, , and  in  both  of 
these  instajices  the  results  predicted  by  the  t]ieor>'  were  verified.  These 
experiments  wore  the  first  to  produce,  in  composites,  steady  shock  waves  at 
svich  low  stresses,  and  the  first  to  display  acceleration  waves  in  composite 
materials. 

1 \' . S Ij- a^t u£e  J'e st  s 

Warnica  and  Charest  1W5J  tested  by  means  of  plate  impact  himinated 
quartz-phenolic  with  the  layers  inclined  at  ‘)0°,  0“,  45°  to  the  direction 
of  wave  propagation.  i\ilse  durations  were  1 to  d psec.  A definite  spall 
threshold  was  obtained  only  for  the  0°  composites,  wliere  spallatioit  occurred 
due  to  deUuninat  ion.  In  the  cases  of  45°  and  90°  angle  inclinations,  dehunina- 
tion  occurred  before  spallation.  In  the  0°  case,  a significant  difference 
in  stress  was  observed  between  the  onset  of  microscopic  cracking  and 
m;icroscopi c spal  1 . 

tireen,  Babcock  and  Perkins  ItlOj  invest igatexi  experimentally  the  problem 
of  degradation  of  mechanical  properties  of  composites.  I'hey  found  that 
huninated  i\uart  ;-phenol  ic  composite  testeil  with  momentum  traps  ujx  to  9 k-bar 
hud  no  degradation  in  compressive  [properties.  I'or  c;gT)on-i'heno  I ic  laminates, 
stress  up  to  T”  k-bar  had  no  effect  hut  at  8 k-bar  a strength  decrease  of 
dO-JPa  was  fi'und.  In  these  experiments,  it  was  felt  that  not  all  extranoxnis 
d;ijnage  was  su|ppresstxl  and  the  specimens  might  have  felt  siwe  tension  as  well 


as  compression. 
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I'orklivs,  HatK'ook,  Stiiorlock  and  .lonos  [I's)  conUuctiVi  an  oxj'cr  imont  a I 
stiuiy  ot'  the  et'feet  of  impulsive  prelo.klinj;s  on  the  behavior  of  nuart:- 
phenolie  composite.  iTiiar  t ;-plienol  ic  himinates  wore  impact-tested  with 
momentum  traps.  lair  comjiression  strength  .md  stress  wave  loadiiij;  parallel 
to  the  lamination,  dojiradation  occurred  at  both  room  and  elev.ited  tempera- 
tures. I'or  tension  strength  and  stress  wave  loading  normal  to  the  lamin.it  ion, 
degradation  occurred  .it  room  temperature.  Uniaxial  strain  and  unia.xi.il 
stress  preloadings  were  performed  on  iiuart r-phenol ic  to  assess  the  effect 
of  an  impulsive  prelo.id  on  tlie  subsequent  uniaxial  stress  behavior  of  a com- 
posite material.  In  these  tests  oliseiTable  damage  w.is  limited  to  the  phenolic 
matrix.  As  .i  result,  juist -jirelo.id  ex]ieriments  on  the  specimen  under  uni.Lxial 
stress  correlated  the  deiiendence  of  the  fracture  strength  on  the  extent  of 
the  d. image  to  the  jihenol  ic  matrix  as  well  as  the  direction  of  the  preload 
rel.itive  to  the  fiber  l.iy-up. 

Schuster  .ind  Retxl  [S.i]  investigated  the  fracture  behavior  of  boron- 
aluminum  composite  m.itei  i.ils  by  means  of  shock  loading.  The  ciniiposite  was 
formed  from  boron  fil.unents  and  aluminum  matrix.  The  boron  filiuwents  ran 
perpendicular  to  the  direction  of  stress-wave  loading  in  two  orthogon.il 
directions  in  the  pl.uie  of  the  siH'cimen.  Two  tvpes  of  aluminum  matrix  were 
used:  one  was  brazed  and  the  other  w.is  diffusion-bonded  to  the  fil.unents. 

Roth  composites  suffered  filimient  cracking  thouglvt  to  result  from  the  cixn- 
pressive  wave  passage.  The  sp.il  1 cracks  formed  on  \ilane  parallel  to  the  free 
surface  of  the  composite. 

Reed  and  Schuster  [Re]  further  studied  the  f i l.uiient  fracture  and  post- 
impact strength  of  boron-aluminum  composites.  The.v  tested  specimens  with 
momentum  tr.ips  to  suppress  sp.ill  .iiul  then  measured  tensile  strength  in  the 
plane  of  the  specimen.  fhey  found  the  degrad.it  ion  began  to  occur  at  a fairlv 
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wol  t'lvor  votucitv  t h.iii  tli.it  ic'i[u  i roil  to  *111.111  .1  him  i mini. 

Rot.iiiioil  stroiii;th  as  low  as  0110  thii\l  ol’  tho  iioniial  st  roiijtt  li  ’.  .is  oht.iii’.oii 
, at  tiio  hik^liost  t'lyor  lolocitios.  It  .iiipoaroil  th.it  t tu-  primari’  moJo  of 

‘ il.uii.ino  in  those  s}'ocinions  was  tho  iloboiulini:  of  tho  fihors  from  tho  siirrouiul- 

ini;  matrix.  I'ho  iiu  ost  ijtat  ion  sop.iratoil  tho  offoot-*  of  oompross  ivo  shook 
loading;  from  offoots  of  sp.illation.  rr.iotnro  of  tho  fil.uiionts  was  shown 
to  ooour  diirini;  tho  initi.il  p.iss.i.oo  of  tho  lo.idiiii;  w.ivo. 

B.irboo,  So.uii.in  .iiul  I'rowilson  [H.s]  m.ido  .in  o\por imonta  1 stnd>  of  d>n.imio 
fr.ioturo  of  .i  qu.irt  :-phonol  io  oomposito  .ind  a si  I ioa-plionol  io  ooinposito. 

I'ho  i|ii.irt  :-phonol  io  oomiHisito  oomposod  of  unidirootion.il  0 loso  1 y- p.iokod 
oarns  iv.is  tostisJ  with  w.ivo  propai;.it  ion  p.ir.il  lol  to  tho  y.irns.  Tho'  found 
th.it  or.ioks  both  parallol  .itul  porpond  iou  l.ir  to  tho  >.irns  oooiirrod  in  tho 
m.itri.x.  Tho  oimiposito  w.is  .1 1 so  tostod  with  momontum  t r.ips  and  thon  or.ioks 
woro  found  only  p.ir.illoi  to  tho  t'ibors.  thaok  w.ts  oomprossivo  in  origin. 

Tor  .1  l.unin.itod  s 1 1 io.i-phonol  io  oomposito.  no  oixiipross  i\ o dam.igo  ooourrod 
whon  tostod  witli  momontum  traps,  but  it  was  found  th.it  whon  tostixl  without 
momontum  traps  .1  sp.ill  lino  was  formod  rouglily  parallol  to  tho  froo  surf.ioo. 
Tho  sp.il  I took  pi. 100  partl>'  b)  dobonding  botwoon  tho  l.iyors  and  p.irt  ly  by 
fr.ioturo  though  tho  l,i>ors,  .is  w.is  roijuirod  to  form  .1  spall  piano  parallol 
to  tho  froo  surfaoi . 

Gorborioh  [tl.'']  oxpor imont .1 1 1 >■  an.ilyoiwl  tho  fr.ioturo  boh.iv  ior  of  ,i  oom- 
posito with  duotilo  fibers.  Aluminum  m.itrix  oompositos  roinforood  with 
unidiroot  ional  stainloss-stool  wiros  ot'  -150  ksi  tonsilo  strength  woro 
an.ilyzod  in  terms  of  strength,  stiffness,  .ind  fr.ioturo  oharaotorist  ios 
along  and  aoross  tho  fiber  arr.iv.  Tonsilo  strengths  .is  high  as  JPO  ksi  woro 


obtainoil  in  tho  fiber  dirootion  while  ISO  ksi  tr.insvorso  strengths  wore  ro 
oordo'd.  Tong i tud ina  1 strength  and  moduli  .igroixl  ro.ison.ibl\-  well  with  rule 
of  mixture  in'odictions  up  to  lO'k  fiber  by  volume.  Intorfaoial  bonding  w.is 


suft’icic'ut  to  provide  transverse  stiffness  enlianccment  but  did  not  contribute 
to  transverse  strength.  Crack  propagation  across  fibers  was  found  to  be  con- 
trolled by  the  very  ductile  higli  strength  fibers.  Crack.  i)ropagation  betwceti 
fibers  was  controlled  by  fiber  spacing.  Critical  stress  intensities  for 
transverse  crack  propagation  were  reported  to  be  as  low  as  one  tliird  tliat 
of  the  matrix. 

Sierakowski,  et  al.  [S12J  [iresentcd  results  of  an  cxj)er imcntal  progr.ain 
systematically  evaluating  the  defoniiation  and  fracture  of  steel  wire-reinforced 
epoxy  composite  systems.  The  program  involved  mecb.an ical  testing  in  the  strain 
rate  range  10  ^ to  lO'^/sec.  and  impact  testing  using  massive  elastic  targets 
at  strain  rate  approximately  loV^t'c.  Specific  results  included  static  and 
dynamic  properties,  strain  rate  sensitivity,  information  on  tlte  nature  and 
character  of  dyi\;imic  fracture,  influence  of  specimen  geometry  and  reinforce- 
ment spacing,  etc.  I•ilrtller,  they  projiosed  a simplified  energy  criterion  for 
predicting  failure  modes  and  critical  velocities  for  composite  specimens  witli 
a brittle  matrix . 

Cohen  and  Berkowitz  [C12]  also  studied  experimentally  dyn;unic  fracture 
of  a quartZ'[)henol  ic  composite  uitder  stress-wave  loading  in  uniaxial  strain. 
They  carried  out  thin-flyer  tests  using  5 mil  and  15  mil  Mylar  flyers 
impacting  0.25  inch  thick  composites.  fhey  found  significant  difference 
between  microsco|iic  and  gross  spall,  and  tluit  spall  occurred  by  del:iminat  ion. 
They  found  secondary  cracks  perpendicular  to  tlte  impact  faces  which  were 
most  intense  near  the  front  and  rear  faces  of  tl\e  material  and  could  have 
been  caused  by  compression  wave  or  late-time  flexure  of  the  specimens.  They 
also  found  that  a constant  stress  tensile  fracture  criteria  apjilied  in  a low 
impulse  region  while  a rate  process  criteria  aiij'licd  in  a Itigh  imjiulse 


region. 


I 
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lloovor  aiul  (iiu'ss  ^U*si^iu\l  oxporimonts  ti>  moasure  t lu*  (.lyi\amic 

I'ractiiro  tonj;(inc'ss  aiui  the  work-of-rracliiro  par;unotor  'i  as  t'uiietious 

of  till'  rate  of  loaJin^;.  k||  is  the  ciMtioal  stress  intensity  factor  at  which 
crack  initiation  occurs  aiul  Yj-j,  is  o measure  of  the  enerjjy  ahsorhed  as  a 
crack  initiates  and  propagates  t hrouy;h  the  matei'ial.  rite  materials  tested 
were  carhon-carhon  composites  havinij  almost  exclusive  application  in  space 
industry.  The  data  oht. lined  slioiild  he  of  jiractical  use,  and  the  techninue 
described  ma\’  he  useful  for  the  dynamic  testing;  of  otlier  materials. 

Drumheller  109)  invest  ir.at ed  the  effect  of  delumdinit  "o  stress  wave 
propagation  in  comiiosite  materials  both  analytically  and  experimentally. 
lUfferin.y  from  Sve's  approach  l-S-lJl  h\’  i nt rivluc in>;  dehondin>j,  he  studied 
the  phase  velocit>’  behavior  in  the  limit  of  zero  wave  number,  lie  found  that  | 

i 

in  the  limit  tliere  are  three  values  of  phase  velocity  instead  of  two  which 

were  found  lyv  Sve  without  account  inp,  tlie  effect  of  debondinp.  lie  confiriiuxl 

these  three  values  by  exi'erinunts  on  a laminated  composite  consistinp  of 

stacks  of  stainless  steel  and  po  1 yaiiet  h>' 1 -met  hacr\  la  t e jilates.  To  obtain 

povxt  apreement  lietween  theory  and  experiment  for  tlie  third  velocity,  it  was  ' 

found  necessarx-  to  employ  h i ph-pressure  iiotvmielliy  1-met  hacrylate  data. 

brumheller  and  Norwood  Ibl.x),  and  brumheller  and  l.underpan  [IHdJ 
further  studied  the  debond inp  effect  on  the  behavior  of  stress  waves  in  com- 
posite materials  theoretically  and  experimentally.  I'he  theoretical  investi- 
pation  was  focused  on  the  problem  of  a deboiuled  comiiosite,  a situation 
commonly  observeil  in  flyer-pl.ite  impact  exj'eriments.  In  contrast  to  the 
fully  bondetl  case,  an  ailditional  stable  wave  mode,  tliey  suppesti'd,  may 
propapate  in  thi'  composite.  Tlius,  anotlier  boundary  coinlition  is  required 


to  obtain  the  solution.  They  postulated  such  a condition.  The  essential  con- 
cept was  to  introiluce  an  interf.ice  warpinp  stiffness  rel.it  inp  the  difference 


1 J5 

hetwi'cii  the  const i tiiciit  stress  aiui  the  average  iiitert'ace  stress,  witl)  ilie 
difference  between  the  corresjiond ing  displacements,  'llie  theory  was  applied 
to  a configuration  modeling  tlie  flyer-jilate  experiment.  l,x]>er  imont  s were 
then  conducted  oi\  a composite  formed  from  laminates  i>f  t>l)(>l-'lt)  aluminum  and 
]>olyinethy  1 -met hacry la t e.  Pressure  were  sufficiently  high  to  cause  dehonding. 
Some  difficulties  arose  in  comiiaring  theory  and  exper imettt s . The  theory 
was  based  on  the  assumption  of  total  dehonding.  However,  tlie  experimental 
results  indicated  significant  bond  strength.  After  they  modified  the 
theoretical  results  by  an  emi>irical  correction  to  calculated  wave  speed, 
ipiite  reasoitable  agreement  was  obtained.  The  tlteor)-  presented  is  unsatis- 
factory in  several  respects  as  they  admitted. 

Baldwin  and  Sierakowski  iBl]  investigated  the  uniaxial  static  and 
dynamic  fracture  characteristics  of  a composite  material  consisting  of  .an 
aluminum  matrix,  A-!;>  casting  alloy,  and  stainless  steel  fibers,  type  .101. 
l)yn;unic  compression  specimens  h^avled  par.tllel  to  the  fiber  direct  ioti  failed 
by  fil.uiient  buckling,  white  loading  transverse  to  the  fihmients  produced 
fiber  matrix  dobonding.  Ihe  com\iosite  system  tested  vtid  not  exhibit  any 
rate-sensitivity  in  its  failure  characteristics. 
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V . c:(>Ml'AR  1 SONS  AN!)  U1h:OMMI:NI)AT  1 ONS 

V . I C’ ompa risoiis  o f_  V ;n-  i mis  IJi foj-  i c s 

riio  of  foct  i voiK'SS  aiul  tlio  appl  icatu  1 ity  of  a tlioory  Jopoiui  on  t ho  kind 
of  solution  ono  is  Uiokinj;  for.  A thoory  may  ho  vory  offootivo  in  fiiulin^i 
tho  solution  to  sonio  aspocts  of  wavo  [iropajtat  i on  hut  is  totally  unsuitahlo 
for  fiiulinx  tho  solution  to  othor  aspocts  of  tho  wavo  phonomona.  Tho  applica- 
bility of  a thoory,  in  most  instances,  roflocts  tho  assiunpt ions  mado  in  tho 
thoory.  In  comiiarins  tho  ajipl  icahi  1 ity  of  various  thoorios,  wo  will  discuss 
some  of  tho  hotter  ostahlishod  thoorios  as  to  what  limitations  aiivl  assum|Uions 
aro  imposed  on  the  theory,  and  the  merits  and  drawbacks  of  tho  thoory.  In 
particular,  wo  will  look  at  the  following  (piestions:  Is  the  theory  applicable 
to  nonlinear  composites  and  non-periodic  composites?  Clan  tho  theory  be  easily 
used  for  solvinjt  the  transient  |U'oblems?  How  jjood  is  tho  approximation  ami 
how  difficult  is  it  to  include  the  hiphor  order  tonns  to  >;ot  a bettor  approx i- 
mat ion? 


V.  1 . 1 T.ffoy‘_t  ivo  modulus  theories 

In  the  effective  moilulus  theories,  a cominisito  is  replaced  by  an  aniso- 
tro\HC  linear  solid  whoso  static  responses  aro  oi|uivalont  to  tho  macroscopic 
responses  of  the  compositi'.  Ttu'  relations  between  the  macroscoi'ic  stresses 
and  strains  are  doriveil  solely  on  tho  bases  ot  static  loailinj;.  IVith  tiu' 
effective  moilulus  theories,  the  wave  speeds  in  a comjiosite  aro  indoi'ondont 
of  tho  f roiiuonc i os . Theroforo,  tho  disi'crsion  phenomena  observed  in  tho 

experiments  are  not  predicted  by  tho  offootivo  moilulus  thoorios.  (lonsoipient  1 > , 
the  theories  are  totally  inailequate  for  predicting  the  dynamic  rosj'onses  of  a 
comiios  i t e. 

I'or  the  static  response  of  a composite,  tho  off'oettvo  iiu'ilulus  theories 
;ire  quite  adequate.  The  theories  can  be,  and  have  boon,  a)qilioTi  to  nonlinear 
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comjiosite  materials  such  as  elastic-plastic  composites.  There  is  no  reason 
why  the  theories  cannot  he  extended  to  non-(>eri<.>d ic  composites,  although  the 
analyses  would  become  complicated.  In  most  effective  modulus  theories,  the 
ratio  of  the  micro-dimension  to  tlie  macro-dimension  is  not  critical.  The 
tlteories  are  exact  rej;ardless  of  the  ratio. 

V.  l.d  1!  f f ec  t i ve  s t i f f tie  s s rheotnes 

There  are  two  crucial  steps  in  the  effective  stiffness  tlieories.  We 
will  use  the  bilaminates  as  an  example.  I'irstly,  the  displacements  within 
each  layer  is  expressed  in  a Taylor  series  (or  otlier  \io lynomi als)  about  the 

midplane  of  the  layer.  The  coefficients  ot'  the  series  are  therefore  defined  1 

I 

only  at  the  discrete  points,  the  midpl.iiu's.  The  continuity  of  displacement  ; 

( 

at  the  layer  interface  yields  a finite  difference  equation.  With  the  assumed  i 

I 

displacement  field,  one  obtains  the  strains,  strain  energy  and  kinetic  eiier)’y  r 

I 

in  each  layer.  The  second  important  step  is  the  smoothinji  operation  in  wliich  j 

the  functions  (ireviouslv  defined  only  on  discrete  points  are  extended  to  f 

I 

defined  for  all  points.  This  is  .iccomp  1 i shed  by  takiny;  a weighted  average  of  ) 

the  strain  energy  in  the  layers.  Assuming  the  smallness  of  the  layer  thickness, 
the  continuity  condition  is  rewritten  in  a differential  form  and  lienee  a con- 

I 

tinuum  theory  is  developed.  i 

i 

Application  of  Hamilton's  princii'le  is  to  obtain  the  best  approximate  J 

solution  for  the  assumed  displacement  field.  This  results  in  a s)';.tem  of 
differential  equations  for  the  d i sji lacement  in  the  composite. 

The  accuracy  of  the  theory  depends  on  the  series  ex]ninsion  of  the  dis- 
placement in  the  layer  which  in  turn  depends  on  the  thickness  of  the  layer. 

Higher  order  approximat ion  can  be  obtained  but  one  has  to  do  the  entire 
derivation  from  the  very  bejjinnin^;.  Applic.it  ion  to  transient  [iroblems  is  .i 


108 


matter  of  solving  the  system  of  ili  fferetit  ial  ec|iiations  ami  there  are  12  of 
them  for  the  lowest  order  of  approximation.  Sinee  the  basic  elements  in 
the  theor>'  are  the  assumption  of  the  displacement  field,  smoothing  operation, 

I 

and  applicatiini  of  Hamilton's  jirinciple,  the  theory  can  he,  in  principle, 
extended  to  nonlinear  system.  However,  solving  12  nonlinear  differential 
equations  is  not  a simple  matter. 

Application  of  the  theory  to  non-]ier i od i c composites  is  not  possible. 

However,  the  theory  can  he  apiilied  to  jier iodi ca  1 1 v layered  composites  in 
which  each  unit  cell  cimsists  of  more  than  two  dissimilar  la>ers. 

V.  I . .>  I'heory  ^f  iiiK'ract  in^j  cont  inua 

In  the  theory  of  interacting  cont inua,  not  only  the  displacement,  hut 
the  St  I'ess  also  are  exj'anded  in  a power  series  in  sjiace  v.ariahles  about  the 
midplane  of  each  layer.  Insteail  of  using  the  Hamilton's  [irinciple  to  obtain 
approximate  differential  equations  for  the  coefficients  of  the  power  series, 
the  entire  power  series  for  the  d i spl  .icement  and  the  stress  are  substituted  ' 
into  tfie  equation  of  motion  and  the  constitutive  ei|uation  to  oluain  differen- 
tial-recurrence relations  for  the  coefficients  of  the  power  .series.  Therefore, 
coefficients  of  the  higlier  order  power  are  exiiresseil  in  terms  of  the  coeffi- 
cients of  the  lowest  order  power  in  the  stress  and  the  di splaceiiK'nt  . The 
original  power  series  expansions  are  rewritten  in  series  of  time  differen- 
ti.-ition  of  the  coefficients  of  the  lowest  order  expansions.  ;\p\il  icat  ion  ot 
the  continuity  in  stress  and  displacement  at  the  layer  interfaces  yields 
four  differential-difference  equations  for  the  stress  and  displacement 
at  the  midplanes  of  the  two-layered  com]iosite.  The  finite  difference 
in  the  d i f ferent  i a 1 -<.1  i f ference  equation  contains  the  thickness  of  a 
unit  cell  in  the  composite.  The  assumption  of  smallness  of 


tho  tlui'ki\oss  allows  oiio  to  oxpaiul  tho  il  i t't'oroiu'i'  oijuation  in  Taylor  series 
and  one  lias  four  partial-differential  ei)uations  of  infinite  order. 

As  it  turns  out,  eaeli  differential  operator  is  aecompaiiied  h>'  the  param- 
eter i whieli  is  the  r.itio  of  the  thiekness  of  a unit  cell  to  the  typical 
macro-dimension  iJ.  Moreover,  the  differential  ojH'rators  are  all  of  even 
order.s.  ihe  ;ero  order  approximation  yields  the  effective  modulus  theory. 

The  first  order  theory  (which  is  first  order  in  and  could  have  been  called 

second  oi\ler)  \ ields  .1  result  siimewhat  better  than  the  effective  stiffness 
theory.  This  is  not  surprisiiii;  since  the  first  order  theory  takes  into 

account  the  continuity  in  displacement  and  stress  at  the  layer  interface 
where.is  the  first  order  theor>  in  the  effective  stiffness  theory  takes  into 
account  the  continuitv  in  displacement  only.  Presumably,  one  could  develop  a 
second  order  effective  stiffness  tlieor>’  which  will  be  equivalent  to  the  first 
order  theory  of  interact  in.c  continua. 

The  theory  of  interacting;  continua  has  m.iny  useful  features,  firstly, 
the  improvement  of  the  accur.icv-  b\’  includini;  the  hifjher  order  terms  is  easier 
to  accomplish,  for  the  bi laminates,  in  particular,  we  have  an  explicit 
closed-form  expression  for  any  order  of  accur.icy  we  wish  to  obtain  by  usin.c 
I'.q.  (.->2)  derived  on  p.  l.^.  .Secondly,  this  theory  can  be  useii  either  for  .1 
steady  state  vibration  of  the  composite  or  for  .1  transient  response  of  the 
composite.  Thirdly,  at  least  for  the  bi  lamin.ites . the  exact  frequenc>  equa- 
tion vs  readily  recovered  from  the  theory. 

There  is  no  reason  wh.v  the  theory  cannot  be  extended  to  nonlinear 
composites.  However,  it  would  cert.iinly  become  unwieldy  if  an  order  of  hi({her 
than  tho  first  one  is  attempted.  As  in  the  effective  stiffness  theory, 
application  to  non- peri  oil  ic  composites  is  not  possible. 


V.1.4  Mixture  theories 


In  the  mixture  theories,  the  constituents  are  assumed  to  coexist  and 
are  allowed  to  have  an  independent  motion,  even  though  perfect  bonding  between 
the  constituents  is  understood.  The  assiunption  on  tlie  relative  motion  of 
the  constituents  is  the  crux  of  the  theories  and  is  also  the  main  roadblock 
for  improving  the  accuracy  of  the  theories.  Although  satisfactory  results 
are  obtained  for  liarmonic  waves  in  a simple  bilamin.ite,  the  problem  of  assuming 
a suitable  interaction  between  the  laminates  still  remains. 

Mixture  theories  can  be  applied  to  both  harmonic  waves  and  transient 
waves.  Extension  to  nonlinear  composites  is  possible  but  again  tlu'  diff.cult> 
lies  in  defining  the  interaction  between  the  constituents. 

V . 1 . 5 0 1 h e ^t  h cojj  vj. 

There  are  other  less  widely  used  theories  wliich  we  will  not  co?npare  in 
detail  here.  However,  a few  words  .about  the  Eloipiet  theory  and  the  varia- 
tional techniques  are  in  order. 

I'loquet  theory  applies  to  a system  of  linear  differential  eipiation  with 
coefficients  which  are  periodic  functions.  Therefore,  it  can  be  readily 
applied  to  composites  whose  constituents  are  arranged  periodically.  The  most 
one  can  get  from  apj-ilying  the  i loquet  theory  to  a composite  is  the  form  of  a 
steady  wave  train  jiropagating  in  the  composite  and  the  relation  between  the 
frequency  and  the  wave  number  of  the  wave  train,  for  a transient  wave,  one 
has  to  superimpose  w.ive  trains  of  all  frequencies  to  achieve  the  specified 
initial  and  boundary  conditions.  Therefore,  I'loquet  theory  is  most  convenient 
Tor  solving  steady  state  motion  in  linear,  periodic  composites.  It  is  not 
suitable  for  transient  wave  motions  except  when  only  an  asvTuplotic  solution 
for  large  time  is  desired. 


Ill 
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111  variational  techniquos,  tho  governing  di  t'foront  ia  I equations  for 
wave  motions  in  a composite  are  written  in  a form  of  integrals  such  tliat, 
by  setting  the  first  variation  of  tlie  integrals  to  zero,  one  recovers  tl\e 
governing  equations.  Tlie  form  of  the  integrals  is  not  unique,  and  hence 
some  variational  methods  .ire  more  effective  than  others.  Unlike  the  I'loquet 
theory,  variational  methods  .ire  not  limited  to  linear  system;.,  altliough  in 
literatures  they  are  applievl  mo'^tlv  to  linear  composites,  liven  for  linear 
composites,  variational  techniques  *re  used  for  solving  steady-state  vibra- 
tions only.  Applications  to  tr.insient  problems  are  possible  theoretically 
but  impractical.  .Some  of  the  ».li  fficul  t ies  in  applying  the  variational 
principles  are  the  .issumption  of  the  test  functions  and  how  one  chooses  the 
next  improved  test  functions.  Therefore,  the  choice  of  test  functions 
requires  subjective  judgements  except  in  simple  cases  where  the  judgements 
a re  trivial. 

V.l.o  Concluding  rem.irks 

Although  v.irious  theories  are  .ivailable  for  t re.it  ing  linear,  periodic 
composites,  they  c.in  be  equivalent  to  one  another  if  a proper  assiunpt  ion  is 
m.ide  in  each  theorv.  An  ex.imiile  was  given  by  llegemier  (1111)  in  which  he 
showed  the  evpi i va 1 ence  of  the  modified  first  order  theory  of  interacting 
continu.i  to  the  bin.iry  mixture  theory  for  waves  propagating  normal  to  the 
layering  of  a composite.  By  combining  the  equations  for  the  momentum  and 
the  constitutive  relation  for  a layered  composite,  he  obtained  the  constitu- 
tive equations  in  binary  mixture  form  in  which  the  "interaction"  term  for 
the  mixture  theory  w.is  deduced. 

One  can  also  compare  various  theories  by  looking  at  the  frequency 
equations  of  each  theory.  It  shouUI  be  noted,  however,  that  if  a theory 
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is  a t'lrst  orilor  thoor>  , oik-  oati  only  look  at  terms  win  oh  arc  of  first 
Ollier.  It  does  not  make  sense  to  incluile  the  higher  order  terms  for  com- 
parison when  the  theories  are  of  first  order.  I'or  instance,  for  liarmonic 
waves  propa.cat  in>;  normal  to  the  layering  of  a com]iosite,  the  first  ovdi  i 
freipiency  equation  based  on  the  theory  of  interact  in.y  coin  inua  can  he 
obtained  from  Iq.  (J(>)  on  p.  IJ  as 


H -a_,e’k-)  - il  - b^t:^k\-Mo‘^  = 0 iSq) 

where  c^^,  k denote  non-dimensional  phase  velocity  and  wave  number. 
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defined  in  liqs.  iJ.'i)  and  l2") . 
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a function 

of  c'k' . If  we  assume  that 
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(OO) 

it  can  bo  sliown  that,  by  substituting;  1U\.  l,oOl  into  tS'.M  and  equat  iiij;  the 
coefficients  of  same  powers  in  ck. 


A = a , - b , 


where  use  has  been  made  of  I q.  i-”') . I'o  comiuire  l.q.  (^oOl  with  the 
equation  based  on  effective  stiffness  theory  (.Sl8|,  we  rewrite  I'q. 
[.SIS]  in  the  followiii);  form  b\  usin>;  the  non-dimensional  notations 
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n = 0 


In  Hq.  li>2),  the  matrix  layer  and  the  reinforcinj;  layer  are  denoted  by 
superscripts  lU  and  l2)  , respectively,  and  k is  defined  by  fq. 

Af;.iin,  a.ssiiminj;  that  c‘  is  jiiven  by  liq.  (oO) , one  obtains 


113 


A 


(63) 


Since  this  is  not  the  same  as  hq.  (61),  the  first  order  theory  of  interact- 
ing continua  and  that  of  effective  stiffness  are  not  equivalent.  However, 
one  could  make  the  two  theories  equivalent  if,  in  the  effective  stiffness 
theory,  a different  assumption  than  that  of  fS18]  is  made  on  the  smoothing 
operation  in  which  the  weighted  average  of  strain  energy  and  kinetic  energy 
were  calculated.  One  could  also  make  a different  assumption  on  the  dis- 
placement field  in  the  layers. 


A remark  on  higher  order  approximation  should  also  be  in  order.  For  the 
theory  of  interacting  continua,  the  frequency  equation  for  higher  order  appro- 
ximation is  obtained  by  adding  terms  of  (ek)*^,  (ek)^,..  in  l:q.(60)  without 
changing  tlie  coefficient  of  the  first  order  term  given  by  Eq.  (61).  For 
the  effective  stiffness  theory,  a second  order  approximation  studied  in  [D12] 
appeared  to  have  changed  the  coefficient  A given  by  Eq.  (63)  to  that  given 
by  Eq.  (61)  because  the  frequency  curve  of  the  second  order  appeared  to  have 
the  same  curvature  as  that  of  the  exact  solution  at  k = 0. 

In  summary,  it  is  seen  that  the  differences  in  the  theories  are  (a)  the 
assumptions  made,  (b)  the  ease  in  improving  the  accuracy,  (c)  the  ease  in 
solving  transient  problems  and  (d)  extensions  to  nonlinear  and  non-periodic 
composites.  The  theory  of  interacting  continua  seems  to  be  the  best  among 
all  theories  when  all  merits  and  drawbacks  are  taken  into  consideration. 

This,  of  course,  does  not  mean  that  one  should  use  the  theory  of  interacting 
continua  for  all  problems. 


V . 2 Recommendations 
V.2.1  lAjwar  composites 

For  the  response  of  composites  whose  governing  differential  equations 
arc  linear,  the  various  theories  presented  here  arc  adequate  for  predicting 
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fill'  first  iikhU'  t'r<,'i)iK'MO>’  I'niiat  ion . lor  the  transiont  proM  oms , t ho  tlioorios 
are  also  adequate  tor  predicting  the  lioad-of-t he-pul so  propagation  at  large 
time  and  at  largo  distance  from  the  point  of  impact.  Ihe  tlu'ories  are  in 
general  less  successful  in  predicting,  say,  the  space-wise  stress  variation 
at  any  ti.xed  time.  In  particular,  tl\e  stress  near  the  impact  end  and  at 
points  which  are  at  finite  distance  from  the  impact  end  has  not  been  i>redict- 
ed  satisfactorily  hy  any  of  the  existing  theories.  It  is  felt  that  a tiu'orv 
based  on  vi scv^elast ic  analogy  between  an  elastic  composite  and  a viscoelastic 
solid  may  be  able  to  accom|ilish  this  objective.  The  approach  may  predict 
not  only  the  short-time  response  but  also  the  long-time  liead-of-the-pulse 
response . 

V . J . J Non  1 i noa  r com^ios  i t os 

.\lthough  some  attempts  have  been  made  to  analyse  tlie  dviiamic  response 
of  nonlinear  composites,  tliere  seems  to  be  no  theory  avail.ible  wliich  can 
predict  the  transient  response  of  nonlinear  composites  sat isfactori 1> . 

Unlike  for  linear  composites,  imin'oving  the  accuracy  of  an  existing  theory 
for  nonlinear  composites  presents  anotlier  difficult  problem.  The  difficulty 
of  analyzing  the  transient  response  of  a nonlinear  composite  can  be  illu- 
strated by  considering  wave  propagation  normal  to  the  layering  of  a 
hi  Laminate  due  to  a step  norm.il  loail  on  the  surf.ace  of  the  comjHisite.  Tor 
this  simple  one-dimensional  problem,  the  waves  »n  the  first  layer  after 
the  apjilication  of  the  normal  load  will  be  simple  waves  or  a shock  wave 
ilepending  on  the  stress-strain  relation  of  the  layer.  Assuming  that  the 
waves  are  simple  waves,  the  reflected  waves  from  the  interface  boundar> 
between  the  first  aiul  the  second  laver  will  be  an  unloading  wave  which  may 
generate  a shock  wave.  The  waves  transmitted  to  the  seconti  layer  will  be 
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simple  waves.  As  one  can  see,  the  solution  is  alreajv  complicated  cnoijj;h 

I 

after  the  first  reflection  and  transmission.  For  a composite  which  contains 
many  layers  and  lienee  will  have  multiple  reflections  and  transmissions,  any 
attempt  to  solve  the  problem  c.\actly  or  nearly  exactly  i.s  not  likely  to 
Succeed . 

What  approaches  one  should  use  to  solve  nonlinear  composite  problems  ] 

] 

require  further  investigations.  It  seems  tliat  the  first  problem  one  should 
study  is  the  one-dimensional  waves  propagating  normal  to  the  layerings  of 
a nonlinear  clastic  bilaminatc.  Clearly,  the  shock  waves  within  each  layer 
have  to  be  smoothed  out  so  that  only  the  macroscoiiic  response  is  obtained. 

If  elastic-plastic  materials  arc  used  in  the  layerings,  one  lias  tlie  additional 
difficulty  of  tracing  the  unloading  and  reloading  boundaries.  As  in  the  case 
of  linear  composites,  a theory  based  on  nonlinear  viscoelastic  modeling  of 
nonlinear  composites  may  prove  to  be  tiic  most  jiractical  approach  for  solving 
transient  wave  problems  in  nonlinear  composites. 
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